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EDITORIAL 


En estos tiempos, la preocupación por los temas relacionados con los cambios climáticos ha 
ido en aumento entre la parte más informada de la sociedad. La difusión de problemas 
ambientales tiene un efecto positivo pues alienta la discusión sobre los efectos a largo plazo 
que estos fenómenos pueden tener sobre el ambiente y, sobre las actividades humanas. 
También, sirve para crear conciencia en buena parte de la población sobre las opciones que 
tenemos los seres humanos y acerca de cuáles medidas se deberían adoptar por la socie- 
dad a fin de atenuar los eventuales efectos negativos. Es de presumirse que, en el futuro 
previsible, nadie podrá argüir que no ha tenido a su alcance noticias acerca de lo que pasa 
en su entorno en relación con los temas y problemas ambientales; en ese escenario, presu- 
mimos que el conjunto de la sociedad tomará medidas informadas, con la aspiración de 
conservar las condiciones del ambiente lo menos alteradas posibles para permitir un desa- 
rrollo armónico del ser humano con la naturaleza. Esta revista tiene como propósito servir de 
foro para los estudiosos del tema forestal que, hoy por hoy, es uno de los que más relevancia 
tienen con respecto al aprovechamiento racional y la conservación de los recursos naturales. 

En este número incluimos seis trabajos de gran importancia. El primero tiene que ver el 
encalado del suelo y el potencial de inoculación de una multicepa con propágalos de hongos 
micorrícicos arbuscular en el desarrollo y acumulación de fósforo de Pinus halepensis Mili., 
cultivado en un suelo andisol. En invernadero se evaluó la aplicación de cal para elevar el pH 
del suelo a diferentes niveles con y sin la adición de micorriza. Los resultados mostraron que 
la inoculación micorrícica o el encalado del suelo para elevar el pH del suelo promovieron el 
mayor desarrollo de los brinzales medido como peso seco aéreo, diámetro del tallo, altura y 
volumen radical, así como en la acumulación de fósforo portas plantas, en comparación con 
el testigo (sin cal agrícola o sin inocular). El encalado estimuló el porcentaje de colonización, 
independientemente de la inoculación. La validación en campo del tratamiento de encalado 
e inoculación con hma promovió 95% de supervivencia de los árboles con incrementos sig- 
nificativos en altura y diámetro de tallo. Los resultados de invernadero y campo sugieren la 
conveniencia de encalar e inocular con micorriza arbuscular como prácticas efectivas para 
el establecimiento comercial de P. halepensis en andisoles del estado de México. 

El segundo artículo describe la patogenicidad de Chrysoporthe cubensis en clones de 
Eucalyptus grandis y E. urophylla en el sureste de México. Estas especies de eucalipto son 
las más promisorias para las plantaciones forestales en el sureste de México. Su potencial 
de producción está amenazado por el patógeno del cancro, Chrysoporthe cubensis (Bruner) 
Gryzenhout & M.J. Wingf, que afecta a la mayoría de los árboles y puede ocasionar daños 
al arbolado de E. grandis y afectar un gran porcentaje del volumen comercial en ambas 
especies. Este estudio se realizó para evaluarla patogenicidad de C. cubensis en clones de 
E. grandis y cuatro de E. urophylla. Los resultados mostraron que los aislamientos de 
C. cubensis son capaces de causar lesiones en las ramas jóvenes de los clones de ambas 
especies, observándose diferencias significativas entre los clones y los aislamientos. 

El tercer manuscrito trata el tema del crecimiento de árboles individuales de Geoffreea 
decorticans Burk. en un bosque del Chaco semiárido, Argentina. El objetivo del trabajo es 
contribuir al conocimiento del crecimiento de árboles individuales de esta especie (chañar). 
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definir en función de la edad la evolución del incremento radial (ir), diámetro a la altura del 
pecho (dap), sección normal (sn) y estimar el incremento medio (ima) y anual (ia) con el fin 
proponer directrices para la correcta gestión de la especie. El estudio fue realizado en un 
bosque del Chaco Semiárido, Parque Los Quebrachos, Alberdi, Santiago del Estero, Argen- 
tina. El chañar presenta anillos de crecimiento bien demarcados con incremento radial 
importante. Este crecimiento permite definir el manejo en rotaciones cortas. 

El siguiente documento plantea los estudios que se hicieron para generar la estimación de 
contenido de carbono de la biomasa aérea en la cuenca del Río Magdalena, México, de seis 
tipos de cobertura vegetal. Esta cuenca está ubicada en el suelo de conservación del Distrito 
Federal. Para ello, se elaboró un mapa de las clases de cobertura vegetal y uso de suelo del 
área de estudio mediante la fotointerpretación de fotografías aéreas. La cobertura vegetal que 
presentó mayor cantidad de carbono almacenado corresponde a la clase de Abies religiosa 
cerrado; le sigue la clase de Abies religiosa-Pinus hartwegii cerrado y la clase de Abies reli- 
giosa abierto. Las clases de Pinus sp. y Pinus hartwegii cerrado, de Pinus hartwegii-Abies 
religiosa cerrado y de Pinus hartwegii abierto, almacenan menor cantidad de carbono. 

Otro artículo describe los modelos desarrollados para estimar las tarifas de biomasa aérea 
para abedul (Betula pubescens Ehrh.) y roble fQuercus robur L.) en Galicia, noroeste de 
España. Se emplearon datos obtenidos con el muestreo destructivo de 50 abedules y 50 
robles, en los que se obtuvo el peso seco de biomasa total y por fracciones. Se seleccionaron 
los modelos que presentaban el mejor ajuste para cada fracción de biomasa arbórea conside- 
rada. Se realizó un ajuste simultáneo de las ecuaciones de estimación de biomasa por fraccio- 
nes junto con la de biomasa total, empleando el procedimiento estadístico denominado NSUR 
(Non linear Seemingly Un relatad Regressions). También fue necesario un ajuste ponderado 
para corregirla existencia de heterocedasticidad. El número de condición verificó que no exis- 
tían problemas graves de multicolinealidad. Al final se obtuvo, para cada especie, un sistema 
de siete ecuaciones de estimación de biomasa aérea para las distintas fracciones y para la 
biomasa total. Estas ecuaciones explicaron gran parte de la variabilidad observada en la bio- 
masa total para abedul y roble. Se recomienda la utilización de las ecuaciones desarrolladas 
en este estudio en sustitución de las ecuaciones de biomasa existentes para la región. 

Finalmente, se presenta un trabajo para hacerla reconstrucción dendrocronológica de patro- 
nes de precipitación de largo plazo en el Parque Nacional de Basaseachi, Chihuahua, 
México. El propósito fue generar una cronología de anillos de crecimiento y usarla como 
método indirecto para el desarrollo de una reconstrucción de precipitación estacional en el 
sitio de estudio. Para el desarrollo se colectaron núcleos de crecimiento de ejemplares de 
tres especies de pino. Con los índices dendrocronológicos y datos de precipitación de la 
estación climática más cercana se generó un modelo de regresión lineal bivariado que tuvo 
una varianza significativa. La precipitación de los meses de octubre, diciembre, febrero y 
abril fue la más correlacionada con el grosor de los anillos. La reconstrucción se extendió 
225 años, de 2007-1782. Durante este tiempo ocurrieron trece sequías severas, diez de las 
cuales afectaron grandes áreas en el norte y centro de México. Estamos seguros que estos 
trabajos serán del mayor interés para los estudiosos de los temas del aprovechamiento 
forestal y la conservación de los recursos forestales del mundo. 


Raymundo Dávalos Sotelo 

Editor 


EDITORIAL 


In these times, the concern for issues related to climate changa has been increasing between 
the more informad parí: of society. The wide spread dissemination of news on environmental 
problems has, in general, a positiva effect because it encourages discussion about the possible 
long-term effects these phenomena and facts can have on the environment and, ultimately, on 
human activities. At the same time, serves to raise awareness in good parí of the population 
about the options that we have as human beings and about which measures should be taken 
by society in order to mitígate any negativa effects. It is presumed that, in the foreseeable 
Mure, nobody can argüe that he/she have not had news about what is happening in the envi- 
ronment in relation to issues and environmental problems; in this scenario, we hypothesize that 
the whole society will take informad action, with the aim ofpreserving the ambient conditions as 
least alterad as possible, to permit harmonious development of human beings with natura. This 
Journal aims to próvida a forum for scholars of the forest subject which, today is one of those 
who have more relevance with regard to the rational use and conservation of natural resources. 

In this number, we inelude six aríicles of gear interest and relevance. The first one explains 
that, in order to evalúate soil liming and the potential of inoculation of a multistrain with pro- 
pagules of arbuscular mycorrhizal fungi (amf) in the development and accumulation of P. 
halepensis Mili phosphorus, grown in an Andisol, it was evaluated, in a greenhouse, the 
application of lime to raise the soil pH, with and without the addition of mycorrhizae. The 
results showed that mycorrhizal inoculation and the soil liming to pH of 6,0 to 6,5 led to a 
furíher development of seedling, measured in dry aerial weight, stem diameter, height, and 
root volume, as well as phosphorus accumulation byplants, com pared with the control (with 
no agricultural liming and without inoculation with amf). Liming stimulated the percentage of 
colonization, regardiess inoculation with amf. The fieid validation of the liming treatment and 
AMF inoculation promoted 95% survivai oftrees with an imporíant increase in height and stem 
diameter, during an evaluation period of five months. The greenhouse and fieid results 
suggest the benefits of liming and inoculation with arbuscular mycorrhizal as effective practi- 
cas for the commercial establishment of P. halepensis in Andisols in the State of México. 

The second aríicle deais with the pathogenicity ofChrysopoúhe cubensis to Eucalyptus gran- 
dis and E. urophylla clones in Southeastern of México. Eucalyptus grandis S.T Blake and 
E. urophylla HUI ex Maiden are the most promissory species for forest plantations in the 
southeastern of México. Their production potential is threatened by the canker pathogen, 
Chrysoporthe cubensis (Bruner) Gryzenhout & M.J. Wingf, which affeets most ofthe traes 
and can cause damage of about 43,6% in E. grandis trees and up to 40%> ofthe merchanta- 
ble volume of both species. This study was carried out to assess the pathogenicity of 
C. cubensis in six E. grandis and four E. urophylla clones. Results showed that the fouriso- 
lates ofC. cubensis are capable ofeausing lesions on young twigs of Eucalyptus clones and 
significant differences were observad between clones and isolates. All testad clones were 
susceptible to isolateA3, but exhibited varying leveis oftolerance to the attack. 

A third papar studies the growth of individual trees of Geoffrea decorticans Burk. In a forest from 
the semiarid Chaco, Argentina. The objective ofthis study i/i/as to improve the knowledge on the 
annual growth retes of Geoffreea decorticans (chañar), to define the evolution ofannual growth 
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increment (ir), diameter at breast height (dar), normal section (dn) with age, to estímate mean 
annual (ima) and annual Increment (ia) In orderto propose guldellnes for proper management of 
the specles. The study i/i/as carried out In a semiarid Chaco forest, at Quebrachos Park, AlberdI, 
Santiago del Estero Province, Argentina. Definitlon of ring boundarles and ring wlth measuring 
were made wlth computerizad equipment aniol and the computar program catras. The data 
were analyzed wlth the statistical software package infostat. Chañar presentad well-defined 
annual growth rings. The age ofcompletlon for dar was determinad to be the five years, and ns 
Is achieved after 37 years. This growth rate can define the management on short rotatlons. 

The following papar presente an estimatlon of carbón content In the Magdalena RIver basin, 
México. Carbón content of aboveground blomass ofslx types ofvegetatlon cover In the Magda- 
lena RIver basin, located In the conservatlon land In the México City i^'as estimated. Henee, a 
map of cover vegetatlon and land use to the study area i/i/as made; a hierarchieal stratlfied sam- 
pling design wlth distributlon systematic wlth a mesh ofsampling sites equidistant of 550 meters 
was used forfieid work. The vegetatlon cover that showed the highest amount of carbón stored 
corresponde to the forest of the f/?e Abies religiosa closed, followed by the Ab\es religiosa - Pinus 
hartwegii closed and Abies religiosa open; next were the covers of Pinus sp, Pinus hartwegii 
closed followed by Pinus hartwegii - Abies religiosa closed and Pinus hartwegii open. 

The next paper presenta modeis developed for above-ground blomass equatlons for birch 
(Betula pubescens Ehrh.) and pedunculate oak fQuercus robur L.) In Galicia (North Western 
Spain). Data from a destructiva sampling of 50 birch traes and 50 pedunculate oak trees, that 
provided total dry blomass weight and dry blomass weight by components: wood, bark, bran- 
chesofdifferentthicknesses, andieaves, was used. DIfferentmodelswerefitted for each above- 
ground component and the best ofthem was chosen In a first step. All the modeis were refitted 
simultaneously by use of Nonllnear Seemingly Unrelated Regressions (nsur) to ensure the 
additivlty of the blomass equatlons. A system ofseven equatlons for the different above-ground 
blomass components, and for the total blomass, i/i/as obtalned for each specles. These equa- 
tlons explalned a great dea! ofthe observad varlablilty for total blomass. The developed blomass 
equatlons are recommended In substitutlon ofthe existing regional blomass equatlons. 

Flnally a paper Is ineluded that describes the fendrochronologlcal reconstructlon ofiong-term 
precipitatlon patterns In BasaseachI National Park, Chihuahua, México. The purpose of this 
study was to reconstruct the precipitatlon at BasaseachI National Park (bnr) In Chihuahua, 
México. Tree-ring samples from pine specles including Pinus durangensis, P. lumholtzii and 
P.engelmannii were collected in and near BNP and and they were cross-dated with existing 
chronologies. Ring widths ofeach sample were measured and modeis to remove non-clima- 
tic trends from the data were applied. The relationship between precipitation from the weather 
station in Yécora, Sonora, and the ring width Indices were modeled using a bivariate linear 
regression. It was found that precipitation from the months of October, December, February 
and April was most highiy correlated with tree ring growth. The reconstructlon extends 225 
years from 2007-1782. During this time period there were thirteen severe droughts, ten of 
which affected large areas in northern and central México. We are sure these papers will be 
very useful for the researchers interested in forest utilization and management, as well as for 
those individuáis and organizations focused in work to preserve the worid's natural resources 


Raymundo Cávalos Sotelo 

Editor 
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ARTÍCULO DE INVESTIGACIÓN 


Encalado y micorriza para corregir deficiencia 
de fósforo en un Andisol cultivado 
con Pinus halepensis 

Lime and mycorrhizae in the phosphorus deficiencies 
correction in an andisol cultivated with Pinus halepensis Mili. 

María E. Álvarez-Sanchez\ Elizabeth Hernández-Acosta^*, 
Ranferi Maldonado-Torres^ y Miguel Rivera-González^ 

RESUMEN 

Con la finalidad de evaluar el encalado del suelo y el potencial de inoculación de una multicepa con 
propágulos de hongos micorrícicos arbuscular (hma) en el desarrollo y acumulación de fósforo de Pinus 
halepensis Mili., cultivado en un Andisol, en invernadero se evaluó la aplicación de cal para elevar el pH 
del suelo de 5,3 a 5,5; 6,0 y 6,5; con y sin la adición de micorriza. Los resultados mostraron que la inocu- 
lación micorrícica o el encalado del suelo para alcanzar un pH de 6,0 a 6,5; promovieron el mayor desa- 
rrollo de los brinzales medido como peso seco aéreo, diámetro del tallo, altura y volumen radical, así 
como en la acumulación de fósforo por las plantas, en comparación con el testigo (sin cal agrícola o sin 
inocular con hma). El encalado estimuló el porcentaje de colonización, independientemente de la inocula- 
ción con HMA. La validación en campo del tratamiento de encalado a pH de 6,5 e inoculación con hma 
promovió 95% de supervivencia de los árboles con un incremento promedio en altura y diámetro de tallo 
de 22,1 cm y 5 cm, respectivamente, durante un período de evaluación de cinco meses. Los resultados 
de invernadero y campo sugieren la conveniencia de encalar e inocular con micorriza arbuscular como 
prácticas efectivas para el establecimiento comercial de P. halepensis en Andisoles del estado de México. 

PALABRAS CLAVE: 

Carbonato de calcio, inoculación micorrícica arbuscular, Pinus halepensis Mili. 


ABSTRACT 

In order to evalúate soil liming and the potential of inoculation of a multistrain with propagules of 
arbuscular mycorrhizal fungi (amf) in the development and accumulation of P. halepensis Mili phospho- 
rus, grown in an Andisol, it was evaluated, in a greenhouse, the application of lime to raise the soil pH 
from 5,3 to 5,5; 6,0 and 6,5; with and without the addition of mycorrhizae. The results showed that 
mycorrhizal inoculation and the soil liming to pH of 6,0 to 6,5; led to a further development of seedling, 
measured in dry aerial weight, stem diameter, height, and root volume, as well as phosphorus accumu- 
lation by plants, compared with the control (with no agricultural liming and without inoculation with amf). 
Liming stimulated the percentage of colonization, regardiess inoculation with amf. The fieid validation of 
the liming treatment to pH 6,5 and amf inoculation prometed 95% survival of trees with an average 
increase in height and stem diameter of 22,1 cm and 5 cm, respectively, during an evaluation period of 
five months. The greenhouse and fieid results suggest the benefits of liming and inoculation with arbus- 


1 Posgrado en Agroforestería para el Desarrollo Sostenible. Universidad Autónoma Chapingo. km 38.5 carre- 
tera México-Texcoco, Texoco, Méx. 56230. C.e.: edna_alvarez30@yahoo.com. mx 

2 Exestudiante de posgrado. Universidad Autónoma Chapingo. km 38.5 carretera México-Texcoco, Texoco, 
Méx. 56230 

* Autor para correspondencia: elizahac@yahoo.com 
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Encalado y micorriza para corregir deficiencia de fósforo en un Andisol 


cular mycorrhizal as effective practicas 
for the commercial establishment of P. 
halepensis in Andisols in the State of 
México. 

KEY WORDS: 

lime, arbuscular mhycorrizae, Pinus hale- 
pensis Mili, Andisoil, phosphorus. 


INTRODUCCIÓN 

En el oriente del estado de México existen 
bosques de clima templado con predomi- 
nio de especies de cipreses, encinos y 
pinos. Los suelos donde se desarrolla 
esta vegetación son principalmente Andi- 
soles, cuyas características distintivas son 
pH ácido y alto contenido de alófano, 
material amorfo que retiene el fósforo dis- 
ponible y lo hace muy poco aprovechable 
para las plantas (Shoji et al., 1993). La 
baja disponibilidad de fósforo en estos 
suelos y el exceso de aluminio activo limi- 
tan el desarrollo y la producción de las 
plantaciones forestales. La reforestación y 
el aprovechamiento comercial de R hale- 
pensis en Andisoles del estado de México 
como alternativa de manejo en estos eco- 
sistemas alterados, requiere en el corto 
plazo después del trasplante, una rápida 
adaptación que maximice la supervivencia 
de los árboles: para esto, es necesario 
generar condiciones óptimas de nutrición 
para el desarrollo de dicha especie fores- 
tal, entre otros factores. El encalado de los 
suelos es una de las prácticas que se 
emplean para corregir deficiencias de fós- 
foro, reducir los niveles de toxicidad del 
aluminio y favorecer la disponibilidad de 
nutrimentos como el calcio, magnesio y 
molibdeno (Alvarado y Fallas, 2004); tam- 
bién se recurre a la inoculación con mico- 
rrizas para aumentar la eficiencia en el 
aprovechamiento de los nutrimentos del 
suelo, en especial de fósforo. 

Hasta hace pocos años, se conside- 
raba que especies de la familia Pinaceae 
sólo desarrollaban simbiosis con los hon- 


gos ectomicorrícicos, sin embargo, traba- 
jos realizados por Smith et al. (1998) y 
Cameron et al. (2011) demuestran que 
especies de la familia Pinaceae, como 
P contorta, P ponderosa y P strobus, for- 
man simbiosis con hongo micorrícico 
arbuscular (hma) de la especie Glomus 
intraradices. La inoculación con hma en 
algunas especies maderables en vivero 
ha mostrado algunos resultados promete- 
dores en la supervivencia de los árboles 
en programas de restauración de comuni- 
dades tropicales (Ramos-Zapata y Gua- 
darrama, 2004) y paisajismo (Appleton et 
al., 2003; Wiseman y Wells, 2005). 

La importancia ecológica de la colo- 
nización con HMA en especies vegetales 
típicamente ectomicorrícicos es escasa y 
amerita determinar si esta relación favo- 
rece la nutrición, desarrollo y supervivencia 
de plántulas de pinos para su producción 
comercial. 


OBJETIVO 

Evaluar el encalado de un Andisol y el 
potencial de la inoculación con micorriza 
arbuscular en el desarrollo y acumulación 
de fósforo de brinzales de pino P halepensis 
Mili, y validar en campo el mejor tratamiento 
en una plantación con dicha especie. 

METODOLOGÍA 

Experimento en invernadero 

Se empleó un suelo con pH moderada- 
mente ácido (5,3) medido en agua, un 
contenido bajo de fósforo disponible 
determinado por el método P-Bray (5 mg 
kg-^), según las metodologías indicadas 
en la Norma Oficial Mexicana nom- 021 -rec- 
nat-2000 (Se-marnat, 2002). Con base en 
sus características, este suelo es clasifi- 
cado como Andisol. 
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Se probaron dosis de 0 g kg-^ ; 0,2 g kg-^ ; 
1 ,0 g kg-^ y 2,0 g kg-^ de suelo de cal agrícola 
(CaCOg) con 90% de pureza, para incre- 
mentar el pH de 5,3 a 5,5; 6,0 y 6,5; con y 
sin inoculación micorrícica. La dosis de cal 
aplicada en cada tratamiento se estimó de 
acuerdo con la técnica de la curva de titu- 
lación descrita por Goijberg y Aguilar 
(1987), para alcanzar el pH deseado. El 
inóculo micorrícico comercial conocido 
como Endospor formado por las especies 
Gigaspora margarita, Glomus mosseae, 
Glomus clarum, Glomus deserticola, Glo- 
mus etunicatum, Glomus brasilianum y 
Glomus intraradices, se aplicó a una dosis 
de 5 g por brinzal, a través del agua de 
riego en el momento de la plantación. 

El experimento correspondió a un 
factorial con el primer factor dado por el 
encalado (con cuatro niveles) y el segundo 
factor inoculación micorrícica con dos 
niveles (sin y con inoculación); su combi- 
nación generó ocho tratamientos. Dadas 
las condiciones de invernadero, se empleó 
un diseño experimental completamente al 
azar con cinco repeticiones. La unidad 
experimental estuvo constituida por un 
recipiente que contenía un brinzal. 

La dosis de cal se aplicó al suelo 15 
días antes del trasplante, tiempo sufi- 
ciente para lograr que se estabilizara el 
pH al valor propuesto; posteriormente, los 
brinzales de R halepensis Mili, de 41 cm 
(desviación estándar {ds) = ± 3,0 cm) de 
altura promedio se colocaron en bolsas de 
polietileno negro con capacidad para 5 kg 
de suelo. El experimento tuvo una dura- 
ción de cinco meses. 


Medición de variabies de crecimiento 

Al final del experimento, en cada uno de 
los brinzales, se midió la altura total de 
planta desde la base hasta el ápice y el 
diámetro del tallo en la base. Una vez 
extraído y lavado el sistema radical, se 


determinó su volumen en una probeta a 
partir del nivel de agua desplazada. Se 
colectaron muestras de raíz para evaluar 
la colonización micorrícica; el resto del 
material vegetal se secó en estufa de aire 
forzado a 70 °C y se registró el peso seco, 
luego se midió la concentración de fósforo 
en el tejido de la parte aérea en un dige- 
rido (mezcla de HCIO^ y HNOg concen- 
trado, relación 1:2), mediante el método 
del fosfovanadomolibdato (Jackson, 1970). 
También se tomaron muestras de suelo de 
cada unidad experimental, para determi- 
nar pH (en agua, relación 1:2) y P-Bray 
(Semarnap, 2002). 

Para evaluar la colonización micorrí- 
cica, se tiñeron raíces finas por el método 
de Phillips y Hayman (1970); el porcentaje 
de longitud radical colonizada se deter- 
minó por el método de las intersecciones. 
También se realizaron cortes histológicos 
y tinción en parafina para observar la pre- 
sencia de las estructuras micorrícicas 
(hifas, esporas y arbúsculos). 

Las variables de respuesta evaluadas 
se sometieron a un análisis de varianza 
acorde con el diseño de tratamientos y 
experimental empleados (factorial comple- 
tamente al azar) y a una prueba de medias 
entre niveles de encalado e inoculación 
micorrícica según el procedimiento de 
Tukey (a=0,05) para las variables diámetro 
de tallo, peso seco aéreo, fósforo acumu- 
lado en planta, volumen de raíz y segmen- 
tos colonizados, en donde la interacción no 
resultó significativa. 


Experimento en campo 

El experimento de campo consistió en 
evaluar el mejor tratamiento obtenido en 
el invernadero, que fue adición de cal 
agrícola para alcanzar un pH de 6,5 e ino- 
culación micorrícica. Éste se instaló en el 
predio forestal El Llano, ubicado en el 
ejido San Miguel Coatlinchán, Texcoco, 
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México, en las coordenadas geográficas 
19° 24’ 20” y 19° 29’ 30” latitud norte y 
98° 47’ 00” y 98°29’ 30” longitud oeste. El 
suelo del predio presentó a una profundi- 
dad de 0 cm a 30 cm, pH de 5,3; 6,0% de 
materia orgánica; 5 mg kg-^ de P-Bray-1 ; 
6,2 cmol kg-^; 1,8 cmol kg-^ y 0,9 cmol kg-^ 
de Ca, Mg y K intercambiables; 31 ,0 mg kg \ 
1,0 mg kg \ 16,2 mg kg-^ y 0,5 mg kg-^ de 
Fe, Cu, Mn, Zn y B; y 42% de saturación 
de aluminio, según las metodologías indi- 
cadas en la Norma Oficial Mexicana 
NOM-021-RECNAT-2000 (Semamap, 2002). 

El experimento de campo consistió 
en la siembra de 40 árboles de R halepen- 
sis, distribuidos a un espaciamiento pro- 
medio entre árboles de 2,5 m por 1 ,8 m. El 
encalado se realizó 60 días antes del tras- 
plante, las plantas se inocularon con la mico- 
rriza y se procedió a la plantación de 20 
árboles con altura inicial promedio de 41 ,8 
cm (desviación estándar {ds) = ± 3,5 cm) y 
un diámetro de tallo de 5,5 cm {ds = ± 0,8 
cm), el resto de la población correspondió 
al tratamiento testigo (sin cal y sin inocula- 
ción con micorriza) que presentaba una 
altura de 49 cm {ds = ± 5,6 cm) y un diá- 


metro de 5,4 cm {ds = ± 0,7 cm). Se eva- 
luó altura de planta y diámetro del tallo, al 
inicio de la plantación y después de cinco 
meses de su establecimiento en campo, 
también se cuantificó la supervivencia de 
brinzales. 


RESULTADOS Y DISCUSIÓN 

Experimento de invernadero 

El encalado afectó positiva y significativa- 
mente el desarrollo de los brinzales en las 
variables altura, diámetro del tallo y volu- 
men de raíz; la acumulación de fósforo en 
la planta y el porcentaje de segmentos de 
raíz colonizados. La inoculación con hma 
influenció significativamente a todas las 
variables de crecimiento; P acumulado en 
la planta y segmentos de raíz colonizados. 
La interacción de estos factores sólo fue 
significativamente positiva en la altura de 
planta. Por otra parte, tanto el encalado 
como la inoculación con hma no influyeron 
en la disponibilidad de fósforo (P-Bray) en 
el suelo (Tabla 1). 


Tabla 1. Pruebas de significancia de los factores niveles de encalado (E), micorriza (hma) 
y su interacción en las variables de crecimiento, fósforo acumulado y segmentos de raíz 

colonizados en brinzales de pino. 


Fuentes de 
variación 

Aitura 

Diámetro taiio 

Peso 

seco aéreo 

Voiumen de 
raíz 

P 

acumuiado 

Segmentos 
de raíz colonizados 

P-Bray 

Encalado 

** 

- cm 

** 

g planta"' 
ns 

cm'^ 

** 

mg planta"' 
** 

% 

** 

mg kg"' 
ns 

HMA 

** 

** 

** 

** 

** 

* 

ns 

Ex HMA 

** 

ns 

ns 

ns 

ns 

ns 

ns 

^CV (%) 

6,2 

8,5 

10,7 

14,9 

18,8 

12,5 

20,9 


ns no significativo, ** significativo a P < 0,01 


significativo a P < 0,05 ^CV coeficiente de variación 
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En la figura 1 se muestra la interac- 
ción de los factores estudiados en altura 
de planta. Los niveles de encalado que 
elevaron el pH del suelo de 5,3 a 6,0 y 6,5 
(niveles 2 y 3, respectivamente), favore- 
cieron significativamente el crecimiento en 
altura de los brinzales con respecto a los 
tratamientos sin encalar y con o sin inocu- 
lar con hma; aunque sin diferencias signifi- 
cativas entre estos niveles. Con el nivel de 
encalado 3 (equivalente a pH 6,5) e inocu- 
lación con HMA, el incremento en altura fue 
de 48% con respecto al tratamiento sin 
encalar y sin inocular, y de 30,1 % con ino- 
culación micorrícica. En la tabla 2 se pre- 
senta el efecto del encalado y el de la 
inoculación micorrícica en las variables 
donde no se manifestó la interacción de 
estos factores. El aumento en el pH del 
suelo con el nivel 3 de encalado también 
favoreció significativamente el desarrollo 
de los árboles en términos de diámetro y 
volumen de raíz, con incrementos de 48,6 % 
y 53,8 %, respectivamente. Asimismo, la 
adición de micorrizas contribuyó en el cre- 
cimiento de los brinzales, aunque sólo 
resultó en un incremento de 9,3% en el 


diámetro del tallo y de 38,5 % en el volu- 
men radical. A diferencia de estas varia- 
bles de crecimiento, el peso seco de la 
planta sólo fue influenciado con la inocula- 
ción HMA, con un incremento significativo 
de 36,2%, con respecto al promedio de los 
tratamientos sin inocular. 

El encalado y la inoculación con hma 
no ejercieron un efecto significativo en el 
fósforo disponible en el suelo (P-Bray), 
éste se incrementó de 5 mg kg-^ a 5,9 mg 
kg \ contenido considerado como defi- 
ciente en el suelo (Castellanos, 2000). 
Sin embargo, tanto la adición de micorri- 
zas como el encalado del suelo se refle- 
jaron en la nutrición con fósforo de los 
brinzales. En la tabla 2 se muestra que 
tanto el aumento en el pH del suelo a 6 
(nivel 2 de encalado) como la adición de 
micorriza favorecieron el fósforo acumu- 
lado en la planta, con un incremento de 
55,9% y de 31,1%, en comparación con 
el promedio de los tratamientos sin enca- 
lar o sin inocular, respectivamente. Los 
niveles de fósforo extractable no necesa- 
riamente reflejan la disponibilidad de fós- 




Nivel de encalado 


Figura 1 . Efecto de la interacción nivel de encalado e inoculación micorrícica 
en la altura de planta de brinzales de pino {R halepensis Mili.). 
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foro para las plantas, ya que las hifas de 
los hongos pueden mejorar el acceso a 
este nutrimento por las raíces, debido a 
su gran capacidad de explorar entre los 
poros más finos del suelo, donde las raí- 
ces de la planta no pueden penetrar 
(Smith y Read, 2008); de incrementar la 
actividad de la rizósfera a través de la 
liberación de fosfatases que mineralizan 
el fósforo orgánico (Feng et al., 2002) y 
de mantener tasas de influjo de fosfato 
hacia la superficie radical superiores al 
de las raíces de la planta (Smith et al., 
2001). Estos efectos pudieron ocurrir en 
los brinzales de P. halepensis, ya que fue 
notable el incremento en el volumen radi- 
cal generado por la micorriza (38,5%), 
con respecto al tratamiento sin encalar o 
sin inocular, lo que demuestra el efecto 
benéfico de aumentar la densidad radical 
y con ello la capacidad de exploración de 
las raíces en el suelo y, por tanto, de 
acceso a los nutrimentos. 

Si bien no hubo interacción entre los 
factores estudiados en la mayoría de las 
variables evaluadas, la inoculación con 
HMA ejerció un efecto sinérgico favorable 
para el desarrollo de P. halepensis. 

Existen pocos estudios del efecto de 
la inoculación con hma en pináceas o espe- 
cies maderables que son típicamente 
ectomicorrícicas, las respuestas en el cre- 
cimiento de las plantas pueden ir desde sin 
efecto (Smith et al., 1998; Appleton et al., 
2003), con algún efecto favorable Morrison 
et al., 1993) o incluso en detrimento del 
crecimiento (Cameron etal., 2011). 

Entre los factores que pueden influir 
en la falta de respuesta a la inoculación con 
HMA están: niveles de alta fertilidad (Morri- 
son et al., 1993) y de fósforo disponible en 
el suelo (Wiseman y Wells, 2005), el género 
dentro de la especie de micorriza empleada 
(Smith, 1998; Cameron, 2011) combinada 
con una adecuada fertilización de las plan- 
tas (Appleton etal., 2003). 


Colonización micorrícica 

Las observaciones realizadas al micros- 
copio permitieron identificar hifas y espo- 
ras, pero no arbúsculos. Estas estructuras 
se han identificado en raíces de pináceas 
inoculadas con especies de las familias 
Glomaceae (Smith et al., 1998; Cameron 
et al., 2011). En la tabla 2 se observa que 
las raíces del tratamiento sin encalar, así 
como las del tratamiento sin inocular estu- 
vieron altamente colonizadas por hma 
(31,1% y 46,9%, respectivamente). Dife- 
rentes investigaciones muestran que un 
suelo forestal (o sustratos compuestos 
con material vegetal) puede contener una 
importante proporción de propágulos de 
HMA (Morrison et al., 1993; Appleton et al., 
2003; Wiseman y Wells, 2005), dicho 
potencial de inoculación depende de un 
nivel deficiente en fósforo extractable 
(Montaño et al., 2001; Wiseman y Wells, 
2005) y que el suelo contenga fuentes 
naturales de propágulos derivados de 
plantas hospederas de hma, aunque la 
especie de planta hospedera puede, a su 
vez, definir el desarrollo de los hma (Carr- 
enho et al., 2002). 

El aumento del pH del suelo influyó 
significativamente en el desarrollo de la 
colonización por hifas y vesículas en las 
raíces de los brinzales de pino. A un pH de 
6,5 el porcentaje de colonización fue de 
70,2% y de 55,7% con la inoculación; en 
ambos casos, fue notable el incremento 
en el porcentaje de segmentos coloniza- 
dos con respecto al promedio de los trata- 
mientos sin encalar (31 ,1 %) o sin inocular 
(46,9%). La importancia de los resultados 
obtenidos en este estudio radican en que 
el incremento en los niveles de coloniza- 
ción alcanzados -aun con competencia 
de micorrizas autóctonas- estuvieron 
asociados de forma significativa en el 
corto tiempo (cinco meses) con el 
aumento en el crecimiento y nutrición con 
fósforo de P. halepensis. Diferentes inves- 
tigaciones realizadas en árboles madera- 
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Tabla 2. Efecto del encalado (medias sobre los tratamientos de inoculación micorrícica) 
y de la inoculación micorrícica (medias sobre los tratamientos de encalado) 

en las variables evaluadas. 


Factor 

Diámetro 

tallo 

Peso 

seco aéreo 

Volumen de 
raíz 

P 

acumulado 

Segmentos 
de raíz 
colonizados 

P-Bray 


cm 

g planta-1 

cm^ 

mg planta-1 

% 

mg kg-1 

Niveles de 
encalado 







0 

7,4c 

27,8a 

26c 

44,7 b 

31,1c 

5,4 

1 

7 Je 

29,5a 

25c 

57,1ab 

52,3b 

5,7 

2 

9,8b 

30,5a 

33,0b 

69,7a 

51,6b 

6,1 

3 

11,0a 

28,4a 

40,0a 

60,4a 

70,2a 

6,4 

"DMS1 

0,93 

3,78 

5,58 

13,22 

14,53 

ns 

Con HMA 

9,4a 

33,5a 

36,0a 

65,8a 

55,7a 

5,9 

Sin HMA 

8,6b 

24,6b 

26,0b 

50,2b 

46,9b 

5,8 

^DMS2 

0,49 

2,01 

2,97 

7,00 

7,40 

ns 


^DMS1 para comparar entre niveles de encalado, 

^DMS2 para comparar entre niveles de inoculación micorrícica, según prueba de Tukey (a=0,05). 

Letras iguales dentro de una columna según el factor de comparación, no son significativamente distintas. 


bles inoculados con hma, con frecuencia, 
no han mostrado beneficios tangibles en 
el crecimiento de los árboles, aunque 
hayan desarrollado un alto porcentaje de 
colonización; en otros casos, el beneficio 
de la inoculación se ha presentado años 
después. Por ejemplo, Morrison et al. 
(1993) inocularon con Glomus intraradi- 
ces o Glomus fasciculatum a las especies 
maderables Acer platanoide L., Sorbus 
aucu paria L., Malus Mili, y Fraxinus penn- 
sylvanica March.; la inoculación incre- 
mentó significativamente la colonización 
en dichas especies, pero sin efecto en el 
crecimiento. Después de la inoculación. 


las especies se trasplantaron a un 
segundo campo con competencia de hon- 
gos micorrícicos nativos; dos años des- 
pués del postrasplante, sólo Glomus 
intraradices mejoró en forma significativa 
el crecimiento de S. aucuparia. Appleton 
etal. (2003) utilizaron un producto comer- 
cial con micorrizas (los géneros predomi- 
nantes fueron Glomus clarum, Glomus 
etunicatum, Glomus intraradices y Entro- 
phospora columbiana) en encino {Quer- 
cus palustris y Quercus phellos) y maple 
{Acer rubrum), después de un año los 
árboles no mostraron respuesta a la ino- 
culación medida por el contenido de cloro- 
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fila o diámetro del tronco; sin embargo, 
cuando se hizo una combinación de inó- 
culo comercial con fertilizante (9% de N, 
45% P 2 O 5 y 1 5% de K 2 O) el peso seco de 
las raíces del encino y del maple se incre- 
mentó significativamente así como el por- 
centaje de colonización. Smith et al. 
(1998) inocularon plántulas de diferentes 
especies de pináceas con Glomus intrara- 
dices y, aunque las plántulas presentaron 
colonización después de ocho semanas 
de desarrollo, no hubo efectos significati- 
vos en el peso seco aéreo y de raíz. 

Por otra parte, estudios realizados 
por Cameron et al. (2011) indican que el 
efecto de la inoculación con hma en el 
crecimiento de plántulas de pino está 
relacionado con la familia a la que perte- 
necen los HMA. Ellos probaron dos espe- 
cies de HMA de las familias Glomaceae 
(G. intraradices y G mosseae), Acaulos- 
poraceae {Entrophospora columbiana y 
Acaulospora morrowiae) y Gigasporineae 
{Gig aspara margarita y Scutellospora 
calospora) en su capacidad para coloni- 
zar y afectar el crecimiento de pino 
blanco {Pinus strobus). 

Después de siete semanas de desa- 
rrollo, Glomus intraradices fue el único 
que colonizó las raíces de los pinos; las 
especies Acaulosporaceae redujeron la 
biomasa de las plántulas, mientras que las 
de Gigasporineae tuvieron el menor efecto 
sobre el peso seco. El producto comercial 
que se utilizó en la presente investigación 
contenía una multicepa compuesta por 
especies de la familia Glomaceae {Glo- 
mus mosseae, G. clarum, G. deserticola, 
G. etunicatum, G. brasilianum, G. intrara- 
dices) y una especie de Gigasporineae 
(Gigaspora margarita), por lo tanto, la 
compatibilidad de las especies de hma 
contenidas en el inóculo comercial, el 
nivel deficiente de fósforo extractable en 
el suelo y la práctica del encalado contri- 
buyeron a la respuesta significativa en el 
desarrollo y nutrición del R halepensis. 


Validación en campo 

La validación en campo, en un período de 
evaluación de cinco meses, con el trata- 
miento de encalado a pH de 6,5 y micorriza 
promovió 95% de supervivencia de los 
árboles, con un incremento promedio en 
altura de 22,1 cm (sd = ± 6,9 cm) y diámetro 
de tallo de 5 cm (sd = ± 1 ,2 cm), en compa- 
ración con el testigo, en el que sólo hubo 
1 0% de supervivencia. El incremento en la 
supervivencia de las plantas micorrizadas 
podría explicarse por una mejora del estado 
nutrimental estimulada por el encalado, en 
el momento de la plantación, en relación 
con las plantas no inoculadas. Los resulta- 
dos de este experimento en el período de 
evaluación (cinco meses) demuestran que 
la inoculación de los brinzales del pino en el 
momento del trasplante, mejoró la supervi- 
vencia y el crecimiento de esta especie en el 
Andisol donde se establecieron los pinos. 

La mayoría de los trabajos de investi- 
gación sobre reforestación con coniferas 
inoculadas con hongos ectomicorrícicos 
han mostrado, desde hace 60 años, que 
mejora los procesos de repoblación y reve- 
getación en distintas situaciones ambienta- 
les y de suelos (Pera y Parladé, 2005). 
Estos antecedentes no existen para conife- 
ras inoculadas con hongos endomicorríci- 
cos. Las interacciones de especies de 
hongos micorrícicos arbusculares con piná- 
ceas son recientes (Cameron etal., 2011). 


CONCLUSIONES 

Los resultados de invernadero y campo 
sugieren la conveniencia del encalado e 
inoculación con endomicorrizas, como 
una práctica con posibilidades de efectivi- 
dad para la producción comercial de R 
halepensis en la zona de estudio, debido 
al incremento significativo en altura, diá- 
metro de tallo, volumen radical coloniza- 
ción micorrícica, fósforo acumulado en los 
brinzales y de supervivencia en campo. 
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Patogenicidad de Chrysoporthe cubensis 
en clones de Eucalyptus grandis 
y E. urophylla en el sureste de México 

Pathogenicity of Chrysoporthe cubensis to Eucalyptus 
grandis and E. urophylla clones in the Southeastern of México 

J. Carlos Juárez-Palacios\ J. Amador Honorato-Salazar^, 

Luis Vázquez-Silva^ y J.F. Conrado Parraguirre-Lezama^ 

RESUMEN 

Eucalyptus urophylla S.T Blake y E. granáis Hill ex Maiden son las especies más promisorias 
para las plantaciones forestales en el Sureste de México. Su potencial de producción está amenazado 
por el patógeno del cancro, Chrysoporthe cubensis (Bruner) Gryzenhout & M.J. Wingf., el cual afecta 
a la mayoría de los árboles y puede ocasionar daños de 43,6 % en arbolado de E. granáis y hasta 40% 
del volumen comercial en ambas especies. Este estudio se realizó para evaluar la patogenicidad de C. 
cubensis en seis clones de E. granáis y cuatro de E. urophylla. Cuatro aislamientos de C. cubensis 
fueron inoculados en ramas jóvenes de los clones de ambos eucaliptos en condiciones de laboratorio 
y el aislamiento altamente patogénico fue seleccionado después de treinta días. Cuatro plantas de 
cada clon fueron subsecuentemente inoculadas en condiciones de invernadero con el aislamiento A3, 
que fue el más virulento, para determinar su tolerancia al patógeno después de trece meses. Los resul- 
tados mostraron que los cuatro aislamientos de C. cubensis son capaces de causar lesiones en las 
ramas jóvenes de los clones de ambas especies, observándose diferencias significativas entre los 
clones y los aislamientos. Todos los clones fueron susceptibles al aislamiento A3, pero mostraron nive- 
les variables de tolerancia al ataque. 

PALABRAS CLAVE: 

Clones, Chrysoporthe cubensis, Eucalyptus, tolerancia, patogenicidad. 


ABSTRACT 

Eucalyptus granáis S.T. Blake and E. urophylla Hill ex Maiden are the most promissory species 
for forest plantations in the southeastern of México. Their production potential is threaten by the canker 
pathogen, Chrysoporthe cubensis (Bruner) Gryzenhout & M.J. Wingf., which affects most of the trees 
and can cause damage of about 43,6% in E. granáis trees and up to 40% of the merchantable volume 
of both species. This study was carried out to assess the pathogenicity of C. cubensis in six E. granáis 
and four E. urophylla clones. Four isolates of C. cubensis were inoculated on young twigs of both 
eucalyptus clones in laboratory conditions and the highiy pathogenic isolate was selected after thirty 
days. Four plants of each clone were subsequently inoculated in greenhouse conditions with the isolate 
A3, which was the most virulent, to determine their tolerance to the pathogen after thirteen months. In 
both triáis, lesión length in the bark was measured. Results showed that the four isolates of C. cubensis 
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2 Autor para correspondencia. INIFAP, Campo Experimental San Martinito. CIR-Golfo Centro. Km. 56.5 Carretera 
Federal México-Puebla, San Martinito, Tlahuapán, Puebla. C. P. 74100. C.e.: honorato.amador@inifap.gob.mx. 
Benemérita Universidad Autónoma de Puebla Unidad Regional Tetela. Calle Universidad s/n. Barrio Benito 
Juárez, Tetela de Ocampo, Puebla. 
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are capable of causing lesiona on young 
twigs of Eucalyptus clones and significant 
differences were observad between clo- 
nes and isolates. All tested clones were 
susceptible to isolate A3, but exhibited 
varying leveis of tolerance to the attack. 

KEY WORDS: 

Clones, Chrysoporthe cubensis, Eucalyp- 
tus, tolerance, pathogenicity. 


INTRODUCCIÓN 

El género Eucalyptus (Myrtales: Myrta- 
ceae), conforma una amplia diversidad 
biológica, con un número de especies 
entre 51 3 y 700 (Chippendale, 1 988; Potts 
y Pederick, 2000). Diferentes especies de 
eucalipto se plantan en muchos lugares 
del mundo para la producción de madera 
sólida y pulpa para producción de papel, 
en varios países de Europa, África, Asia, 
América y en algunos otros lugares como 
las islas del pacifico, lo cual constituye un 
componente muy importante para la 
industria silvícola. En México, 202 000 ha 
han sido plantadas con diferentes espe- 
cies forestales (Conafor, 2010), de las 
cuales Eucalyptus urophylla S.T. Blake y 
Eucalyptus grandis Hill ex Maiden son las 
más promisorias en el sureste de México, 
principalmente en los estados de Chiapas, 
Oaxaca, Tabasco y Veracruz. La produc- 
ción de estas especies es para pasta celu- 
lósica para producción de papel, madera 
de aserrío, pisos de parquet, postes para 
comunicaciones y construcción, leña y 
carbón vegetal (Ceccon y Martínez, 1999). 

El hongo Chrysoporthe cubensis 
(Bruner) Gryzenhout & M. J. Wingf. [previa- 
mente catalogado como Cryphonectria 
cubensis (Bruner) C. S. Hodges (Gryzen- 
hout et al., 2004)] es causante de la enfer- 
medad conocida como cancro del eucalipto. 
Este patógeno también ataca a otras espe- 
cies de árboles en las áreas tropicales y 
sub-tropicales del mundo (Wingfieid, 2003). 
Los daños causados por el hongo pueden 


ser hasta de 43,6% en arbolado de E. gran- 
dis (Paredes et al., 2010) y mayor de 40% 
en el volumen de producción de madera 
(Venter et al., 2001), sobre todo en planta- 
ciones comerciales. 

En años recientes, numerosos aisla- 
mientos de C. cubensis han sido colecta- 
dos de las diferentes plantaciones de 
eucalipto del mundo, para evaluar su pato- 
genicidad tanto en árboles como en clones 
de eucaliptos (Myburg et al., 1999; Myburg 
etal., 2002; 2003; Gryzenhout etal., 2005). 

Generalmente la patogenicidad de un 
organismo está influenciada por varios fac- 
tores, como la temperatura, disponibilidad 
de nutrientes, las sequías, la defoliación y 
los daños mecánicos, así como aspectos 
genéticos, tanto del patógeno como del 
hospedero (Oíd et al., 1986; 1990; Oíd y 
Davison, 2000; Guimaráes etal., 2010). En 
algunos casos, Chrysoporthe, está aso- 
ciado con cancros muy grandes como res- 
puesta del hospedero (Wardiaw, 1999; 
Nakabonge, 2006). Los cancros que se 
presentan anualmente en los árboles de 
eucalipto causan daños severos en el cam- 
bium, que pueden causar la muerte de los 
árboles, debido a la reacción virulenta del 
patógeno (Yuan y Mohammed, 2000). En 
otras partes del mundo, como en Sudá- 
frica, los cancros se desarrollan en la base 
de los tallos muy superficialmente, mien- 
tras que en Sudamérica y el sureste de 
Asia, éstos crecen a diferentes alturas de 
los tallos y en las ramas (Conradie et al., 
1990), lo que sugiere que puede existir 
cierta resistencia o tolerancia al ataque del 
patógeno en muchos individuos. 

En México, se ha encontrado el 
hongo de C. cubensis fructificando en 
plantas de Sericea (Clidemia: Melastoma- 
taceae), estas plantas son comunes como 
malezas que se encuentran en las orillas 
de los caminos y en las proximidades de 
las plantaciones de eucalipto (Ferreira, 
1 999). También se han encontrado estruc- 
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turas reproductoras de Chrysoporthe en 
corteza de eucalipto en plantaciones de 
Chiapas y Oaxaca, que están asociadas a 
los cancros encontrados en la base o en 
otras partes de los tallos (Juárez et al., 
2006; Bernal etai, 2009). 

En México se han realizado pruebas 
de patogenicidad con algunas especies de 
eucalipto con el hongo clasificado anterior- 
mente como Cryphonectria cubensis 
(Cruz, 2000); sin embargo, no existen 
reportes de este tipo de pruebas para E. 
grandisy E. urophylla, que cubren la mayor 
superficie de plantaciones forestales en el 
sureste de México, por lo que el objetivo 
del presente trabajo fue probar diferentes 
aislamientos de Chrysoporthe cubensis en 
clones de eucalipto para evaluar la patoge- 
nicidad del hongo y la resistencia de algún 
genotipo de Eucalyptus grandis y E. uro- 
phylla seleccionado y clonado. 

METODOLOGÍA 


Selección de clones 

Los clones de E. grandis y E. urophylla, 
fueron obtenidos a partir de genotipos 
seleccionados fenotípicamente en dos 
ensayos de fuentes de semillas, ubicados 
en las localidades de Aguascalientes, 
Chiapas (17°28’ LN, 92°14’ LO) y de La 
Gasolinera, Oaxaca (17°20’ LN, 95°23’ 
LO) (Tabla 1). El material vegetativo fue 
seleccionado directamente en las planta- 
ciones, seleccionando árboles fenotípica- 
mente superiores y con sanidad probada 
visualmente. Los árboles seleccionados 
se derribaron dejando un tocón a una 
altura de 0,30 cm para esperar sus rebro- 
tes. Una vez que hubo rebrotes, de cada 
tocón se seleccionaron ramillas tiernas de 
aproximadamente 0,5 cm de diámetro, las 
cuales fueron protegidas con papel 
secante humedecido con agua comercial 
de garrafón y colocadas en una hielera 


para evitar su deshidratación. Éstas, se 
llevaron a las instalaciones de la empresa 
Forestaciones Operativas de México, S.A. 
de C.V., ubicada en Las Choapas, Ver., en 
donde se enraizaron partes apicales de 
las ramas tiernas de cada fuente traída de 
campo, formando así los clones de cada 
genotipo seleccionado. Una vez enraiza- 
dos y con un tamaño de altura de 50 cm a 
60 cm, los clones se trasladaron al inver- 
nadero del C.E. San Martinito, Puebla, del 
Instituto Nacional de Investigaciones 
Forestales, Agrícolas y Pecuarias (inifap), 
en donde se realizaron las pruebas de 
patogenicidad. El número total de genoti- 
pos fue de seis para E. grandis y cuatro 
para E. urophylla. 

Los experimentos se realizaron en el 
laboratorio (a 19 °C ± 1,5 °C) e inverna- 
dero (a 30 °C ± 2 °C) estableciéndose dos 
ensayos; en el primero se realizó la inocu- 
lación de cuatro aislamientos de C. cuben- 
sis en porciones de ramas jóvenes y en el 
segundo la inoculación del aislamiento A3 
en plantas clonadas. En cada ensayo se 
utilizaron cuatro plantas de clones prove- 
nientes de los genotipos seleccionados en 
campo, en un diseño experimental com- 
pletamente al azar. 


Origen y selección de los aislamientos 

En las plantaciones experimentales de 
E. grandis y E. urophylla en Chiapas y 
Oaxaca, se recolectaron muestras en 
árboles que presentaban corteza con can- 
cros típicos y que contenían cuerpos de 
reproducción del hongo. Las muestras del 
hongo fueron aisladas en extracto de 
Malta agar al 2% y fueron preservadas a 
4 °C, obteniéndose cuatro aislamientos Al 
(S,A,,RJ, A2 (S 3 A,RJ, A3 (S 3 ^ 3 R,) y A4 
(S 3 A.,gR 3 ) de C. cubensis (Tabla 2), a los 
cuales se les realizó una identificación 
morfológica y molecular (Juárez et al., 
2006; Santos etal., 2006). 
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Tabla 1 . Selección de clones de Eucalyptus urophylla y E. grandis de cinco fuentes 
establecidas en dos sitios experimentales (Aguascalientes, Chis, y La Gasolinera, Oax.) 

del sureste de México. 


C/on® 

Especie 

Fuente^ 

Plantación 

CL1.- 2211(t) 

E. urophylla 

U-Pantar 

Aguascalientes, Chiapas 

CL2.- 2212(s) 

E. urophylla 

U-Pantar 

Aguascalientes, Chiapas 

CL3.- 2221 (t) 

E. urophylla 

U-Egon 

Aguascalientes, Chiapas 

CL4.- 2222(s) 

E. urophylla 

U-Egon 

Aguascalientes, Chiapas 

CL5.- 2231 (t) 

E. grandis 

G-386 

Aguascalientes, Chiapas 

CL6.- 2232(s) 

E. grandis 

G-386 

Aguascalientes, Chiapas 

CL7.- 2201 (t) 

E. grandis 

G-2464 

Aguascalientes, Chiapas 

CL8.- 2202(s) 

E. grandis 

G-2464 

Aguascalientes, Chiapas 

CL9.- 3231 (t) 

E. grandis 

G-5091 

La Gasolinera, Oaxaca 

CL10.- 3232(s) 

E. grandis 

G-5091 

La Gasolinera, Oaxaca 


Clones obtenidos en plantaciones experimentales de Eucalyptus spp. (Tolerantes (t) y susceptibles (s) a 
Chrysoporthe cubensis) en el sureste de México. 

Las fuentes de la semilla de eucalipto fueron: U-Pantar y U-Egon de Indonesia, G-386 de Cali Colombia, 
G-2464 de WaterVal, Sudáfrica y G-5091 de Waterhoutboom, Sudáfrica. 


Previo a las pruebas de patogenicidad, se 
evaluó la agresividad de los aislamientos 
para lo cual se inocularon ramas de euca- 
liptos de 13 meses de edad, tomando 
como criterio de mayor agresividad del 
patógeno, el mayor tamaño del cancro 
desarrollado en las ramas inoculadas 
(Juárez et al., 2006). 


Inoculación de ramas (varetas) 

En condiciones estériles de laboratorio se 
inocularon cuatro ramas jóvenes de cada 
una de las plantas de eucalipto crecidas 
en invernadero. Las ramas tenían 1 cm de 
diámetro y 15 cm de largo que se inocula- 


ron con los cuatro aislamientos del hongo 
C. cubensis, previamente identificado 
(Fig. 1). Para la inoculación se hizo una 
perforación en la corteza de la rama con 
un sacabocado de 3 mm de diámetro en la 
parte central de la rama. La inoculación se 
realizó tomando un disco de micelio pre- 
viamente desarrollado con conidias del 
hongo de los aislamientos seleccionados. 
Las ramas inoculadas se depositaron en 
cajas esterilizadas de plástico de 30 cm x 
25 cm X 20 cm con una rejilla y cerrado 
hermético, en donde se colocó agua des- 
tilada, colocándose entre la rama y el 
líquido una rejilla de plástico para evitar el 
contacto con el agua. Las cajas se coloca- 
ron en una estufa incubadora a 28 °C. La 
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Tabla 2. Selección de aislamientos de Chrysoporthe cubensis, provenientes de los 
ensayos de fuentes de Eucalyptus del sureste de México. 


Aislamiento 

Hospedero 

Origen 

A1 .- S 2 A. 1 . 1 R 4 

E. urophylla 

Aguascalientes, Chiapas 

A2.- S 3 A,R, 

E. urophylla 

La Gasolinera, Oaxaca 

A3.- S 3 A.|gR 2 

E. grandis 

La Gasolinera, Oaxaca 

A4.- S 3 A^gR 3 

E. grandis 

La Gasolinera, Oaxaca 


S=sitio, A=árbol, R=repetición. 


medición de los cancros se llevó a cabo 
con un calibrador (vernier) digital marca 
“Mitutoyo” Modelo CDN-P12. Las medicio- 
nes se realizaron cada cinco días en un 
lapso de 30 días. 


Inoculación en plantas jóvenes 

Para cada ensayo, se prepararon cajas de 
Petri (de 9 cm de diámetro) con extracto 
de Malta Agar (MYEA; 2% la malta extra y 
2% agar), en donde se desarrollaron los 
aislamientos del hongo por 10 días a 28 
°C creciendo en ellas colonias de micelio 
de color amarillo naranja. 

En cuatro plantas de cada clon de 
eucalipto de 18 meses de edad y de 1,5 
cm a 2,0 cm de diámetro, se inoculó el 
aislamiento A3 del hongo C. cubensis, 
identificado como el que mostró mayor 
agresividad en las pruebas de las vare- 
tas. Las inoculaciones se realizaron en 
los tallos a una altura de 20 cm, perfo- 
rando su corteza con un sacabocado de 
3 mm de diámetro, tomando después un 
disco del mismo diámetro de micelio con 
conidias del aislamiento del hongo y 
colocándolo en la corteza perforada con 
un trozo de gasa estéril impregnada con 
agua destilada y fijándolo con cinta de 


Parafilm “M” para evitar la pérdida de 
humedad (Fig. 1). Las plantas se distribu- 
yeron al azar sobre mesas en el inverna- 
dero a una temperatura de 30 °C ± 2 °C y 
se regaron cada tres días. La infección 
del hongo y desarrollo del cancro se 
evaluaron cada 10 días midiendo la longi- 
tud del cancro. 


Análisis del desarrollo de cancros 

La caracterización del desarrollo de los 
cancros, ocasionados por la inoculación 
del hongo, se realizó con base en las 
variables descriptivas c y b'^ del modelo 
de Weibull, considerando también la seve- 
ridad inicial (YJ y final (Y^), así como el 
tiempo que tardó en desarrollarse el can- 
cro como enfermedad (T^). La función de 
distribución de Weibull (Pennypacker et 
al., 1980) se seleccionó previa explora- 
ción gráfica de un total de 40 curvas de 
progreso inducidas por los cuatro aisla- 
mientos durante el tiempo de su desarro- 
llo, lo cual permitió obtener un modelo 
flexible para describir las diversas formas 
de las curvas del progreso de crecimiento 
de los cancros. 

El modelo se expresa de la manera 
siguiente: 
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Y=^-e[-{{/bY]^>0 (1) 

donde: 

Y = proporción de la severidad del cancro 
e = base del logaritmo neperiano (2,1718) 
t = tiempo en días 

b = parámetro asociado en su forma inversa 
con la tasa de velocidad del desarrollo 
del cancro 

c = parámetro asociado a la forma de la 
curva de progreso del cancro. 

El parámetro b está inversamente relacio- 
nado con la tasa de infección aparente (r) 
y el c está relacionado con la forma de la 
función de densidad (dy/d,) de la curva de 
progreso del cancro. Los valores b y c se 


determinaron mediante el ajuste de regre- 
sión no lineal, usando el procedimiento 
NLiN del paquete estadístico sas (sas. Ins- 
tituto Inc., 1991). Después de obtenerse 
los valores de b, se obtuvo su valor inverso 
(b'^), para su uso en los análisis posterio- 
res. Las otras variables descriptivas fue- 
ron: el tiempo total de la duración del 
desarrollo del cancro en la rama en días 
(T,), medido desde que se inició el desa- 
rrollo del cancro (YJ hasta que éste ter- 
minó en sus mediciones. 

Finalmente para conocer la diferen- 
cia estadística entre clones, agresividad 
de aislamientos, susceptibilidad de clones 
y diferencia de tiempo del desarrollo de 
cancros, se realizó un análisis de varianza 



Figura 1. Inoculación de plantas de eucalipto de 18 meses de edad con los aislamientos 
obtenidos. Ay D) desinfección y perforación del tallo con sacabocados, B y E) inóculo del 
hongo colocado en la perforación del tallo, C y F) cubrimiento del inóculo con gasa 
húmeda y Parafilm e inoculación de aislamientos en ramas y plantas de eucalipto de 1,0 

cm de diámetro. 
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de factores completamente al azar 
mediante el modelo general siguiente 
(Martínez, 1996; Wapole etal., 2007): 

y¡iki = ^iiki (2) 

Siendo el crecimiento del cancro, el 
efecto de los factores: tiempo (i), aisla- 
miento (j) y clon (k), para el ensayo de 
ramas y clon (i), resistencia (j) y etapa de 
desarrollo (k), para el ensayo de plantas y 
□,^^/el error experimental. Posteriormente 
se realizó una comparación múltiple de 
medias utilizando la prueba de Tukey a un 
nivel de significancia de a=0,05 (SAS, Ins- 
titute Inc., 1991). 


RESULTADOS 


Severidad de la infección en ramas de 
plantas jóvenes 

El desarrollo de los cancros inició a los 
cinco días después de la inoculación, pre- 
sentando diversidad en la forma e intensi- 
dad de desarrollo de los cancros en los 10 
clones evaluados. En su mayoría se 
encontraron curvas de desarrollo de tipo 
exponencial, sigmoidal, simétricas y asi- 
métricas, mostrando diversos grados de 
intensidad de la enfermedad (tasas epide- 
miológicas): también se observaron varia- 
ciones en el inicio y la duración del 
desarrollo del cancro. El modelo de Wei- 
bull describió adecuadamente el desarro- 
llo de la enfermedad a pesar de la 
heterogeneidad de las curvas de progreso 
del daño causado. Así, para los diez clo- 
nes, los coeficientes de determinación 
(R2) estuvieron en un intervalo entre 0,90 
y 0,98, superando la mayoría el 90 % de la 
variación explicada (Tabla 3). Por otra 
parte, en el análisis de varianza, las varia- 
bles estudiadas (clon, aislamiento y 
tiempo de desarrollo de la lesión y su inte- 
racción), fueron altamente significativas 
(Tabla 4). 


La mejor definición de los 10 clones 
de eucalipto por nivel de resistencia, fue a 
una distancia euclidiana promedio de 0,7 
en los cuatro aislamientos evaluados y las 
curvas de progreso promedio de similitud 
de resistencia de cada grupo, se presen- 
tan en forma gráfica para cada aisla- 
miento. 

El aislamiento Al (S 2 A.,.,R 4 ) presentó 
una severidad inicial (YJ que fluctuó entre 
0,48 cm y 3,53 cm, la cual se observó a 
los cinco días después de la inoculación 
en la totalidad de los clones. La mayor 
severidad inicial se presentó en los clo- 
nes CL2 y CL4, con valores promedios 
respectivos de 3,53 cm y 3,36 cm. La 
severidad final mayor (Y^) a los 30 días se 
observó en el clon CL2, con un desarrollo 
de 16,04 cm; mientras que la menor fue 
en el CL9 con 9,34 cm (Fig. 2). Los valo- 
res del estimador del parámetro epide- 
miológico b'^ , de la velocidad de desarrollo 
del cancro, varió de 0,0025 a 0,012, el 
valor menor fue en el clon CL9 y el mayor 
fue en el clon CL1 . El estimador c, tuvo 
valores que fluctuaron de 0,883 a 1,899, 
presentándose el mayor valor en el clon 
CL1 y en menor en el clon CL2. Los valo- 
res de ajuste R^ en este aislamiento, 
variaron de 0,959 a 0,986 (Tabla 3). 

Las curvas de progreso del aisla- 
miento A2 (S 3 A 4 R 4 ), presentó una severi- 
dad inicial (YJ que fluctuó entre 0,28 cm 
a 4,02 cm, la cual se observó a los cinco 
días después de la inoculación en la tota- 
lidad de los clones. La mayor severidad 
inicial fue en el clon CL 8 y la menor en el 
clon CL3. La severidad final mayor (Y^) se 
observó a los 30 días en el clon CL4, con 
un desarrollo promedio de 16,36 cm, 
mientras que el clon CL3 presentó la 
menor severidad final con 3,49 cm (Fig. 
3). Los valores del estimador del paráme- 
tro epidemiológico b'^ fueron entre 0,0005 
y 0 , 01 01 ; con el menor valor para el clon 
CL5 y mayor para el clon CL 6 . 
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Tabla 3. Valores promedio de las variables epidemiológicas asociadas a las curvas de 
progreso de la enfermedad inducida por Chrysophorte cubensis en 10 clones de 

eucalipto. 


Aislamiento 

Clon 

c 

b-1 

Yo 

Yf 

R2 


CU 

1,90 

0,0120 

0,48 

13,41 

0,983 


CL2 

0,88 

0,0046 

3,53 

16,04 

0,963 


CL3 

1,75 

0,0113 

0,65 

13,97 

0,959 


CL4 

0,91 

0,0049 

3,36 

15,95 

0,962 

Al 

CL5 

1,19 

0,0055 

1,38 

11,07 

0,984 

CL6 

1,30 

0,0078 

1,47 

14,02 

0,986 


CL7 

1,43 

0,0080 

0,98 

12,07 

0,967 


CL8 

0,92 

0,0044 

2,98 

14,48 

0,975 


CL9 

0,90 

0,0025 

1,92 

9,34 

0,963 


CL10 

1,12 

0,0056 

1,81 

12,73 

0,966 


CU 

0,93 

0,0021 

1,43 

7,40 

0,943 


CL2 

1,21 

0,0079 

1,98 

16,05 

0,975 


CL3 

1,43 

0,0032 

0,28 

3,49 

0,956 


CL4 

0,91 

0,0050 

3,43 

16,36 

0,955 

A2 

CL5 

0,74 

0,0005 

1,11 

4,13 

0,974 

CL6 

1,46 

0,0101 

1,28 

16,13 

0,969 


CL7 

1,42 

0,0055 

0,59 

7,35 

0,961 


CL8 

0,74 

0,0027 

4,02 

14,40 

0,935 


CL9 

1,29 

0,0048 

0,82 

7,88 

0,966 


CL10 

1,21 

0,0072 

1,78 

14,46 

0,946 


CU 

1,75 

0,0112 

0,65 

13,87 

0,919 


CL2 

0,81 

0,0038 

3,98 

15,94 

0,926 


CL3 

1,75 

0,0116 

0,69 

14,58 

0,945 


CL4 

0,85 

0,0043 

3,80 

16,22 

0,956 

A3 

CL5 

1,12 

0,0056 

1,83 

12,77 

0,942 

CL6 

0,98 

0,0057 

3,09 

16,47 

0,925 


CL7 

1,26 

0,0078 

1,69 

14,96 

0,914 


CL8 

0,69 

0,0027 

5,13 

16,48 

0,939 


CL9 

0,97 

0,0042 

2,35 

12,62 

0,948 


CL10 

0,93 

0,0054 

3,41 

16,74 

0,926 


CU 

1,72 

0,0103 

0,60 

12,45 

0,906 


CL2 

0,87 

0,0040 

3,24 

14,52 

0,968 


CL3 

1,81 

0,0112 

0,55 

13,03 

0,939 


CL4 

0,83 

0,0037 

3,54 

14,79 

0,937 

A4 

CL5 

0,88 

0,0022 

1,87 

8,69 

0,906 

CL6 

1,02 

0,0058 

2,71 

15,62 

0,904 


CL7 

1,25 

0,0057 

1,14 

10,28 

0,913 


CL8 

0,76 

0,0031 

4,09 

15,10 

0,921 


CL9 

0,89 

0,0027 

2,13 

10,15 

0,954 


CL10 

1,05 

0,0063 

2,61 

15,90 

0,940 


c y b"'' Forma de la curva y el estimado de la tasa de infección aparente, respectivamente, estimados por el modelo 
de Weibull. 

Coeficiente de determinación. 

Yg Severidad inicial. 

Yf Severidad final. 

Los valores son el promedio de cuatro repeticiones. 
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Tabla 4. Análisis de varianza de las variables estudiadas: clon, aislamiento y tiempo de 
desarrollo de la lesión ocasionada por los aislamientos del hongo C. cubensis en ramas 

de eucalipto. 


Fuente 

GL 

Suma de 
Cuadrados 

Cuadrado 

Medio 

Valor de F 

Pr>F 

Modelo 

279 

28368,625 

101,680 

313,88 

<0,0001 

Ti 

6 

21574,519 

3595,753 

11099,70 

<0,0001 

Ai 

3 

745,640 

248,547 

767,24 

<0,0001 

CL 

9 

3167,189 

351,910 

1086,31 

<0,0001 

Ti*Ai*CL 

261 

2881 ,277 

11,039 

34,08 

<0,0001 

Error 

840 

272,118 

0,324 



Total Corregido 

1119 

28640,742 





R2 

C. V. (%) 

Raíz del CME 

Promedio 



0,9905 

8,88 

0,569 

6,41 




Tiempo (días) 


Figura 2. Curvas de desarrollo de cancros ocasionados por el aislamiento A1 del hongo 
Chrysoporthe cubensis. Los números de cada curva de progreso representan los 
nombres de los clones que tuvieron el mismo comportamiento de resistencia a la 

infección del hongo. 
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Tiempo (días) 


Figura 3. Curvas de desarrollo de cancros ocasionados por el aislamiento A2 del hongo 

Chrysoporthe cubensis. 


El estimador c, tuvo valores que variaron 
entre 0,74 y 1,46, que correspondieron a 
los clones CL5 y CL6, respectivamente. 
En este aislamiento, los valores de ajuste 
de R2 variaron de 0,943 a 0,975 (Tabla 3) 
e impactó a cinco grupos de resistencia. 

Las curvas de progreso del aisla- 
miento A3 (S 3 A.,gR 2 ), presentaron una 
severidad inicial (YJ que fluctuó entre 
0,65 cm y 5,1 3 cm la cual se observó a los 
cinco días después de la inoculación en la 
totalidad de los clones, en donde la mayor 
severidad se presentó en el clon CL8, con 
un promedio de 5,13 cm y la menor en el 
clon CU con 0,65 cm. La severidad final 
mayor (Y^) se manifestó a los 30 días en el 
clon CL10, con un desarrollo de 16,74 cm 
y la menor en el clon CL9 con 12,62 cm 
(Fig. 4). 

Los valores del estimador del pará- 
metro epidemiológico b'^ resultaron ser 
entre 0,0027 y 0,0116, siendo el menor 
valor para el clon CL8 y el mayor para el 
clon CL3. Los valores del estimador c fue- 
ron de 0,69 a 1,75, presentándose el 
mayor valor en el clon CL1 y el menor en 


el clon CL8. Los valores de ajuste de R^, 
para este aislamiento, fluctuaron de 0,914 
a 0,956 (Tabla 3). 

Las curvas de progreso del aisla- 
miento A4 (SgA^igRg), mostraroo una 
severidad inicial (YJ que fluctuó entre 
0,55 cm a 4,09 cm, observándose a los 
cinco días después de la inoculación en 
la totalidad de los clones. La mayor seve- 
ridad inicial se presentó en el clon CL8 
con un promedio de 4,09 cm y la menor 
severidad fue para el clon CL3 (Fig. 5). 
La mayor severidad final (Y^) se observó 
a los 30 días en el clon CL10 con un 
desarrollo de 15,90 cm, mientras que la 
menor severidad final fue en el clon CL5. 
Los valores del estimador epidemioló- 
gico b \ de la velocidad de desarrollo del 
cancro, fluctuaron de 0,0022 a 0,0112, 
siendo el valor menor para el clon CL5y 
el valor mayor para el clon CL3. Los valo- 
res de c fueron entre 0,76 y 1 ,81 , presen- 
tándose el mayor valor en el clon CL3 y 
en menor en el clon CL8. Los valores de 
ajuste de R^en este aislamiento, variaron 
de 0,904 a 0,968 y se conformaron cuatro 
grupos de resistencia. 
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Figura 4. Curvas de desarrollo de cancros ocasionados por el aislamiento A3 del hongo 

Chrysoporthe cubensis. 


El análisis de varianza de las ramas 
indicó que hubo diferencias significativas 
(p < 0,05) entre clones, aislamientos y 
tiempo de desarrollo de los cancros 
(Tabla 4). La comparación de medias con 
la prueba de Tukey indicó que los clones 
se agruparon en siete grupos. Los clones 


CL8 y CL4 presentaron el mayor desarro- 
llo de cancro, mientras que el clon CL3 
fue el menos afectado. Asimismo, los clo- 
nes considerados como susceptibles 
tuvieron mayor desarrollo de cancro que 
los clones considerados como tolerantes 
(Tabla 5). 



Figura 5. Curvas de desarrollo de cancros ocasionados por el aislamiento A4 del hongo 

Chrysoporthe cubensis. 
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Tabla 5. Comparación de medias de las variables evaluadas 
(clones, aislamientos y tiempo) en el desarrollo de cancro causado por el hongo 

C. cubensis. 


CLONES 

Promedio 

AISLAMIENTO 

Promedio 

TIEMPO 

(días) 

Promedio 


CL8 

8,57 

A 

A3 

7,55 

A 

30 

12,27 

A 

CL4 

8,55 

A 

A4 

6,48 

B 

25 

11,08 

B 

CL2 

8,23 

B 

A1 

6,35 

C 

20 

9,70 

C 

CL6 

7,46 

C 

A2 

5,24 

D 

15 

5,96 

D 

CL10 

7,43 

C 




10 

3,77 

E 

CL9 

5,09 

D 




5 

2,05 

F 

CL7 

4,93 

E 




1 

0,00 

G 

CU 

4,81 

E 







CL5 

4,60 

F 







CL3 

4,46 

G 








Variables con la misma letra no son significativamente diferentes Tukey (a =0,05). 


Con respecto a la comparación de medias 
de los aislamientos, los resultados indica- 
ron que el aislamiento A3 fue estadística- 
mente (p < 0,05) diferente a los demás 
aislamientos y mostró mayor agresividad, 
seguido de los aislamientos A4, A1 y A2. 
El aislamiento que mostró menor agresivi- 
dad fue el A2. 

Hubo diferencias significativas (p < 
0,05) durante el desarrollo de cancros, 
aumentando gradualmente con el tiempo 
y mostrando el mayor desarrollo a 30 días 
después de la inoculación. 


Severidad de la infección 
de los cancros en plantas jóvenes 

El desarrollo de los cancros inició a los 15 
días después de la inoculación, 10 días 
más que en las ramas en condiciones 
controladas, presentando diversidad en la 
forma e intensidad de desarrollo de los 
cancros en los 10 clones evaluados. La 
dimensión de los cancros en plantas fue 
variable, alcanzando dimensiones de 3 
cm a 16 cm, por lo que en algunos casos, 
las plantas formaron una especie de callo 
en la lesión (Fig. 6. K, L). 
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Figura 6. Dimensión de los cancros en ramas y plantas de eucalipto. (G y L) cancro 
definido en las ramas y tallo de la planta, (H y K) tallo con formación de cancro con 
bordes suculentos en ramas y tallo de eucalipto (I) formación de cancro con presencia 
de cuerpos de reproducción del hongo (picnidios), (J) formación de cancros en tallos de 

plantas en invernadero después de la inoculación. 



Figura 7. Desarrollo de cancros del hongo Chrisoporthe cubensis en tallos 
de diferentes clones de eucalipto tolerantes y susceptibles al hongo, inoculados 
con el aislamiento A3 a los 15 (inicial) y 395 días (final) después de la inoculación. 
Los valores en paréntesis son el incremento del desarrollo. 
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El desarrollo del cancro inicial más bajo 
se presentó en el clon CL8, con un creci- 
miento promedio de 3,58 cm y el más alto 
en el clon CL3 con una dimensión de 
6,48 cm; mientras que en el desarrollo 
final del cancro, el clon CL3 mostró un 
crecimiento menor de 6,65 cm y el mayor 
fue en el clon CL8 con 17,48 cm. El incre- 
mento en el desarrollo de los cancros 
entre el crecimiento inicial y final (des- 
pués de 395 días) fue menor en los clo- 
nes tolerantes que en los susceptibles. El 
incremento en crecimiento fue de 0,18 
cm a 3,20 cm para los clones tolerantes y 
de 5,10 cm a 13,90 cm para los clones 
susceptibles (Fig. 7). 


El análisis de varianza de las varia- 
bles estudiadas en las plantas (clon, gru- 
pos de resistencia, etapas de desarrollo y 
su interacción) indicó diferencias alta- 
mente significativas (p < 0,05, Tabla 6). El 
clon CL8 presentó el mayor desarrollo del 
cancro, mientras que el clon CL7 mostró 
el menor desarrollo. Al igual que en las 
ramas, las plantas de los clones consi- 
derados como susceptibles presentaron 
mayor desarrollo de cancro que las 
plantas de los clones considerados como 
tolerantes (Tabla 7). Diferencias estadísti- 
cas significativas (p < 0,05) se encontra- 
ron entre las etapas inicial (15 días) y final 
(395 días). 


Tabla 6. Análisis de varianza de las variables: clon (CL), resistencia (res) 
y etapa (eta), del desarrollo de la lesión ocasionada por la inoculación del aislamiento A3 
del hongo C. cubensis en plantas jóvenes de eucalipto. 


Fuente 

GL 

Suma de 
Cuadrados 

Cuadrado 

Medio 

Valor de F 

Pr>F 

Modelo 

19 

1130,35 

59,49 

30,72 

<0,0001 

CL 

9 

223,45 

24,83 

12,82 

<0,0001 

RES 

1 

153,46 

153,46 

79,23 

<0,0001 

ETA 

1 

590,26 

590,26 

304,75 

<0,0001 

CL*RES*ETA 

8 

163,18 

20,40 

10,53 

<0,0001 

Error 

60 

116,21 

1,94 



Total Corregido 

79 

1246,56 





R2 

C. V. (%) 

Raíz del CME 

Promedio 



0,9068 

19,11 

1,392 

7,28 
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Tabla 7. Comparación de medias del comportamiento de las variables evaluadas 
de la enfermedad cancro del eucalipto ocasionada por C. cubensis, 

con respecto a su desarrollo. 


CLONES 

Promedio 

RESISTENCIA 

Promedio 

ETAPA 

Promedio 

CL8 

10,53 

A 

Susceptible 

8,67 

A 

Final 

9,93 A 

CL10 

9,15 

B 

Tolerante 

5,90 

B 

Inicial 

4,50 B 


CL4 9,11 B 

CL6 7,30 C 

CL2 7,25 C 


CL3 6,60 C D 

CL5 6,24 C D 

CL9 5,81 D 

CU 5,45 D 

CL7 5,43 D 

Variables con la misma letra no son significativamente diferentes Tukey (a =0,05). 


DISCUSIÓN 

Los estudios de resistencia inducidos 
con diversos aislamientos fungosos de 
Chrysoporthe spp. en especies de euca- 
lipto {Eucalyptus spp.) no son abundan- 
tes (e.g., Paredes et al., 2010); sin 
embargo, los estudios reportados indi- 
can que el eucalipto muestra una sus- 
ceptibilidad relativa a la enfermedad del 
cancro causada por diferentes aisla- 
mientos y especies de Chrysoporthe] 
además muestran que la severidad de 
daño causado varía entre especies y 
aislamientos fungosos de la misma 
especie, siendo algunos aislamientos, 
de la misma especie, más patogénicos 
que otros (van Heerden y Wingfieid, 
2001, 2002; Gryzenhout et al., 2003; 


Roux et al., 2003; van Heerden et al., 
2005, Chen et al., 2010; Chungu et al., 
2010; Guimaráes et al., 2010). 

Al igual que en los estudios de pato- 
genicidad reportados, en este estudio se 
observó que C. cubensis es capaz de 
causar lesiones en clones de eucaliptos. 
Las lesiones causadas en las ramas y 
árboles jóvenes de eucalipto, en el 
tiempo evaluado, sugieren que la patoge- 
nicidad de C. cubensis puede ser capaz 
de causar serios daños, como ha suce- 
dido con otras especies del patógeno 
(Gryzenhout et al., 2003; Chungu et al., 
2010) y encontrado en plantaciones 
comerciales del sureste de México (Pare- 
des et al., 2010). 
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En el ensayo de pruebas de patogeni- 
cidad de aislamientos se produjeron lesio- 
nes en los 10 clones probados. La 
respuesta fue similar a la reportada por 
Chen et al. (2010) en donde diversos aisla- 
mientos de C. cubensis causaron daños de 
diferente magnitud a clones de Eucalyptus 
grandis. También Roux et al. (2003) y van 
Heerden et al. (2005) reportaron variación 
en las lesiones causadas por diferentes 
aislamientos de Cryphonectria cubensis en 
diversos clones de eucalipto, lo cual con- 
cuerda con los resultados obtenidos en el 
presente trabajo, con la diferencia de que 
en este estudio se utilizaron varetas para 
determinar la patogenicidad de los aisla- 
mientos y en los estudios reportados se 
utilizaron árboles de 6 a 36 meses de edad 
bajo condiciones de invernadero. 

La longitud del daño obtenido en las 
varetas por aislamiento fue en promedio de 
52,4 mm a 75,5 mm, lo cual está dentro de 
los resultados reportados por Chen et al. 
(2010), quienes obtuvieron, después de 
seis semanas, un promedio de 22 mm a 62 
mm para sus seis aislamientos. Sin 
embargo, el daño en longitud es mucho 
menor a lo reportado por Roux et al. (2003) 
para Cryphonectria cubensis, en donde se 
registró un daño promedio de 61,90 mm a 
156,75 mm para tres aislamientos después 
de cuatro semanas de la inoculación. 

El comportamiento variable del daño 
causado por el aislamiento más virulento 
entre los clones, es similar a los resulta- 
dos obtenidos por Chen et al. (2010) para 
Chrysoporthe cubensis, por Conradie et 
al. (1992), van Heerden y Wingfieid (2002) 
y van Heerden et al. (2005) para Crypho- 
nectria cubensis y Gryzenhout et al. 
(2003) para Cryphonectria eucalypti, aun- 
que la extensión de daño es distinto por 
las diferencias en los tipos de aislamien- 
tos y el tiempo de medición. 

En la mayoría de los estudios de 
patogenicidad, el enfoque analítico ha 


sido univariado, esto es, que utiliza un 
solo parámetro de tipo epidemiológico, en 
la caracterización de clones o genotipos 
por su nivel de resistencia en condiciones 
de campo. Este enfoque; sin embargo, 
para la evaluación extensiva de clones o 
procedencias, no permite una completa 
caracterización de la resistencia (Ander- 
son et al., 1 990). En este estudio se utilizó 
un enfoque multivariado para caracterizar 
los 10 clones de eucalipto por su nivel de 
resistencia a los cuatro aislamientos. Este 
enfoque tiene como inconveniente la eva- 
luación exhaustiva de la intensidad de la 
enfermedad para construir las curvas de 
progreso del daño y posteriormente esti- 
mar parámetros epidemiológicos asocia- 
dos a estas curvas para su uso con fines 
comparativos (Kranz, 1983). La integra- 
ción en el análisis de varios parámetros 
permite explicar aspectos distintos de este 
fenómeno, así permite determinar la 
posible existencia de la resistencia en par- 
tes de plantas jóvenes en algunos clones, 
mientras que la tasa de intensidad de la 
epidemia b’^ permite tipificar el nivel de 
resistencia cuantitativa. En este trabajo se 
decidió incluir, para fines comparativos, el 
análisis de conglomerados. Los paráme- 
tros de tasa epidemiológica b’^ y forma de 
la curva (c) fueron incluidas en estas 
variables por su importancia en la caracte- 
rización de la resistencia (Parleviiet y van 
Ommeerent, 1975). 

El estudio de inoculación de aisla- 
mientos de C. cubensis en ramas y árboles 
jóvenes proporcionan datos importantes 
para la selección de plantas tolerantes al 
cancro causado por este patógeno, antes 
de realizar la plantación y tienen implicacio- 
nes importantes para los programas de 
mejoramiento de especies de eucalipto 
que tengan como objetivo reducir el 
impacto de C. cubensis. Asimismo, la ino- 
culación artificial en invernadero es de utili- 
dad potencial para seleccionar las plantas 
en un tiempo relativamente corto con un 
número reducido de plantas que puede 
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extrapolarse a pruebas de campo (Conra- 
die et al., 1 992), por lo que se ha utilizado en 
pruebas de exploración para la selección 
de especies de eucaliptos resistentes y el 
manejo efectivo del cancro causado, no 
solo por C. cubensis, sino por otras espe- 
cies (van Heerden y Wingfieid, 2002; Gryz- 
enhout et al., 2003; Roux et al. 2003; van 
Heerden et al., 2005; Chen et al., 2010; 
Guimaráes et al., 2010). 


CONCLUSIONES 

El estudio de patogenicidad usando dife- 
rentes aislamientos de Chrysoporthe 
cubensis en diferentes clones de Eucalyp- 
tus grandis y E. urophylla del Sureste de 
México indicó diferentes grados de tole- 
rancia al patógeno entre los clones, tanto 
en ramas como en los árboles jóvenes. 

Los clones susceptibles CL8 (Euca- 
lyptus grandis) y CL4 (E. urophylla) de las 
fuentes G-2464 y U-Egon, respectiva- 
mente, fueron los más afectados por el 
ataque de Chrysoporthe cubensis en 
ramas y plantas jóvenes. 

El modelo de Weibull describió ade- 
cuadamente el desarrollo de la enferme- 
dad considerando la severidad inicial (YJ 
y final (Y^), la tasa de infección aparente 
(b’^), el valor de la forma de la curva (c) y 
el tiempo de desarrollo del cancro (T^). 

Los clones inoculados, tanto en las 
varetas como en los árboles jóvenes mos- 
traron diferencias en su nivel de tolerancia 
al daño causado por C. cubensis, siendo 
mayor el daño en los clones susceptibles 
que en los clones tolerantes. 
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Crecimiento de árboles individuales 
de Geoffreea decorticans Burk., en un bosque 
del Chaco semiárido, Argentina 

Growth of individual trees of Geoffreea decorticans Burk. 
in a forest from the Semiarid Chaco, Argentina 

Ana María Giménez^*, Patricia Hernández\ 

Norfol Ríos^ y Federico Calatayu^ 

RESUMEN 

El objetivo del trabajo es contribuir al conocimiento del crecimiento de árboles individuales de 
Geoffreea decorticans (chañar), definir en función de la edad la evolución del incremento radial (IR), 
diámetro a la altura del pecho (DAR), sección normal (SN) y estimar el incremento medio (IMA) y anual 
(lA) con el fin proponer directrices para la correcta gestión de la especie. El estudio fue realizado en un 
bosque del Chaco Semiárido, Parque Los Quebrachos, Alberdi, Santiago del Estero, Argentina. Para el 
análisis epidométrico de fuste se derribaron 16 individuos que cumplían con la condición de: DAP 
superior a 15 cm; fuste sin defectos; de buena vitalidad. De los árboles derribados se extrajeron rodajas 
de 5 cm de espesor, a las alturas de 0,30 m; 1 ,30 m y final del fuste. La marcación y medición de anillos 
se efectuó con el Equipo Computarizado ANIOL y el programa CATRAS. Los datos fueron analizados 
con el software estadístico IFOSTAT. El chañar presenta anillos de crecimiento bien demarcados. El 
incremento radial promedio es 5,90 mm (1 ,79-16,55). La edad de culminación del crecimiento en DAP 
es a los 5 años y en SN se logra a partir de los 37 años. Este crecimiento permite definir el manejo en 
rotaciones cortas. 

PALABRAS CLAVE: 

Anillos de crecimiento, chañar, crecimiento radial, edad del árbol, manejo forestal. 


ABSTRACT 

The objective of this study was to improve the knowledge on the annual growth rates of Geoffreea 
decorticans (chañar), to define the evolution of annual growth increment (IR), diameter at breast height 
(DAP), normal section (DN) with age, to estímate mean annual (IMA) and annual increment (lA) in order 
to propose guidelines for proper management of the species. The study was carried out in a semiarid 
Chaco forest, at Quebrachos Park, Alberdi, Santiago del Estero Province, Argentina. For epidometric 
stem analysis, 16 vigorous trees were selected based on DAP greater than 15 cm and olean bole. From 
each tree cross-sections of 5 cm thickness were removed at 0,30 m; 1 ,30 m above ground, and at the 
end of the stem. Definition of ring boundaries and ring with measuring were made with computerized 
equipment ANIQL and the Computer program CATRAS. The data were analyzed with the statistical 
software package INFQSTAT. Chañar presentad well-defined annual growth rings. The average thick- 
ness of rings is 5,90 mm (1 ,79 - 16,55). The age of completion for DAP was determinad to be the five 
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years, and NS is achieved after 37 years. 
This growth rate can define the manage- 
ment on short rotations. 

KEY WORDS: 

Growth rings, chañar, radial growth, age 
of tree, forest management. 


INTRODUCCIÓN 

Las especies arbóreas del género Geo- 
ffroea que vegetan en Argentina son; 
Geoffroea striata Morong (manduvi-rá) 
y Geoffroea decorticans Burkart (chañar). 
Ambas tienen importancia florística y son 
usadas localmente, ya que sus maderas 
tienen escaso valor comercial. La madera 
y la corteza de esta especie han sido 
descritos en forma comparada y se con- 
sideran de alta afinidad (Giménez, 2004). 
El chañar es una especie arbórea secun- 
daria de amplia distribución geográfica, 
se extiende desde el sur de Perú, norte 
de Chile, Bolivia, Chaco Paraguayo y 
oeste de Uruguay, hasta la Patagonia 
Argentina. En este último país ocupa 
principalmente las provincias fitogeográ- 
ficas del Chaco, Monte y Espinal, desde 
Jujuy a Río Negro y Neuquén (Coscaron 
Arias y Ganmdullo, 2004). Puede encon- 
trarse como un integrante secundario del 
bosque chaqueño, o formando rodales 
puros en zonas abiertas, inundables, con 
niveles freáticos poco profundos y sali- 
nos. Es una especie colonizadora de 
suelos húmedos, ya consolidados y ricos 
en limos, muy agresiva en el Chaco 
Boreal (Mereles y Degen, 1997). En el 
Chaco semiárido es una de las especies 
de mayor frecuencia, presente en 80% 
de las áreas de estudio junto a Cercidium 
praecox, Prosopis nigra y Prosopis alba 
(Giménez y Hernández, 2008). 

El chañar se reproduce por semillas 
o por vía vegetativa mediante raíces 
gemíferas y cepa. Las semillas poseen 
diseminación endozoica. La activación de 
yemas en las ramas, corona del tallo o las 


raíces pueden ser provocadas, entre otros 
factores, por el fuego. Es alta la capacidad 
de rebrote después de sufrir incendios y 
cortes severos. 

Esta especie tiene dos aspectos a 
considerar, en algunas zonas de Argen- 
tina es reconocida como una leñosa 
arbustiva invasora; mientras que en bos- 
ques heterogéneos del Chaco semiárido 
es un componente secundario. Es apre- 
ciada localmente por el uso de su madera 
y otros valores (Giménez y Moglia, 2003). 
En las zonas semiáridas y subhúmedas 
forma densas colonias, debido a que su 
sistema radical da origen a un conjunto de 
plantas dispuesta en forma circular u oval, 
llamadas isletas, con plantas de mayor 
tamaño y edad, próximas al centro. En 
algunas regiones, se relaciona a la espe- 
cie con perjuicios que ocasiona al compor- 
tarse como invasora en pastizales 
degradados, sin considerar los importan- 
tes productos que ofrece (Anderson, 
1976). La concepción que la biomasa de 
chañar puede tener uso industrial, no ha 
sido debidamente desarrollada, con 
aspectos positivos complementarios, 
obteniendo ventajas productivas (gana- 
dera y forestal) económicas y ecológicas. 
En la implementación de sistemas silvo- 
pastoriles debe tomarse en cuenta la 
cobertura y la densidad de individuos, 
regulando mediante podas y raleos la 
entrada de luz y la circulación del ganado 
(Echeverría y Giulietti, 2001). 

Es de interés regional, investigar 
sobre el crecimiento y ecología del chañar 
con el fin de darle un adecuado manejo. 
Una herramienta válida para tal análisis es 
el estudio de crecimiento a partir del con- 
teo y medición de anillos. El ancho de los 
anillos de crecimiento es el parámetro que 
tradicionalmente se emplea para inferir 
cambios en la historia de vida de los árbo- 
les y la relación con su entorno. El empleo 
de los anillos de crecimiento de las leño- 
sas permite cuantificar en forma precisa 
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las velocidades de crecimiento radial de 
los árboles en distintos ambientes, así 
como determinar diferencias en las estruc- 
turas poblacionales del bosque en rela- 
ción con factores tales como el clima, el 
suelo, los disturbios naturales y la historia 
de uso (Villalba, 2000). Estudios que 
requieran información sobre la estructura 
de edad y la dinámica de los bosques, la 
determinación de sus turnos de aprove- 
chamiento, el establecimiento de las rela- 
ciones entre crecimiento leñoso, o 
simplemente la reacción del bosque a un 
tratamiento silvícola determinado, necesi- 
tan del fechado y la evaluación de los rit- 
mos pasados y presentes del crecimiento 
(Schweingruber, 1996). 

El conocimiento de la estructura y 
dinámica de estos bosques permitirá plan- 
tear estrategias de manejo y recuperación 
de los mismos (Villagra et al., 2002). La 
aplicación del análisis epidométrico en el 
manejo es directa y representa una herra- 
mienta indispensable en la planificación 
de aprovechamientos sustentables de 
bosques. Las tasas de crecimiento obser- 
vadas en las diferentes especies nativas 
estudiadas, indican que los turnos biológi- 
cos de corta son más largos de lo que se 
ha pensado y se viene aplicando hasta la 
actualidad (Giménez y Ríos, 1999; Gimé- 
nez et al., 2000; Giménez et al., 2003). La 
modelación del crecimiento del leño de 
especies nativas en su ambiente, permite 
la estimación de los momentos en que los 
árboles alcanzan el máximo crecimiento 
medio, información muy útil para una 
apropiada corta, asegurando la renova- 
ción del monte (Juárez de Galíndez y Bal- 
zarini, 2003). 


OBJETIVOS 

El objetivo del presente trabajo es contri- 
buir al conocimiento del crecimiento de 
árboles individuales de Geoffreea decortl- 
cans] definir la evolución del incremento 


radial (ir), diámetro a la altura del pecho 
(DAP), sección normal (sn) con la edad, esti- 
mar el incremento medio (ima) y el incre- 
mento anual (ia), evaluando las 
perspectivas de la especie para un futuro 
manejo sustentable. 

METODOLOGÍA 

El estudio fue realizado en el Parque Los 
Quebrachos, Departamento Alberdi, Pro- 
vincia de Santiago del Estero, Argentina. 
El área corresponde al Distrito Chaqueño 
Occidental y su posición geográfica es 
26°04’1 9,62” Lat. S.; 63°35’41 ,40” Long. O 
(Fig. 1). El clima es subtropical, continen- 
tal, de carácter semiárido, con precipita- 
ciones escasas y estivales alternando con 
sequías en períodos bien delimitados. La 
precipitación media anual es de 600 mm 
(para un intervalo de 500 mm a 700 mm). 

El área de estudio se localiza a esca- 
sos metros al este de la llanura de escurri- 
miento impedido del río Salado (bañados 
activos). Su localización geomorfológica 
coincide con la antigua llanura de inunda- 
ción, la cual fue cubierta por un manto 
superficial de sedimentos eólicos. Como 
consecuencia el material originario de los 
suelos se halla íntimamente relacionado 
al comportamiento del río Salado y a la 
posterior cobertura superficial. Los suelos 
son planos, suavemente elevados con 
vegetación boscosa (Piarfon, 2005). 

El predio en estudio es un bosque 
secundario de dos quebrachos, clausu- 
rado hace 20 años, lo que permitió a la 
vegetación expresarse con gran exube- 
rancia (Giménez y Hernández, 2008). La 
estructura vertical diferencia 4 estratos: 

1) herbáceo y gramíneo, de cobertura 
completa con predominio de Trichioris crl- 
nita; Solanum argentinum, Heliotropium 
sp., Setaria fiebrigl; Bromelia hieren y; 

2) arbustivo, de 3 m de altura compuesto 
por Celtis pallida, Acacia preccz; 3) arbó- 
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reo inferior, formado por especies secun- 
darias: Ziziphus mistol, Prosopis ruscifolia, 
Cercidium praecox, Prosopis nigra, Geof- 
froea decorticans (con alturas de 6 m a 8 
m); 4) arbóreo superior formado por Schi- 
nopsis lorenzii (quebracho colorado) y 
Aspidosperma quebracho-blanco (que- 
bracho blanco) con alturas de 12 m a 20 
m. Las leñosas más abundantes son: 
Schinopsis lorenzii (44%), Zizipphus mis- 
tol, Celtis tala, Prosopis nigra, Geoffreea 
decorticans (1 0%). Según datos del Inven- 
tario forestal, el chañar está representado 
en el rodal con 28 árboles/ha, volumen de 
fuste de 3,335 m%a y volumen total de 
6,11 m^/ha (Piarfon, 2005). 

El chañar está distribuido en peque- 
ños grupos. Es una especie de cuarta 


magnitud, corteza verdosa y flores amari- 
llas muy llamativas (Fig. 2). (Giménez y 
Mogia, 2003). La distribución por clases 
diamétricas indica un rodal joven, con 
escasos pies de diámetros superiores a 20 
cm, en clases de edad entre 5 y 20 años, 
con predominio de la clase diamétrica de 5 
cm a 10 cm. Para el análisis epidométrico 
de fuste, se derribaron 16 individuos de la 
especie Geoffreea decorticans (Gilí., ex 
Hook.& Arn.) Burk., Familia: Fabaceae 
(Leguminosae: Papilionoideae). 

De una muestra de árboles de diferen- 
tes edades tomada al azar, se escogieron 
los que cumplían la condición de: DAP 
superior a 15 cm; fuste sin defectos; de 
buena vitalidad. A cada árbol seleccionado 
se midió el dap a 1 ,30 m; altura de fuste (m) 


Madera y Bosques 19(1), 2013:37-51 


41 




Figura 2. Detalles dendrológicos del chañar. 


y altura total (m), marcándose la dirección 
norte en el fuste. De los árboles derribados 
se extrajeron secciones trasversales de 5 
cm de espesor a una altura sobre el suelo 
de 0,30 m; 1 ,30 m y al final del fuste, acorde 
con metodología de Helinska (1991). Las 
muestras se prepararon con cepilladora 
(cepillo), lijadora de banda y orbital con 
juego de lijas de granulometría de 600 gra- 
nos por pulgada cuadrada a 1000 granos 
por pulgada cuadrada, hasta tener una 
superficie lisa, que demarcara con alta defi- 
nición los anillos de crecimiento. 

La marcación y medición de anillos 
se efectuó con el Equipo Computerizado 
ANiOL y el programa cairas (Aniol, 1991). 
Con este procedimiento se consignó el 
número de años y el espesor del anillo con 
una resolución de 0,01 mm. Las medicio- 
nes se efectuaron de corteza a médula. 
Se sincronizaron los anillos manualmente 
y por cuadrantes, a fin de tener diferencia 
de edad (± 2 años). 


Se analizaron las siguientes varia- 
bles en función de la edad: 

ir: Incremento radial o espesor de anillo 

de crecimiento 

dar: Diámetro a la altura de 1 ,30 m 
sn: Sección normal del árbol (sección a 
la altura de 1,30 m). 


Se calculó el crecimiento medio anual 
(ima) y crecimiento anual (ia) para las 
variables dap y sn. El crecimiento medio 
anual (ima) resulta del cociente entre el 
valor total acumulado de la variable consi- 
derada, en un determinado tiempo de su 
evolución, desde su origen y su edad en 
dicho momento. Incremento anual (ia) es 
la expresión del crecimiento anual que 
resulta del crecimiento en un año. El creci- 
miento anual y el crecimiento medio anual 
coinciden en un determinado momento de 
la vida del árbol o de la masa, en el cual el 
crecimiento medio anual es máximo. La 
edad de culminación biológica se obtuvo a 
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partir del punto de encuentro de la curva 
del iMA y del ia. 

IMA= Variable(n) / n 
IA= Vn-V(n-1) 

donde: 

ima: incremento medio anual 

(n): año 

ia: Incremento anual 

Vn: variable en el año n 

V(n-1): variable en el año n-1 

Análisis estadístico 

Los datos fueron analizados con el soft- 
ware estadístico programa ifostat 2010. 
Para el estudio de las tendencias en creci- 
miento se usó el método basado en la edad 
biológica de los árboles (Komin, 1967), tra- 
bajándose con individuos de edades dife- 
rentes. Se empleó la técnica de suavizado 
de las series individuales de ancho de ani- 
llos para eliminar la variación climática de 
las series de IR, previa modelación del cre- 
cimiento. Se usó un filtro basado en medias 
móviles de ventana igual a 5 años (Juárez 
de Galíndez et al., 2007). 

El espesor de los anillos de cada uno 
de ellos, se alineó con los de otros indivi- 
duos según la edad biológica. Así se 


expresaron biológicamente las tendencias 
de crecimiento de la especie en determi- 
nado sitio. La estandarización de la curva 
de crecimiento medio presume que la 
forma de la estructura a cualquier edad 
biológica es independiente del período de 
tiempo durante el cual se produce. 

Fueron analizados los modelos de 
regresión para las propiedades estudiadas 
y la edad. Se trabajó con un modelo polino- 
mial de segundo orden para modelar la 
estructura de la media (promedio) de la 
variable ir en función de la edad del árbol. 

Modelo de crecimiento: 

y= Bo+B.,.Edad+B 2 *Edad 2 + E^ 

RESULTADOS 

El chañar es una de las pocas especies 
arbóreas del Chaco argentino que pre- 
senta anillos de crecimiento bien defini- 
dos, lo que permite realizar el análisis 
epidométrico del fuste por conteo y medi- 
ción de anillos. A nivel macroscópico se 
demarca una banda oscura correspon- 
diente a la sección de madera tardía, 
donde el tejido es más denso, seguida de 
una banda ancha y clara de madera tem- 
prana al inicio del anillo (Fig. 3). 



Figura 3. Vista macroscópica de los de anillos de chañar 
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Figura 4. Vista microscópica del inicio de un anillo de chañar. 


La definición del anillo está dada por 
la conjunción de dos aspectos anatómi- 
cos: la presencia de una banda de parén- 
quima terminal, más la distribución de los 
vasos a lo largo del anillo en porosidad 
semicircular (Giménez, 2004) (Fig. 4). 

En la tabla 1 se indican los estadísti- 
cos descriptivos de las series promedio de 
IR de 16 individuos. El espesor medio de 


anillos (IR) fue de 5,90 mm, con una inter- 
valo de 1,79 mm a 16,55 mm, desviación 
estándar de 2,54 y el coeficiente de varia- 
bilidad de 43%. Los árboles estudiados 
mostraron edades que oscilan entre 9 a 
33 años. 

Los valores ir fueron variables entre 
los individuos de un mismo rodal (Fig. 5) 
con un valor máximo a los 5 años de edad. 



Figura 5. Relación entre las series IR y edad para el chañar. 
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Tabla 1. Estadística Descriptiva de las series de IR de 16 árboles de chañar. 


IR 

Edad del árbol 
(años) 

Media 

(mm) 

Desviación 

estándar 

(mm) 

Valor mínimo 
(mm) 

Valor máximo 
(mm) 

Árbol 1 

23 

5,05 

2,31 

2,13 

10,93 

Árbol 2 

19 

5,23 

1,63 

2,84 

8,88 

Árbol 3 

20 

4,87 

1,27 

2,63 

7,30 

Árbol 4 

33 

4,87 

1,56 

2,49 

8,99 

Árbol 5 

17 

9,75 

3,35 

5,06 

16,52 

Árbol 6 

15 

7,45 

3,2 

2,97 

12,61 

Árbol 7 

21 

6,27 

1,55 

4,06 

10,47 

Árbol 8 

16 

6,13 

2,27 

1,79 

10,49 

Árbol 9 

16 

5,3 

2,05 

2,77 

9,50 

Árbol 10 

29 

5,14 

1,47 

2,52 

7,82 

Árbol 11 

21 

4,45 

2,05 

2,07 

10,19 

Árbol 12 

10 

11,09 

3,49 

6,89 

16,55 

Árbol 13 

10 

7,48 

1,95 

4,15 

10,63 

Árbol 14 

9 

8,57 

2,27 

4,48 

11,98 

Árbol 15 

26 

5,43 

1,84 

2,72 

9,77 

Árbol 16 

20 

5,73 

1,52 

2,81 

9,06 


En la figura 6 se indica el histograma 
de frecuencia relativa de ir. 

En la tabla 2 se indican las funciones 
de ajuste de las variables ir, dar y sn 

Dada la variabilidad del ir, las series 
se suavizaron con medias móviles en 
escala de 5 años. Se ajustaron a una 


curva explicada por una función polinó- 
mica de segundo grado con un de 0,51 
(Fig.7). 

La evolución DAP/edad (Fig. 8) se 
ajustó a una ecuación polinómica de 
segundo grado con un R^ de 0,84 (30 cm 
DAP a los 31 años). 


Madera y Bosques 19(1), 2013:37-51 


45 



IR mm 


Figura 6. Frecuencia relativa (%) del ir (mm). 


Tabla 

2. 


Estadística 

descriptiva y función de ajuste 

Variable 


N 


Función de ajuste 

IR 


329 

0,51 

IR=8,40E-04 edad^-0,16 edad +7,55 

DAP 


329 

0,84 

DAP= -0,013 edad^ + 13,57 edad + 4,87 

SN 


329 

0,79 

SN = 38,4 edad^ + 1222,07 edad -1732,33 



Edad (años) 

Figura 7. Ajuste del ir suavizado. 
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El incremento promedio en diámetro 
para un intervalo de 33 años fue de 0,879 
cm. Este dato es excepcional en las leño- 
sas arbóreas del Chaco Semiárido, ya que 
se alcanza el diámetro mínimo de corta a 
tan temprana edad. Se estimó el ima e ia 
en función del dap (Fig. 9), a los 33 años; 
el IMA es de 9,16 mm y lAde 5,82 mm. Se 


encontró que el ima e ia se interponen tem- 
pranamente a los 5 años a un valor de 
11,81 mm. 

La sección normal en relación con 
la edad se ajustó a una función polinó- 
mica de segundo grado con de 0,79 
(Fig. 10). 



Figura 8. Evolución del dap y la edad. 
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Figura 9. Culminación del incremento en dap según ima e ia. 
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Se obtuvieron valores calculados y 
ajustados de ima, 2738 mmm, y de ia, 
2176 mm^, a la edad de 33 años (Fig. 
11). Las curvas están prontas a intercep- 
tarse, lo que indica una culminación 
temprana del crecimiento expresado en 
sección, con máximos entre los 20 y 25 
años. 


Para definir el turno biológico de 
corta de una población, se determinó la 
edad de la planta donde el crecimiento 
anual (crecimiento promedio en una edad 
determinada) es menor al crecimiento 
anual promedio alcanzado hasta ese 
momento (crecimiento acumulado dividido 
por la edad de la planta). 



Fig. 10. Evolución de sn y edad. 
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Figura 11 . ima e ia en función de la sn. 
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Figura 12. Tiempo de tránsito por clases diamétricas. 


Como la variable que mejor se corre- 
laciona con el volumen económico es la 
sección normal (sn), se estimó el incre- 
mento anual en sn para cada edad a partir 
del incremento diamétrico medido. Con la 
ecuación de ajuste, se proyectó la evolu- 
ción de la SN a través de sus incrementos, 
lográndose interceptar las curvas de ima e 
lA a la edad de 37 años para un valor de 
2240 mm^.Esto permitió especular sobre 
los diámetros máximos citados por la lite- 
ratura y el punto de culminación del creci- 
miento en SN. La edad de culminación 
expresada en sn fue de 37 años en el 
caso de este estudio. Se alcanzó el diá- 
metro mínimo de corta a los 25 años. 

Con la curva de dap y edad se estimó 
el tiempo de tránsito entre cada clase dia- 
métrica de 5 cm, que varió entre 6 a 10 
años según las clases (Fig.12). 

DISCUSIÓN 

Existe una literatura rica en relación con 
estudios del espesor medio de los anillos 
de chañar. Roth y Giménez (2006) indi- 
can un valor de 3,1 mm (1 mm a 5 mm) 
para el Chaco Semiárido. Echeverría y 
Giulietti (2001) citan un espesor medio 


de 3,4 mm/ año para San Luis, con creci- 
miento en diámetro de 0,65 cm/año. En 
chañarales de Chile, Serra (1998) revela 
un crecimiento en diámetro de 0,7 cm. En 
este estudio el espesor medio de anillos 
para 33 años de intervalo fue de 5,90 
mm/año.Se notó una marcada diferencia 
en la velocidad de crecimiento con la 
edad. El crecimiento fue muy rápido 
durante los primeros 5 años de vida y 
alcanzó un máximo entre los 2 a 5 años 
(5,9 mm) para decrecer a los 30 años, lo 
que indicó la culminación temprana de 
crecimiento. 

Juárez de Galíndez y Balzarini 
(2003), trabajando con especies nativas, 
sugieren la técnica de suavizado aplicada 
a las series de anchos de anillos para eli- 
minar la variación climática de las series 
individuales, previa modelación del creci- 
miento. Los árboles objeto de estudio indi- 
caron un comportamiento de una especie 
de rápido crecimiento y longevidad baja 
para las condiciones ecológicas del 
Chaco. Este resultado implica la necesi- 
dad de replantear la mayoría de los 
emprendimientos forestales existentes, 
incluso la legislación forestal. Tradicional- 
mente se ha aconsejado cortar árboles 
mayores a 25 cm de diámetro para espe- 
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des secundarias, esto significa cortar 
árboles antes de que alcancen la mayor 
tasa de crecimiento forestal. El diámetro 
mínimo de corta se logra a los 25 años. 

En referencia a los anillos de creci- 
miento, las características de demarca- 
ción y el espesor varía según las especies 
de la región. En Prosopis alba el espesor 
medio de anillos es 4,05 mm (2,11 mm a 
5,18 mm) (Giménez et al., 1998); Prosopis 
kuntzel, 3,26 mm (1,9 mm a 4,5 mm) 
(Giménez et al., 1997); Prosopis nigra, 
3,29 mm (1,14 mm a 8,59 mm) (Giménez 
et al., 2003); Prosopis ruscifolla, 5,66 mm 
(3,3 a 6,8 mm) (Giménez et al., 2009); 
Schinopsis quebracho-colorado, 2,19 mm 
(1,6 mm a 4,5 mm) (Giménez y Ríos, 
1 999) y Aspidosperma quebracho-blanco 
con 2,20 mm (Moglia, 2000).Juárez et al. 
(2007) indican para vinal {Prosopis rusel- 
folla) una edad de culminación del incre- 
mento radial a los 11 años de edad, 
mientras que para la especie de este estu- 
dio fue de cinco años. 

Como se evidencia, Geoffroea decor- 
tlcans es la especie de mayor IR promedio 
de las estudiadas para el Chaco Semiá- 
rido, lo cual representa una aliciente muy 
importante. Hay que destacar que el diá- 
metro máximo de 30 cm, citado por Torto- 
relli (2009), es indicativo de la baja 
longevidad de la especie. 


CONCLUSIONES 

Geoffreea decortieans es una especie de 
crecimiento medio y de muy buenas pers- 
pectivas de aprovechamiento en relación 
con el tiempo de tránsito de las clases dia- 
métricas.EI Incremento radial promedio 
(ir) es de 5,90 mm (1,79 mm -16,55 mm), 
con un máximo a edades tempranas (5 
años), para un intervalo de 35 años. La 
edad de culminación expresada en DAR 
fue a los 5 años, en SN se logra a partir de 
los 37 años; el diámetro mínimo de corta 


se logra a los 25 años. Este crecimiento 
permite el manejo en rotaciones cortas. 
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Estimación de contenido de carbono 
en la cuenca del río Magdalena, México 

Estimation of carbón contení 
in the Magdalena River basin, México 

José Mauricio Galeana-Pizaña\ 

José Antonio Benjamín Ordóñez-Díaz^ y Nirani Corona-Romero^ 

RESUMEN 

En el presente estudio se estimó el contenido de carbono de la biomasa aérea de seis tipos de 
cobertura vegetal en la cuenca del río Magdalena, ubicada en el suelo de conservación del Distrito 
Federal. Para ello, se elaboró un mapa de las clases de cobertura vegetal y uso de suelo del área de 
estudio, mediante la fotointerpretación de fotografías aéreas. Se utilizó un diseño de muestreo jerár- 
quico estratificado con distribución sistemática, los sitios resultantes se desplegaron en un enmallado 
con una equidistancia de 550 metros. La colecta de los datos en campo se basó en un diseño anidado; 
con una intensidad de muestreo de 0,64%. Se establecieron 52 sitios de muestreo de 3000 m^ cada 
uno. La cobertura vegetal que presentó mayor cantidad de carbono almacenado corresponde a la clase 
de Abies religiosa cerrado, con 163 Mg C ha'^ promedio; le sigue la clase de Abies religiosa - Pinus 
hartwegii cerrado con 107 Mg C ha'^ y la clase de Abies religiosa abierto con 95 Mg C ha T Las clases 
de Pinus sp. y Pinus hartwegii cerrado tuvieron 47 Mg C ha'^; Pinus hartwegii - Abies religiosa cerrado, 

39 Mg C ha-1 y Pinus hartwegii abierto, 21 Mg C ha T 

PALABRAS CLAVE: 

Abies religiosa, biomasa aérea, bosques, estructura diamétrica, Pinus hartwegii. 

ABSTRACT 

Carbón content of aboveground biomass of six types of vegetation cover in the Magdalena River 
basin, located in the conservation land in the México City, was estimated. Henee, a map of cover vege- 
tation and land use to the study area was made; a hierarchical stratified sampling design with distribu- 
tion systematic with a mesh of sampling sites equidistant of 550 meters was used for fieid work. The 
method used for fieid data collection was based on a nested design, with an intensity of sample of 
0,64%. There were generated 52 sites of 3000 m^. The vegetation cover that showed the highest 
amount of carbón stored corresponde to the forest of the the Abies religiosa closed with 163 Mg C ha'^ 
in average; followed by the Abies religiosa - Pinus hartwegii closed with 107 Mg C ha \ and Abies reli- 
giosa opened with 95 Mg C ha'^; about the covers of Pinus sp., Pinus hartwegii closed with 47 Mg C ha'^; 
Pinus hartwegii - Abies religiosa closed with 39 Mg C ha'^ and Pinus hartwegii opened with 
21 Mg C ha T 

KEYWORDS: 

Abies religiosa, above ground biomass, forest, diametrical structure, Pinus hartwegii. 
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INTRODUCCIÓN 

El cambio climático se define como las 
modificaciones termodinámicas que ocu- 
rren en la atmósfera a escala global, como 
consecuencia del incremento de las con- 
centraciones de gases de efecto inverna- 
dero (GEi), provocando el calentamiento 
del planeta (ipcc, 2007). Esto afecta direc- 
tamente al hombre y a las actividades 
económicas que tienen una correlación 
con el estado del tiempo (e.g.: tempera- 
tura, presión atmosférica, precipitación, 
etc.) produciendo modificaciones en su 
medio (e.g. aumento del nivel de mar), lo 
que a largo plazo provocará una crisis 
social (Galeana, 2008). 

En el año 1976, se reconoció a los 
bosques como uno de los principales 
sumideros de carbono, junto con los océa- 
nos, el suelo y los pastizales (Mena, 
2004). Actualmente, los ecosistemas bos- 
cosos están siendo amenazados por el 
cambio de uso de suelo, la deforestación y 
la quema de combustibles fósiles, siendo 
éstas las principales causas a nivel mun- 
dial de la liberación de dióxido de carbono 
(COJ a la atmósfera (Schneider, 1989; 
Goudie, 1990; Dixon et al., 1994). Los 
cambios en el uso del suelo y cobertura 
vegetal asociados a la fragmentación de 
hábitats son los disturbios más severos de 
las actividades antrópicas, causando un 
aumento de las tasas de extinción de 
especies, sobre todo en las décadas 
recientes (Soulé y Orians, 2001). Lo ante- 
rior indica un severo problema que se está 
acentuando y que se relaciona directa- 
mente con el cambio en la cobertura vege- 
tal, lo cual repercute en el deterioro del 
ambiente (Lambin, 1994; Ordóñez, 2008). 

La deforestación genera, a nivel 
regional, pérdida de la cubierta vegetal, 
detrimento de los recursos forestales, 
deterioro físico y químico del suelo, altera- 
ción del balance hídrico y desestabiliza- 
ción de cuencas; a nivel global, afecta el 


albedo y el balance de agua atmosférica, 
perturbando los patrones climáticos, con- 
tribuyendo así al calentamiento global 
(Houghton, 1991; Johns et al., 1 992; Lund, 
2006). Otra consecuencia es la reducción 
y fragmentación del hábitat (Murcia, 1995) 
lo que se traduce en pérdida de la biodi- 
versidad y eliminación de la variabilidad 
genética de poblaciones y especies 
(Johns et al., 1992), sumado a una pér- 
dida de los servicios ambientales y biocul- 
turales (Otegui, 2007). Es por ello que 
estudios relacionados a los cambios en la 
cobertura vegetal están teniendo mayor 
atención (Brown y Lugo, 1994). 

Históricamente, en una escala global, 
los cambios en la cobertura forestal han 
sido portadores de CO 2 a la atmósfera 
(Houghton, 1991; Brown, 2001), por esta 
razón en la actualidad se considera a los 
bosques como una opción para mitigar las 
emisiones de gei, ya que secuestran y 
mantienen más CO 2 que otros ecosistemas 
terrestres, aparte de participar con 90% del 
flujo anual de carbono en el sistema atmós- 
fera-tierra (Masera, 1995). El ipcc (2007) 
estima que combinando estrategias de 
conservación forestal con proyectos de 
reforestación en el mundo, los bosques 
podrían resultar en un sumidero neto de 
carbono durante los próximos 100 años, 
permitiendo reducir de 20% a 50% las emi- 
siones de CO 2 a la atmósfera. 

Las opciones de mitigación de los gei 
pueden ser agrupadas en tres: La primera 
categoría incluye a las actividades que evi- 
tan la liberación del carbono secuestrado 
(e.g.: conservación y protección a los eco- 
sistemas forestales); la segunda contiene 
actividades que capturan CO 2 (e.g.: afores- 
tación, reforestación y agroforestería); y la 
tercera opción encierra actividades que 
sustituyan el uso de combustibles fósiles 
(Sathaye et al., 2001). Los procesos de 
captura-emisión son parte de un sistema 
con cuatro tipos generales de reservones 
de carbono (vegetación aérea y subterrá- 
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nea, materia en descomposición, suelos y 
productos forestales), con tiempos de resi- 
dencia y flujos asociados muy diferentes 
(Ordóñez, 1999). 

Actualmente es creciente la atención 
global sobre el estimar con mayor certeza 
y precisión la dinámica del carbono 
secuestrado y emitido (balance emisión- 
captura) entre el bosque y la atmósfera. 
Reconocer el papel que juegan los bos- 
ques en el ciclo global del carbono, parti- 
cularmente en la mitigación de los gei, es 
una prioridad (Brown, 2001). La estima- 
ción de la captura de carbono es el punto 
de partida para hacer efectivas las opcio- 
nes de mitigación de los gei en el marco 
del Protocolo de Kyoto, además de permi- 
tir establecer un criterio de base de 
cálculo de la dinámica de este gas cuando 
no se ejecuta una medida de mitigación 
(Masera et al., 2000); es decir, generar un 
caso de referencia con el cual se pueda 
comparar el cambio en las emisiones de 
los GEI y verificar los beneficios adiciona- 
les posteriores (ipcc, 2007). 

La cuenca del río Magdalena brinda 
servicios ambientales a la ciudad de 
México, como son captura de carbono, 
recarga de mantos acuíferos, producción 
de oxígeno, fijación del sustrato y belleza 
escénica, principalmente; lo anterior 
sugiere acciones inmediatas dirigidas a la 
conservación y adecuado manejo de sus 
recursos. Sin embargo, la zona ha recibido 
fuertes presiones antropogénicas como 
asentamientos humanos, procesos de 
cambio de uso de suelo, instalaciones de 
infraestructura para diversos fines, extrac- 
ción de recursos naturales (madera y 
fauna nativa, principalmente) y turismo no 
regulado, lo que ha ocasionado, entre 
otros problemas: deforestación y degrada- 
ción forestal, establecimiento y auge de 
plagas forestales, incendios, contamina- 
ción del acuífero, inseguridad y erosión. 
Por ello, es primordial el desarrollo de una 
línea base confiable para iniciar proyectos 


de captura de carbono en los bosques del 
río Magdalena, lo que permita entender el 
proceso del cambio en la cobertura vege- 
tal y uso del suelo y así poder ofrecer alter- 
nativas de mitigación del cambio climático 
(Fernández et al., 2002; Galeana, 2008). 


OBJETIVOS 

Los objetivos de la presente investigación 
son , estimar el contenido de carbono en la 
biomasa aérea de las coniferas de la 
cuenca del río Magdalena, Distrito Fede- 
ral, México; elaborar un mapa de cober- 
tura vegetal y uso del suelo del área de 
estudio y, cuantificar el contenido de car- 
bono por tipo de cobertura vegetal. 

METODOLOGÍA 
Área de estudio 

La cuenca del río Magdalena se localiza 
dentro del suelo de conservación de las 
delegaciones Magdalena Contreras, Alvaro 
Obregón y Cuajimalpa, en el Distrito Fede- 
ral (zona centro de México) (Fig. 1). 

El área de estudio se encuentra ubi- 
cada en el suroeste de la cuenca de México, 
sobre la Sierra de las Cruces, misma que 
pertenece a la región fisiográfica del Sis- 
tema Volcánico Transversal, dentro de las 
subprovincia de Lagos y Volcanes de Aná- 
huac (iNEGi, 2005a, b,c,d). Sus coordenadas 
extremas son: al Norte 19°17'52”, al Sur 
19°14'35”, al Oeste 99°20'30” y al Este 
99°15'06”(gdf, 2000). 

La cuenca tiene una superficie de 
2981,3 ha de relieve accidentado (2480 
msnm-3860 msnm) con bosques de Ables 
religiosa, PInus hartwegll, Quercus spp. y 
mixto; sus climas son templado subhú- 
medo y semifrío subhúmedo, ambos con 
lluvias en verano y suelos de tipo andosol 
húmico (GDF, 2000). 
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Figura 1. Localización de la cuenca del río Magdalena. 


En la cuenca se forma el río Magdalena, 
que es uno de los cuerpos de agua más 
importantes de la Ciudad de México, con 
un volumen aproximado de agua perma- 
nente de 1 m^s'^ y que es empleado como 
fuente de abastecimiento de agua potable 
(CNA, 1997). 


Elaboración del mapa de vegetación 
y uso de suelo 

Para realizar la estimación de captura 
de carbono se determinaron la clases de 
cobertura vegetal y uso de suelo para la 
elaboración del mapa, las cuales se defi- 
nieron con base en los criterios de 
Reporte Técnico del Inventario Nacional 
Forestal (inf) de la Subsecretaría Fores- 
tal y Fauna Silvestre (sffs, 1994), iden- 
tificando las mismas categorías de 
cobertura vegetal y uso de suelo. Para 
la elaboración del mapa se reagruparon 
las clases de cobertura vegetal en 26 
(Tabla 1). 

Una vez determinadas las clases de 
cobertura vegetal, se realizó la interpreta- 


ción de fotografías aéreas de escala 
1:10 000, tomadas en enero del año 2005 
con dirección norte-sur, y proporcionadas 
por el Instituto de Geografía de la Univer- 
sidad Nacional Autónoma de México. La 
digitación del mapa se generó en pantalla 
usando el software Autocad Map 2002 y la 
base de datos en ArcView 3.2a. 


Diseño de muestreo 

Se establecieron ocho clases de vegeta- 
ción a muestrear para los fines de este 
estudio (Tabla 2). 

Para estimar el contenido de carbono 
se utilizó un diseño de muestreo sistemá- 
tico simple. Esto se hizo sobreponiendo 
una red cuadrangular de puntos equidis- 
tantes a 550 m (Fig. 2) a la capa de vege- 
tación y uso de suelo, con el programa Are 
GIS 9.2. El tamaño de muestra se deter- 
minó con 1% de intensidad de muestreo. 
Los puntos estuvieron definidos por coor- 
denadas Universal Transversal de Merca- 
tor (UTM), representando el centro de las 
unidades a muestrear. 
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Tabla 1 . Relación de las clases de cobertura vegetal con las establecidas por el inf, 1 994. 


Clases en este estudio Claves Clases del INF 


Agricultura 

Asentamientos humanos 

Bosque de Abies religiosa abierto 

Bosque de Abies religiosa cerrado 

Bosque de Abies reiigiosa-Pinus hartwegii abierto 

Bosque de Abies reiigiosa-Pinus hartwegii cerrado 

Bosque de Abies reiigiosa-Pinus sp. abierto 

Bosque de Abies reiigiosa-Pinus sp. cerrado 

Bosque de Abies religiosa-Quercus sp. abierto 

Bosque de Abies religiosa-Quercus sp. cerrado 

Bosque de Quercus sp. abierto 

Bosque de Quercus sp. cerrado 

Bosque de Quercus sp.-Abies religiosa abierto 

Bosque de Quercus sp.-Abies religiosa cerrado 

Bosque de Pinus hartwegii abierto 

Bosque de Pinus hartwegii cerrado 

Bosque de Pinus sp. abierto 

Bosque de Pinus sp. cerrado 

Bosque de Pinus hartwegii -Abies religiosa abierto 

Bosque de Pinus hartwegii -Abies religiosa cerrado 

Bosque de Pinus sp.- Abies religiosa abierto 

Bosques de otras coniferas abierto 

Bosque de Pinus sp.- Abies religiosa cerrado 

Bosque mixto 

Cuerpos de agua 

Pastizal 

Plantación forestal 


AGR 

Agricultura de riego 
Agricultura de temporal 

AS 

Zona urbana 

BA^ 

Bosque de oyamel abierto 

BA 2 

Bosque de oyamel cerrado 

BA-Ph^ 

Bosques de otras coniferas abierto 

BA-Ph2 

Bosques de otras coniferas cerrado 

BA-Psp.i 

Bosques de otras coniferas abierto 

BA-Psp.i 

Bosques de otras coniferas cerrado 

BA-Q^ 

Bosques de otras coniferas abierto 

BA-Q 2 

Bosques de otras coniferas cerrado 

BQ^ 

Bosque de encino abierto 

BQ 2 

Bosque de encino cerrado 

BQ-A^ 

Bosques de otras coniferas abierto 

BQ-A 2 

Bosques de otras coniferas cerrado 

BPh^ 

Bosque de pino abierto 

BPh2 

Bosque de pino cerrado 

BPsp., 

Bosque de pino abierto 

BPSP 2 

Bosque de pino cerrado 

BPh-A^ 

Bosques de otras coniferas abierto 

BPh-A2 

Bosques de otras coniferas cerrado 

BPsp-A.i 

Bosques de otras coniferas abierto 

BPSP-A 2 

Bosques de otras coniferas cerrado 

BM 

No considerado 

H 2 O 

Cuerpos de agua 

PAS 

Pastizal inducido 

REF 

Plantación forestal 


Tabla 2. Clases de cobertura vegetal utilizadas para la estimación de carbono. 


Clases a muestrear 

Claves 

Clases del INF 

Bosque de Abies religiosa abierto 

BA, 

Bosque de oyamel abierto 

Bosque de Abies religiosa cerrado 

BA 2 

Bosque de oyamel cerrado 

Bosque de Abies reiigiosa-Pinus hartwegii cerrado 

BA-PÓ2 

Bosques de otras coniferas 
cerrado 

Bosque de Pinus hartwegii abierto 

BPh^ 

Bosque de pino abierto 

Bosque de Pinus hartwegii cerrado 

BPÓ2 

Bosque de pino cerrado 

Bosque de Pinus hartwegii-Abies religiosa cerrado 

BPh-A2 

No considerado 
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Figura 2. Red cuadrangular de puntos, obtenida del muestreo sistemático simple 

para la zona. 


Unidad de muestreo 

La metodología empleada en campo para la 
colecta de los datos de biomasa aérea se 
basó en un diseño anidado (Ordóñez, 2008; 
Fig. 3), que permitió la colecta de muestras 
y sus réplicas. De acuerdo con una intensi- 
dad de muestreo de 0,64%, se generaron 
52 sitios de 3000 m^, los cuales estuvieron 
conformados por tres parcelas circulares de 
1000 m^ con un radio de 17,84 m. 


Estimación del contenido de carbono 
en biomasa aérea 

Para estimar la biomasa aérea (ba) de los 
árboles se establecieron parcelas circula- 
res de 1 000 m^. En cada parcela se eva- 
lúo la vegetación arbórea. A cada individuo 
(h>1 ,30m) se le midió el diámetro normali- 
zado (dn) y la altura total. Los individuos 
con una altura total inferior a los 1 ,30 m se 


contabilizaron e identificaron por género. 
Todos los datos fueron registrados por 
género y por clase diamétrica, con interva- 
los de 5 centímetros. 

Con los datos recolectados sobre la 
especie Abies religiosa se elaboró la 
siguiente ecuación alométrica para las 
alturas: 

Abies religiosa 

y = -0.004(DN^) + 0,8427(DN) 
r^: 0,8427 

Para las alturas de la especie Pmus hart- 
wegii se utilizó la ecuación alométrica pro- 
puesta por Espinoza (2005): 

Pin US hartwegii 
y = 0,001 5(DN^) + 0.3956(DN) 

P: 0,8883 
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Figura 3. Diseño anidado. 


En la estimación de la biomasa aérea para 
ambas especies se utilizó la siguiente 
ecuación (Brown y Lugo, 1994): 

BA: volumen x densidad de la madera 

donde: 

Volumen: Área basal x altura 

(Grijpma, 2001) 

Área basal: n/4 (DN^) 

(Mostacedo y Fredericksen, 2000) 


Para fines de este estudio se tomó el dato 
de densidad de la madera (t/m^) promedio de 
las monografías expedidas por el Instituto 
Nacional de Investigaciones Forestales, 
Agrícolas y Pecuarias (Musáiem y Ramírez, 
2003; Musáiem y Solís, 2000). 


Todos los valores de biomasa se 
pasaron a carbono mediante un factor de 
0,5 MgDM / MGC (IPCC, 1995). 


RESULTADOS 

Mapa de vegetación 

Las superficies cubiertas por cada una de 
las clases de cobertura vegetal y uso 
de suelo se muestran en el tabla 3. 

El mapa de las diferentes clases de 
cobertura vegetal y uso de suelo (cvus) 
se muestra en la figura 4. Para validar 
cada una de las diferentes cvus se toma- 
ron 266 puntos de verificación bis a bis 
(usando tres sistemas de posiciona- 
miento global). 
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Tabla 3. Superficie que ocupan las diferentes clases de cobertura vegetal y uso de suelo 

en la crm. 


Clases 

Área (ha) 

% 

Agricultura 

31,98 

1,07 

Asentamientos humanos 

99,60 

3,34 

Bosque de Abies religiosa abierto 

183,08 

6,14 

Bosque de Abies religiosa cerrado 

1012,09 

33,95 

Bosque de Abies religiosa-Pinus hartwegii abierto 

19,69 

0,66 

Bosque de Abies religiosa -Pinus hartwegii cerrado 

386,21 

12,95 

Bosque de Abies religiosa-Pinus sp. abierto 

16,49 

0,55 

Bosque de Abies religiosa-Quercus sp. cerrado 

11,38 

0,38 

Bosque de Quercus sp. cerrado 

59,73 

2,00 

Bosque de Quercus sp.-Pinus sp. cerrado 

17,27 

0,58 

Bosque de Pinus hartwegii abierto 

322,14 

10,80 

Bosque de Pinus hartwegii cerrado 

288,56 

9,68 

Bosque de Pinus hartwegii-Abies religiosa cerrado 

192,19 

6,45 

Bosque de Pinus sp. abierto 

20,57 

0,69 

Bosque de Pinus sp. cerrado 

6,47 

0,22 

Bosque de Pinus sp.-Abies religiosa cerrado 

12,21 

0,41 

Bosque mixto 

143,24 

4,80 

Cuerpos de agua 

0,17 

0,01 

Pastizal 

164,15 

5,51 

Plantación forestal 

6,51 

0,22 

Total 

2981,53 

100 


Estructura diamétrica del bosque 

Se definieron 11 clases diamétricas para 
cada uno de los seis tipos de coberturas 
vegetales muestreadas, a fin de describir 
la estructura horizontal de los bosques 
de la cuenca del río Magdalena 

Para describir la estructura del bos- 
que, fueron medidos los diámetros de 
13,863 árboles y las alturas de 300 de 
ellos, de los cuales, 61 % fueron Pinus sp., 
28% Abies religiosa, 7% hojosas, 3% 
Cupresus sp. y 1% Quercus sp. (Fig. 5). 

En la figura 5 se observa que la clase 
diamétrica con más individuos es < 5 cm 
(50%), teniendo una mayor abundancia 
de renuevos de Pinus sp. (72%), seguido 


de Abies religiosa (22%) y en menor 
medida las hojosas (3%), Cupressus sp. 
(2%) y Quercus sp. (1%). La clase diamé- 
trica con menor número de individuos es 
42,6 cm -47,5 cm (1 ,8%) y está compuesta 
principalmente por Abies religiosa (71 %) y 
Pinus sp. (23%) y en menor medida por 
Cupressus sp. (4,5%), hojosas (0,5%) y 
Quercus sp. (0,5%). 

Treinta y cinco por ciento de los indi- 
viduos se encuentra en las clases diamé- 
tricas de 5 cm a 22,5 cm (Fig. 5), donde 
se presenta una mayor abundancia de 
Pinus sp. (60%), seguida de la especie 
Abies religiosa (21%), hojosas (14%) y 
en menor medida Cupressus sp. (4%) 
y Quercus sp. (1%). 
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Figura 5. Número de individuos totales, por clase diamétrica, de los bosques de la 
cuenca del río Magdalena, derivado de 52 sitios de muestreo. 


Las clases diamétricas restantes agrupan 
a 1 5% de los individuos. Abies religiosa es 
la de mayor abundancia con 63%, seguida 
de Pinus sp. (30%), Cupressus sp. (4%), 
hojosas (2,5%) y Quercus sp. (0,5%). 


Contenido de carbono 

Se determinó el carbono almacenado en 
cada una de las coberturas vegetales 
muestreadas (Tabla 4). 


Tabla 4. Contenido de carbono por cobertura vegetal. 


Cobertura vegetal Superficie (ha) n Media y EE Intervalo 

(tC/ha) (tC/ha) 


Bosque de Abies religiosa abierto 
Bosque de Abies religiosa cerrado 
Bosque de Abies - Pinus sp. cerrado 
Bosque de Pinus hartwegii abierto 
Bosque de Pinus hartwegii cerrado 
Bosque de Pinus sp. - Abies religiosa cerrado 


183,08 

4 

95,49 ± 24,51 

40,71 - 

155,82 

1012,09 

22 

162,91 ± 16,69 

72,37 - 

362,20 

341,7 

7 

106,92 ±25 

40,30 - 

205,39 

322,14 

5 

21,21 ± 5,62 

6,25- 

- 32,29 

288,56 

10 

47,37 ± 5,64 

18,37- 

- 67,84 

179,98 

4 

39,39 ± 14,88 

27,27 - 

108,73 


Superficie = área en hectáreas que ocupa la clase vegetal, 
n = número de sitios (cada sitios consta de 3 levantamientos de 0,1 ha). 

Media = Contenido de carbono promedio por cobertura vegetal. 

EE = Error estándar. 

Intervalo = Valores mínimos y máximos de contenido de carbono por cobertura vegetal. 
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DISCUSIÓN 

Estructura diamétrica del bosque 

Los bosques de la Magdalena Contreras 
son irregulares en cuanto a la edad (en la 
clasificación forestal corresponde al tér- 
mino multietáneo y homogéneo; Grijpma, 
2001), concentrándose la mayor cantidad 
de individuos en las etapas de regenera- 
ción inicial y avanzada (< 10 cm de DN; 
Fig. 5). 

Se observa que, derivado de las 
campañas de reforestación, se ha favore- 
cido a diferentes especies del género 
Pinus, cuya distribución natural no corres- 
ponde a la vegetación primaria de los bos- 
ques de la Magdalena Contreras. Los 
Pinus hartwegii son dominantes en la 
estructura diamétrica < 5 cm a 22,5 cm de 
DN, esto indica que se está desarrollando 
un proceso de sucesión dirigida y de com- 
petencia inducida que afecta a los rodales 
puros, tanto de Abies religiosa como a los 
de Pinus hartwegii. 

Según Pascual (2003), el Abies reli- 
giosa es una especie cuyo crecimiento 
está en función de la sombra, el suelo y 
la competencia con otros Abies sp. 
y especies del sotobosque. Dicha com- 
petencia puede provocar la mortandad 
de las plántulas de Abies religiosa, pues 
sus raíces son menos capaces que las 
de las angiospermas (sotobosque) de 
aprovechar el abastecimiento de agua y 
nutrimentos en el suelo, en sus etapas 
iniciales. 

Las hojosas representan 7% del 
total de la población de las plantas regis- 
tradas, con dominancia de arbustos en el 
sotobosque, que parecen ser exógenos y 
que compiten con las plántulas del 
estrato arbóreo y que según Nieto (op. 
cit.), también pueden favorecer el creci- 
miento de las plántulas de Abies reli- 
giosa, al generar un microclima, además. 


como son estacionales, les provee 
nutrientes para su crecimiento. Las plan- 
taciones de Cupressus sp. ocupan 3% 
del total de la población inmersa en las 
clases que fueron muestreadas y los 
Quercus sp. apenas registran 1%. Tam- 
bién se observa la presencia de rudera- 
les que indican perturbación (Acaena 
elongata). 

De acuerdo con los datos registra- 
dos, se cree que en un futuro los bosques 
de la cuenca del río Magdalena cambiarán 
su dominancia de Abies religiosa a Pinus 
sp. Posiblemente las condiciones del 
cambio de clima favorezcan esta transi- 
ción, ya que de haber un incremento en la 
temperatura, el género Pinus tendrá 
mayores probabilidades de sobrevivir, ya 
que presenta mejor adaptación al 
ambiente. 


Contenido de carbono 

La diferencia de valores de contenido de 
carbono almacenado, por tipo de cober- 
tura vegetal, presenta variaciones deriva- 
das de: 

La edad del arbolado, es decir, los 
bosques son multietáneos. El crecimiento 
diferencial derivado de la calidad de sitio 
refleja diferentes densidades de biomasa 
arbórea y por ende almacenes de carbono 
e incrementos disímiles en la misma 
cobertura vegetal. 

El estado de conservación de la 
masa forestal, la densidad del arbolado y 
las condiciones actuales de la estructura 
y composición de cada cobertura vegetal. 

La infraestructura de los caminos y 
las pendientes (accesibilidad). 

Otros procesos que afectan a la 
cuenca y por tanto al contenido de car- 
bono almacenado son; 
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Contaminación atmosférica. La ubi- 
cación espacial de la cuenca del Río 
Magdalena es el suroeste de la cuenca 
de México, dirección a la cual viajan los 
vientos dominantes, arrastrando la conta- 
minación generada en la ciudad, dando 
como consecuencia; clorosis, exposición 
a concentraciones de ozono, afectación 
de la fitomasa, alteración entre los ciclos 
reproductivos de diferentes especies, llu- 
via ácida, retirada de polinizadores y dis- 
persadores, entre otros problemas; esta 
contaminación también genera condicio- 
nes microclimáticas que favorecen el 
desarrollo y distribución vertical de pla- 
gas forestales. 

Plagas. La veda forestal en el suelo 
de conservación del D.F., impide la remo- 
ción de individuos que ya están enfermos 
o que son susceptibles de ser invadidos 
por organismos patógenos, propiciando el 
desarrollo de plagas, dando lugar a una 
disminución en la captura potencial de 
carbono, la purificación del aire, la genera- 
ción de oxígeno y un incremento en la dis- 
tribución de plagas (e.g.: muérdago y 
deseo rtezador). 

Incendios naturales e inducidos. De 
1995 a 2004 se han registrado 157 incen- 
dios (Flores, 2006) asociados a las 
sequías y a las actividades humanas; en 
1998 se registró uno de los incendios más 
severos, afectando la copa de los bos- 
ques de Abies religiosa, principalmente. 
Debido a la falta de comunicación de auto- 
ridades y de la comunidad Magdalena Atli- 
tic, no se han abierto los canales para 
darle un saneamiento adecuado a las 
zonas afectadas; esto también merma de 
forma directa el potencial de captura de 
carbono y sus almacenes, así como los 
continuos de bosque que permiten el des- 
plazamiento de las especies. 


Comparación con otros estudios de 
caso del área de la cuenca del río 
Magdalena 

Con el fin de conocer los resultados de los 
estudios relacionados con la estimación 
del contenido de carbono en el área de la 
Magdalena Contreras, se realizó un cua- 
dro comparativo (Tabla 5), analizando la 
superficie de la unidad de muestren, 
la superficie total muestreada y el carbono 
estimado: 

Se registraron tres estudios en el 
caso de Abies religiosa, en los que los 
autores indican valores que van de las 
58 tC/ha a las 207 tC/ha. Este estudio pre- 
senta resultados basados en la densidad 
del arbolado con valores que van de 
95 tC/ha a 162 tC/ha, lo que permite apre- 
ciar diferencias significativas respecto a 
las sobreestimaciones y subestimaciones 
en los almacenes de carbono en esta 
cobertura vegetal, aunado a la intensidad 
y tipo de muestren. 

Para el bosque de Pinus hartwegii, 
se comparan los resultados con dos estu- 
dios de caso, encontrando valores de 
carbono almacenado que van de 44 tC/ 
ha a 57,6 tC/ha, en ambos casos no se 
clasifica esta cobertura vegetal tomando 
como base la densidad del arbolado y la 
superficie que ocupa (CORENA, 2001 y 
Espinoza, 2005). 

Este estudio, separa esta clase con 
base en la densidad del arbolado, obser- 
vando una diferencia de casi el doble del 
carbono almacenado entre ambas; es 
por ello que se sugiere que las interpola- 
ciones y extrapolaciones del carbono 
almacenado por clase de cobertura vege- 
tal y uso de suelo, se deben hacer con 
conocimiento de causa. 
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Tabla 5. Comparación de los estudios de estimación de contenido de carbono en el suelo 
de conservación del sur del Distrito Federal y en la cuenca del río Magdalena. 


Cobertura vegetal 

Unidad 
de muestreo 
(ha) 

Superficie 

(ha) 

Carbono 

(tC/ha) 

Fuente 

Bosque de Ables religiosa 


0,1 

7,4 

144,07 

CORENA (2001) 


0,06 

11,64 

207,98 

Valenzuela (2001) 


0,0625 

1,75 

58 

Nava (2006) 

Densidad cerrada 

0,3 

6,6 

162,91 

Este estudio 

Densidad abierta 

0,3 

1.2 

95,49 

Este estudio 

Bosque de Pinus hartwegii 


0,1 

7,4 

57,61 

CORENA (2001) 


0,0625 

1,44 

44 

Espinoza (2005) 

Densidad cerrada 

0,3 

3 

47,37 

Este estudio 

Densidad abierta 

0,3 

1,5 

21,21 

Este estudio 


CONCLUSIONES 

En la cuenca del río Magdalena, la cober- 
tura vegetal que registró la mayor cantidad 
de carbono almacenado en la biomasa 
aérea fue Abies religiosa, en su densidad 
cerrada, lo que hace trascendental el diseño 
de políticas públicas, sobre todo en zonas 
cercanas a estos tipos de cobertura, dada 
su relevancia por los servicios ambientales 
que estos brindan a la Ciudad de México. 
Además son necesarias campañas de refo- 
restación adecuadas a las condiciones de 
composición actuales, ya que de seguir las 
tendencias de dichos programas en la 
cuenca, se originará una transición en la 
composición y estructura, principalmente, 
en los bosques de Abies religiosa, puesto 
que se registró que la especie más plantada 
en todas la coberturas vegetales son del 
genero Pinus, posiblemente derivado de la 
falta de información técnica-científica, o a su 
producción en los viveros cercanos. 

La variación en las densidades de 
carbono registradas por tipo de cobertura 


vegetal, obedecen a múltiples factores 
(e.g.: calidad de sitio, accesibilidad, cerca- 
nía a áreas urbanas, presión de los recur- 
sos, presencia de actividades no aptas 
para este tipo de vocación), pero el princi- 
pal factor de afectación y permanencia 
son las actividades antrópicas. Es nece- 
sario hacer comparaciones con estudios 
de caso in situ y ex situ, a fin de enrique- 
cer la posible toma de decisiones sobre el 
recurso, determinar las técnicas que expli- 
quen mejor el fenómeno, los alcances, las 
limitaciones y la precisión de los resulta- 
dos. 

El detalle fino del presente estudio 
permitió disminuir incertidumbres res- 
pecto a la distribución de la cobertura 
vegetal en la cuenca. La escala de análi- 
sis del presente estudio permite compara- 
ciones con los mapas temáticos generados 
por otras instituciones oficiales (e.g.: inegi, 
CORENA, Conabio, entre otros), ya que se 
puede realizar una generalización de las 
coberturas vegetales para hacerlas com- 
patibles. 
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Tarifas de biomasa aérea para abedul 
{Betula pubescens Ehrh.) y roble 
{Quercus robur L.) en el noroeste de España 

Above-ground biomass equations for birch 
(Betula pubescens Ehrh.) and pedunculate oak 
(Quercus robur L.) in north western Spain 

Esteban Gómez-García^*, Felipe Crecente-Campo^ 
y Ulises Diéguez-Aranda^ 

RESUMEN 

El objetivo de este trabajo fue el desarrollo de modelos de estimación de biomasa aérea para 
abedul {Betula pubescens Ehrh.) y roble {Quercus robur L.) en Galicia (noroeste de España). Para ello 
se emplearon datos obtenidos mediante el muestreo destructivo de 50 abedules y 50 robles, en los que 
se obtuvo el peso seco de biomasa total y por fracciones: madera del fuste, corteza del fuste, ramas 
mayores de 2 cm, ramas de 0,5 cm a 2 cm, ramas menores de 0,5 cm, y hojas. En un primer paso se 
seleccionaron los modelos que presentaban el mejor ajuste para cada fracción de biomasa arbórea 
considerada. Posteriormente, y para asegurar la aditividad, se realizó un ajuste simultáneo de las ecua- 
ciones de estimación de biomasa por fracciones junto con la de biomasa total, empleando el procedi- 
miento estadístico denominado NSUR {Nonlinear Seemingly Unrelated Regressions). También fue 
necesario un ajuste ponderado para corregir la existencia de heterocedasticidad. El número de condi- 
ción verificó que no existían problemas graves de multicolinealidad. Al final se obtuvo, para cada espe- 
cie, un sistema de siete ecuaciones de estimación de biomasa aérea para las distintas fracciones y 
para la biomasa total. Estas ecuaciones explicaron como mínimo 79% de la variabilidad observada, 
y en el caso de las ecuaciones de biomasa total 98% para abedul y 97% para roble. Se recomienda la 
utilización de las ecuaciones desarrolladas en este estudio en sustitución de las ecuaciones de bio- 
masa existentes para la región. 

PALABRAS CLAVE: 

Aditividad, ajuste simultáneo, Galicia, heterocedasticidad, multicolinealidad. 

ABSTRACT 

The aim of this study was to develop modeis for above-ground biomass estimation for birch 
{Betula pubescens Ehrh.) and pedunculate oak {Quercus robur L.) in Galicia (northwest of Spain). Data 
from a destructive sampling of 50 birch trees and 50 pedunculate oak trees, that provided total dry bio- 
mass weight and dry biomass weight by components: wood, bark, branches over 2 cm, branches 
between 0,5 cm and 2 cm, branches under 2 cm, and leaves, was used. Different modeis were fitted for 
each above-ground component and the best of them was chosen in a first step. All the modeis were 
refitted simultaneously by use of Nonlinear Seemingly Unrelated Regressions (NSUR) to ensure the 
additivity of the biomass equations. Weighted regression was necessary to correct the heterocedasticity 
of residuals. Not serious problems with multicollineality were observed. Finally, a System of seven equa- 
tions for the different above-ground biomass components, and for the total biomass, was obtained for 
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each species. These equations explained 
at least 79% of the observad variability. 
For total biomass the equations explained 
98% and 97% of the observad variability 
for birch and oak, respectively. The devel- 
oped biomass equations are recom- 
mended in substitution of the existing 
regional biomass equations. 

KEYWORDS: 

Addivity, simultaneóos fitting, Galicia, 
heterocedasticity, multicollineality. 


INTRODUCCIÓN 

La estimación de la biomasa forestal, así 
como la acumulación de carbono en la 
misma, ha ganado importancia en los últi- 
mos años como resultado del Protocolo 
de Kyoto (unfccc, 1997), que reconoce 
que los ecosistemas forestales pueden 
contribuir a la mitigación del efecto inver- 
nadero inducido por las actividades huma- 
nas de extracción y emisión de carbono 
(entre otros elementos) a la atmósfera. 
Por otra parte, la necesidad de reducir la 
dependencia de los combustibles fósiles 
ha traído como consecuencia la conside- 
ración de los residuos de las operaciones 
forestales como una fuente renovable de 
energía. Las mediciones necesarias para 
cuantificar la biomasa de un árbol o de un 
rodal con exactitud son lentas y costosas, 
por lo que se deben desarrollar ecuacio- 
nes que proporcionen estimaciones fia- 
bles de la biomasa y que sean aplicables 
en la gestión forestal práctica. Estas ecua- 
ciones relacionan la biomasa total o por 
fracciones con variables fáciles de medir 
en campo. 

Los primeros estudios sobre cuantifi- 
cación de biomasa arbórea surgieron para 
estimar la producción en peso seco de dife- 
rentes especies, destacando entre ellos los 
relacionados con Larix decidua Mili, y 
Picea abies (L.) Karst, en Suiza (Burger, 
1945; Burger, 1953). En las últimas déca- 
das han aumentado los estudios sobre la 
estimación de biomasa forestal, alcan- 


zando cada vez mayor relevancia y abar- 
cando un gran número de especies y 
estructuras de masa diferentes (p. ej., 
Pardé, 1980; Cunia, 1986; Zeide, 1987; 
Waring y Running, 1998; Johansson, 2000; 
Merino et al., 2003; Cienciala et al., 2005). 

A escala europea, Zianis etal. (2005) 
llevaron a cabo una revisión de las ecua- 
ciones de biomasa desarrolladas para 
especies forestales de los géneros Abies, 
Acer, AInus, Arbutus, Betula, Eucalyptus, 
Fagus, Fraxinus, Larix, Picea, Pinus, 
Popuius, Pseudotsuga, Quercus y Tilia. 
Estos autores encontraron 607 ecuacio- 
nes de biomasa, una minoría de ellas para 
el sur del continente. Montero et al. (2005) 
desarrollaron tarifas de biomasa de árbol 
individual para algunas de las principales 
especies forestales españolas (géneros 
Abies, AInus, Betula, Castanea, Erica, 
Ceratonia, Eucalyptus, Fagus, Fraxinus, 
//ex, Juniperus, Laurus, Myrica, Olea, 
Pinus, Popuius y Quercus), empleando 
modelos alométricos y el diámetro normal 
como variable independiente. En Galicia 
(nororeste de España, Fig. 1), Diéguez- 
Aranda et al. (2009) realizaron una revi- 
sión de las ecuaciones de biomasa 
desarrolladas para la región. Esta revisión 
incluía tarifas de biomasa de árbol indivi- 
dual para abedul {Betula pubescens 
Ehrh.) y roble {Quercus roburL.), desarro- 
lladas por Balboa-Muñas et al. (2005, 
2006, respectivamente); y también para 
Eucalyptus globulus Labill., E. nitens H. 
Deane & Maiden, Pinus radiata D. Don y 
P. p/nasfer Ait. En el caso del abedul usa- 
ron únicamente el diámetro normal como 
variable independiente, mientras que para 
roble se usó el diámetro normal y la altura 
total como variables independientes. El 
principal problema de las tarifas de bio- 
masa existentes para abedul y roble en 
Galicia radica en que se ajustaron 
empleando un reducido número de datos, 
recogidos únicamente en la provincia de 
Lugo y en zonas bastante próximas entre 
si, por lo que su representatividad a nivel 
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gallego es limitada. Además, el muestreo 
de biomasa fue realizado en invierno, por 
lo que los árboles muestreados no dispo- 
nían de hojas, y esta fracción la ajustaron 
con valores estimados. 


OBJETIVOS 

El objetivo de este estudio es el desarrollo 
de nuevas tarifas de biomasa de árbol 
individual para abedul y roble en Galicia, 
para lo que se dispuso de una base de 
datos más completa, que reflejaba mejor 
la realidad gallega y que contaba con 
observaciones en la fracción hojas para 
ambas especies. 

METODOLOGÍA 

Datos 

Para la elaboración de ecuaciones de 
estimación de biomasa aérea se mues- 


trearon destructivamente un total de 50 
abedules y 50 robles, procedentes tanto 
de regeneración natural como de planta- 
ción, repartidos en 11 y 9 localizaciones 
diferentes, respectivamente, en las pro- 
vincias de Lugo, Ourense y Pontevedra 
(Fig. 1). En el invierno de 2002-2003 se 
cortaron 16 abedules en 2 localizacio- 
nes y 31 robles en 4 localizaciones, 
estos datos son los que emplearon Bal- 
boa-Murias et al. (2005, 2006) en el 
desarrollo de sus tarifas de biomasa. En 
el verano de 2009 se cortaron 34 abedu- 
les en 9 localizaciones y 19 robles en 5 
localizaciones, con la intención de com- 
pletar la muestra existente y llegar a los 
50 árboles muestreados para cada 
especie. Se cortaron árboles de diferen- 
tes clases sociológicas (dominante, 
codominante, intermedio y sumergido) 
que no presentasen malformaciones en 
el tronco y la copa y sin desarrollos anó- 
malos o desequilibrios motivados por su 
proximidad a claros o bordes de masa o 
caminos. 



Figura 1 . Localización geográfica de Galicia y lugares donde se realizaron los muéstreos 
de biomasa para abedul (cuadrados) y roble (círculos). También se representan los 
bosques de la región donde dominan estas especies. 
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Antes de cortar cada árbol se midió su 
diámetro normal (a 1 ,3 m sobre el suelo) en 
dos direcciones (una siguiendo curvas de 
nivel y otra en dirección perpendicular a 
ésta) mediante forcípula (modelo Mantax 
Blue de 800 mm, Haglóf, Suecia; resolu- 
ción de 0,1 cm), y su altura total mediante 
hipsómetro (modelo Vertex IV, Haglóf, Sue- 
cia; resolución de 0,1 m). Tras cortar cada 
árbol se midieron la altura del tocón por la 
parte ladera arriba del árbol (siendo éste 
siempre inferior a 10 cm) y la longitud del 
tronco, ambas con cinta métrica. Posterior- 
mente, se dividió la biomasa arbórea aérea 
en las siguientes fracciones: fuste (hasta 7 
cm en punta delgada), ramas mayores de 
7 cm (diámetro mínimo en punta delgada 
de 7 cm), ramas gruesas (diámetro máximo 
en punta gruesa de 7 cm y mínimo en 
punta delgada de 2 cm) y ramas menores 
de 2 cm (incluyendo hojas). Estas fraccio- 
nes se pesaron en campo con una báscula 
GRAM modelo K3-300L (plataforma de 600 
mm X 600 mm, 300 kg de capacidad y 50 g 
de resolución). 

A continuación se recogieron mues- 
tras de cada fracción para llevar a labora- 
torio. En el fuste se cortaron tres discos 
transversales de 4 cm de grosor (el fuste 
se consideró dividido en tres trozas de 
igual longitud y se seleccionó una rodaja 
en el centro de cada una de ellas). Se 
recogió una muestra aleatoria de ramas 
mayores de 7 cm y una muestra aleatoria 
de ramas gruesas (aproximadamente 1 ,5 
kg en cada una de estas fracciones). En 
las ramas menores de 2 cm se tomó al 
azar una muestra de 6 kg (2 kg en la parte 
baja de la copa, 2 kg en la parte media y 2 
kg en la parte alta). Para evitar en lo posi- 
ble la pérdida de humedad, estas mues- 
tras fueron transportadas en sacos de 
plástico cerrados. 

Las muestras se procesaron en labo- 
ratorio el mismo día de la corta de cada 
pie. La muestra de ramas menores de 2 
cm se dividió en ramas finas (diámetro 


máximo en punta gruesa de 2 cm y mínimo 
en punta delgada de 0,5 cm), ramillos 
(diámetro máximo en punta gruesa de 0,5 
cm) y hojas (si las hubiese). Las rodajas 
del fuste se dividieron en madera y cor- 
teza. La proporción, en peso verde, de 
estas nuevas fracciones se aplicó al árbol 
completo, por lo que se consideraron siete 
fracciones: madera (del fuste), corteza (del 
fuste), ramas mayores de 7 cm, ramas grue- 
sas (2 cm - 7 cm), ramas finas (0,5 cm - 2 cm), 
ramillos (< 0,5 cm) y hojas. 

El contenido en humedad de las siete 
fracciones consideradas se determinó en 
laboratorio a partir de las muestras recogi- 
das en campo. Se empleó la totalidad de la 
cantidad disponible en cada fracción, salvo 
en las ramas finas, ramillos y hojas, consi- 
derando suficiente secar aproximadamente 
1,5 kg, 1,0 kg y 0,5 kg respectivamente. 
Las muestras fueron pesadas en verde con 
balanza de precisión gram modelo sbz- 
1500 (1500 g de capacidad y 0,01 g de 
resolución). Posteriormente se secaron en 
una estufa labolan-indelab modelo eco- 
LAN IDL.FI720 (720 I de Capacidad y dimen- 
siones interiores de 120 cm x 100 cm x 60 
cm) a 65 °C hasta peso constante (cuando 
dos mediciones realizadas con una sepa- 
ración de 12 horas difieren menos de 0,5% 
en peso). El contenido en humedad se cal- 
culó en función de los pesos verde y seco. 
A partir del contenido en humedad se 
obtuvo el peso seco total de cada fracción, 
que fue el utilizado para la elaboración de 
las ecuaciones de biomasa. 

El muestreo de biomasa del año 
2002-2003 fue realizado en invierno, por 
lo que los árboles muestreados no dispo- 
nían de hojas. Por tal motivo, el ajuste de 
la ecuación de estimación de biomasa 
para esta fracción se realizó únicamente a 
partir de los datos de los 34 abedules y de 
los 1 9 robles muestreados en el verano de 
2009. Con respecto a las fracciones de 
ramas finas y ramillos de abedul, los datos 
correspondientes al primer muestreo mos- 
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Figura 2: Distribución altura-diámetro de los árboles de B. pubescens (izquierda) 
y Q. robur (derecha) empleados en los ajustes. 


traron una tendencia anómala, siendo 
imposible investigar la causa, aunque 
podría ser motivado por deficiencias en el 
proceso de secado; por este motivo, el 
ajuste de las ecuaciones de estimación de 
biomasa de estas dos fracciones se rea- 
lizó a partir de los datos de los 34 abedu- 
les muestreados en 2009. En la figura 2 se 
representa, para cada especie, la distribu- 
ción altura-diámetro de los árboles 
empleados en los ajustes. En la tabla 1 se 
muestran, para cada especie, los estadís- 
ticos descriptivos de las variables diáme- 
tro normal, altura total y peso seco de 
biomasa aérea total y por fracciones. 


Ajuste de las ecuaciones de biomasa 

La metodología empleada en el ajuste de 
las ecuaciones de biomasa fue similar a la 
que utilizaron Balboa-Murias et al. (2006) 
para Q. robur en Galicia. Una vez determi- 
nada la biomasa seca de las diferentes 
fracciones de los árboles muestreados, se 
desarrollaron ecuaciones que relacionen 
dichos pesos con una o más variables 
independientes. Estas variables indepen- 
dientes deben ser: (i) poco numerosas y 
fáciles de medir para que el modelo tenga 
una amplia gama de aplicaciones y su 
empleo sea sencillo, (ii) variables fuerte- 


mente correlacionadas con la biomasa, y 
(iii) variables poco correlacionadas entre sí 
para evitar problemas de multicolinealidad. 

El diámetro normal es la variable 
más utilizada en las ecuaciones de esti- 
mación de biomasa arbórea {w = f{d)), 
aunque añadiendo la altura del árbol como 
segunda variable explicativa {w = f{d, h)) 
suelen mejorar significativamente los 
ajustes (Wang, 2006). Como variables 
independientes adicionales al diámetro y 
a la altura se pueden emplear otras varia- 
bles de árbol, como la longitud de la copa 
viva (diferencia entre la altura total y la 
altura a la que se localiza la primera rama 
viva) o el diámetro en la base de la copa, 
que han resultado especialmente útiles 
para mejorar las estimaciones de biomasa 
de las fracciones de copa (p. ej., Clark, 
1982; Satoo y Madgwick, 1982). También 
se pueden emplear variables de rodal 
como la edad, el área basimétrica, la 
altura dominante o la altura media (p. ej., 
Satoo y Madgwick, 1982; Bond-Lamberty 
et al., 2002; Porté et al., 2002). Estas 
opciones presentan como inconveniente 
que para emplear las ecuaciones de bio- 
masa es necesario un mayor esfuerzo de 
inventario y, además, las ecuaciones 
desarrolladas presentan un intervalo de 
aplicación más reducido. 


Tabla 1 . Estadísticos descriptivos de los datos de B. pubescens y Q. robar empleados en los ajustes. 
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La selección del modelo a emplear 
para la estimación de la biomasa de cada 
fracción ha de buscar un equilibrio entre la 
precisión, la sencillez y la aplicación prác- 
tica (Wang, 2006). Aunque existe una 
gran variedad de modelos válidos para 
estimar la biomasa, todos ellos derivan de 
alguna de las tres formas matemáticas 
siguientes (Pardé, 1980; Snowdon, 1985; 
Parresol, 1999): 

Lineal (error aditivo): 

7 = ySo + +... + i8yX^- + e (1) 

No lineal (error aditivo): 

Y = X* . . . X(‘ + £ (2) 

No lineal (error multiplicativo): 

Y = (3) 

donde; 

F: biomasa total o de una fracción determi- 
nada 

X:. variables explicativas 
¡3.: parámetros del modelo 
£• término del error 

El modelo (1) se ajusta por regresión 
lineal, el modelo (2) requiere emplear una 
técnica iterativa de regresión no lineal y el 
modelo (3) se ajusta por regresión lineal 
realizando una transformación logarítmica 
de ambos términos de la igualdad. 

En este trabajo, en un primer paso, se 
seleccionó para cada especie y para cada 
fracción de biomasa el modelo que presen- 
taba una mejor capacidad de ajuste. Los 
modelos evaluados emplean como varia- 
bles independientes únicamente variables 
de árbol individual (diámetro normal y, en 
ocasiones, altura total). Se ajustaron mode- 
los lineales y alométricos (modelos (1) y 
(2)). Los modelos lineales se ajustaron por 
mínimos cuadrados empleando el procedi- 


miento REG del programa sas/stat® (sas 
Instituto Inc., 2009) y los modelos no linea- 
les se ajustaron por mínimos cuadrados 
empleando el procedimiento nlin del 
mismo programa, usando el método itera- 
tivo de Gauss-Newton. En un segundo 
paso las ecuaciones de estimación de bio- 
masa por fracciones y de biomasa total se 
ajustaron simultáneamente para garantizar 
la aditividad. 

Inicialmente se ajustaron modelos 
para cada fracción muestreada, obtenién- 
dose resultados pocos satisfactorios para 
la fracción de ramas mayores de 7 cm en 
cualquiera de las dos especies, presen- 
tando además el ajuste simultáneo proble- 
mas de convergencia. Por ello, y teniendo 
en cuenta que no todos los pies disponían 
de datos de esta fracción, las ramas 
mayores de 7 cm se integraron en la frac- 
ción de ramas gruesas. Así las ramas 
gruesas pasaron a ser ramas con un diá- 
metro mayor de 2 cm. 


Heterocedasticidad y multicolinealidad 

Es frecuente que en el desarrollo de tari- 
fas de biomasa las estimaciones presen- 
ten problemas de heterocedasticidad, es 
decir que la varianza de los errores no sea 
constante (Parresol, 2001); y este fenó- 
meno invalida los contrastes de significa- 
ción de los parámetros. Para analizar la 
presencia de heterocedasticidad se anali- 
zaron los gráficos de residuos estudenti- 
zados (Mendoza et al., 2002) frente a 
valores estimados y se calculó el con- 
traste de White (1980), que es indepen- 
diente de la forma funcional del modelo y 
del tipo de distribución de datos. Para 
corregir la presencia de heterocedastici- 
dad se empleó regresión ponderada, 
como en otros trabajos de estimación de 
biomasa arbórea (p. ej., Schiaegel, 1982; 
Cunia, 1986; Parresol, 1999). Así, en la 
fase de ajuste, es necesario ponderar 
cada observación con un peso igual a la 
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inversa de la varianza de su error (a^J. 
Dicha varianza es desconocida, pero con 
frecuencia al examinar los gráficos de 
residuos frente a las variables indepen- 
dientes se observa que esta varianza se 
incrementa cuando se incrementa el valor 
de alguna o todas las variables indepen- 
dientes (X); por tanto se puede estimar 
utilizando una función potencial del tipo 
of=x^ (Neter et al., 1996). El valor del 
exponente k se puede determinar em- 
pleando la metodología propuesta por 
Park (1966), que consiste en utilizar los 
errores del modelo ajustado sin pesos e- 
como variable dependiente en el modelo 
potencial de varianza del error, es decir: 

ef = yX^ o bien (4) 

Inéf = \ny + k\nX. 

En función de las variables independien- 
tes que interviniesen en cada modelo se 
probaron distintos factores de pondera- 
ción (p. ej., etc.). 

Otro problema frecuente en el ajuste 
de ecuaciones de biomasa es la presencia 
de multicolinealidad (fuerte correlación 
entre las variables independientes). Sus 
efectos son similares a los indicados para 
la heterocedasticidad, invalidando tam- 
bién los contrastes de significación. 
Cuando los problemas de multicolineali- 
dad son graves el modelo no es robusto y 
la estimación de parámetros depende en 
exceso del conjunto particular de datos 
utilizado en el ajuste (Myers, 1990). Se 
puede evaluar la presencia de multicoli- 
nealidad en cada uno de los modelos 
ajustados calculando el número de condi- 
ción (NC): 


Cuando las variables definen un sis- 
tema ortogonal todos los autovalores son 
iguales a 1, mientras que cuando existe 
un alto grado de correlación entre varia- 
bles algún autovalor es próximo a cero. 
De acuerdo con Belsley (1991, pp. 139- 
141), si NC < 10 no hay problemas de 
multicolinealidad, si 30 < NC < 100 existen 
problemas de multicolinealidad y si 1000 < 
NC < 3000 los problemas son graves. 


Aditividad y ajuste simultáneo 

Una de las propiedades más importantes 
que deben cumplir las ecuaciones de esti- 
mación de biomasa de las distintas frac- 
ciones en que se divide un árbol es que la 
suma de sus estimaciones sea igual a la 
estimación obtenida por la ecuación de 
estimación de biomasa total, propiedad 
que se denomina aditividad. Aunque exis- 
ten varios procedimientos que se pueden 
utilizar para forzar esta propiedad en un 
sistema de ecuaciones (Parresol, 1999; 
Balboa-Murias et al., 2004), el ajuste 
simultáneo es el más flexible y general; 
puesto que permite usar diferentes mode- 
los matemáticos para cada fracción, con 
distintas variables independientes y fun- 
ciones de ponderación. Se basa en el 
ajuste de un sistema de ecuaciones, apa- 
rentemente no relacionadas, formado por 
las k fracciones arbóreas consideradas 
junto con la de biomasa total. 

^2 =72(^2) 

( 6 ) 


NC 


^ máximo autovalor matriz R '' 


mínimo autovalor matriz R 

\ / 


( 5 ) 


donde R es la matriz de correlaciones 
entre las variables independientes del 
modelo. 




) 
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donde; 

Wj^ ■ biomasa estimada 
: conjunto de variables explicativas 

En este caso, las variables independien- 
tes del modelo de biomasa total son todas 
las variables independientes que apare- 
cen en los restantes modelos, y la aditivi- 
dad está garantizada imponiendo 
restricciones a los parámetros del modelo 
de biomasa total. Este sistema, sin rela- 
ciones analíticas entre ecuaciones, se 
suele resolver empleando regresión sur 
{Seemingly Unrelated Regression) tam- 
bién conocida por mínimos cuadrados 
generalizados conjuntos o regresión Zell- 
ner (Zeliner, 1962). Para sistemas de 
ecuaciones no lineales este método se 
denomina nsur y es igual de válido que 
para ecuaciones lineales. Otra ventaja de 
esta metodología radica en corregir el pro- 
blema de dependencia de errores entre 
las estimaciones de cada fracción de bio- 
masa, considerando que no es realista 
asumir que los errores de cada ecuación 
no estén correlacionados (Borders, 1989; 
Parresol, 1999, 2001). 

Los parámetros de todas las ecua- 
ciones fueron obtenidos de forma simultá- 
nea utilizando la metodología nsur, para 
lo que se utilizó el procedimiento model 
del programa sas/ets® (sas Instituto Inc., 
2008). Se realizó un ajuste simultáneo ini- 
cial y a partir de los residuos obtenidos de 
éste se calcularon los pesos a utilizaren la 
ponderación del ajuste simultáneo defini- 
tivo, con la finalidad de corregir la hetero- 
cedasticidad y garantizar la aditividad. 


Ajuste en dos pasos 

Como ya se ha mencionado, en el caso 
del abedul algunos árboles no disponían 
de datos de las fracciones de ramas finas, 
ramillos y hojas (lo que en conjunto se 
denominó “copa fina”), debido a que no 


existían o a que fueron descartados por 
presentar anomalías. Sin embargo, el 
resto de fracciones presentaban tenden- 
cias lógicas, por lo que se decidió utilizar- 
las en los ajustes. Para solucionar este 
inconveniente, en un paso previo se rea- 
lizó un ajuste simultáneo de las ecuacio- 
nes de estas tres fracciones, más una 
cuarta que las engloba (biomasa de copa 
fina), a partir de los datos de los 34 árbo- 
les muestreados en el verano de 2009. 
Para obtener las estimaciones de los 
parámetros del resto de ecuaciones se 
realizó un ajuste simultáneo del sistema 
constituido por las ecuaciones de biomasa 
de fuste, de corteza, de ramas gruesas y 
de biomasa total (incluyendo en ésta la 
biomasa estimada de copa fina). 

En el caso del roble, al no existir la 
fracción de hojas para algunos árboles, se 
empleó un procedimiento similar. En pri- 
mer lugar se ajustó esta fracción de forma 
independiente, a partir de los 19 árboles 
muestreados en el verano de 2009. Poste- 
riormente se ajustaron el resto de ecua- 
ciones de forma simultánea (fuste, 
corteza, ramas gruesas, ramas finas, 
ramillos y biomasa total -incluyendo en 
ésta la biomasa estimada de las hojas). 


Criterios de selección de modelos 

La selección del mejor modelo para cada 
fracción se ha basado en el análisis de la 
capacidad de ajuste, tanto por métodos 
analíticos como gráficos. Así, a partir de 
los residuos, se calcularon dos estadísti- 
cos utilizados con frecuencia en modeliza- 
ción forestal: el coeficiente de 

determinación (R^) y la raíz del error medio 
cuadrático (remc); cuyas expresiones son 
las siguientes: 


= \- 




( 7 ) 
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( 8 ) 

n — p 


donde: 

Y- : valor observado de la variable depen- 
^ diente 

Y- : valor estimado de la variable depen- 
_ diente 

Y : valor promedio de la variable depen- 
diente 

n ; número total de observaciones utili- 
zado para ajustar el modelo 
p : número de parámetros a estimar 


Aunque existen opiniones fundamenta- 
das que plantean dudas en relación al 
empleo del en la selección de mode- 
los, este estadístico da una idea bastante 
intuitiva de la variabilidad que explican 
los modelos, y la utilidad general de 
alguna medida global de la adecuación 
del modelo parece anular alguna de 
estas limitaciones (Ryan, 1997, p. 424); 
no obstante no debe utilizarse como 
único criterio para elegir el mejor modelo 
(Myers, 1990, p. 166). El estadístico 
REMC resulta útil porque está expresado 
en las mismas unidades que la variable 
dependiente, por lo que da una idea del 
error medio que se comete con el modelo; 
además penaliza los modelos con mayor 
número de parámetros, de acuerdo con 
el principio general de simplicidad cientí- 
fica (Peña, 2002, p. 570). 

Además del cálculo de estadísticos 
de ajuste, una herramienta muy útil para 
determinar si los ajustes se adecúan a los 
datos son los análisis gráficos (Huang, 
2002). Se realizó una representación grá- 
fica de los valores reales de las variables 
frente a los predichos, a los que se super- 
puso la recta ideal 1:1, con el fin de ins- 
peccionar de manera visual si las 


predicciones eran acordes con los valores 
reales. Por último, se analizaron también 
los gráficos de residuos estudentizados 
frente a valores predichos con la finalidad 
de detectar valores atípicos, tendencias 
extrañas de los datos o la existencia de 
problemas de heterocedasticidad. 

RESULTADOS Y DISCUSIÓN 

Los modelos que mostraron los mejores 
resultados para la estimación de la bio- 
masa aérea total y por fracciones para 
abedul y roble se muestran en la tabla 2. 
En la tabla 3 se muestran las estimacio- 
nes de los parámetros obtenidas mediante 
el ajuste simultáneo ponderado, así como 
sus correspondientes errores estándar 
aproximados y contrastes de significación. 
Todos los parámetros resultaron significa- 
tivos al nivel de 5%, salvo el parámetro 
en roble, aunque de forma visual se com- 
probó la idoneidad del modelo elegido y la 
necesidad del parámetro en cuestión, por 
lo que éste se mantuvo en el modelo. 
También se muestran los pesos utilizados 
en la ponderación, los estadísticos de 
bondad de ajuste y el valor del contraste 
de White, que indica que en todos los 
casos el ajuste se realizó en condiciones 
de homogeneidad de la varianza de los 
errores. El análisis gráfico de los residuos 
estudentizados (no mostrado) corroboró 
estos resultados. 

En general, las ecuaciones mostraron 
un ajuste satisfactorio, explicando como 
mínimo 79% de la variabilidad observada. 
En el caso del abedul, cinco de las siete 
ecuaciones desarrolladas explicaron más 
de 90% de la variabilidad de los datos 
observados. Con respecto al roble, sólo 
dos ecuaciones superaron este nivel. Cabe 
destacar, en ambos casos, el buen ajuste 
que presentan las ecuaciones de biomasa 
de madera de fuste y de biomasa total, con 
valores del en torno a 95%-98% para 
ambas especies. 
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Tabla 2. Modelos seleccionados para la estimación de la biomasa aérea total y por 
fracciones para 8. pubescensy Q. robaren Galicia. 


Especie 

Fracción 

Modeio 


Madera fuste 

- Ao + A.í/'/í 


Corteza fuste 

u. //^ 


Ramas gruesas 

"V.: - A,</" ■ 

B. pubescens 

Ramas finas 



Ramillos 



Hojas 

- AmíA// 


Copa fina 



Total 

- Alo + + A>i‘A‘ //" " + 


Madera fuste 

u„=/í,„ + //,,cA/; 


Corteza fuste 

m; = /y,,cA/; 


Ramas gruesas 

»‘V>> - 

Q. robur 

Ramas finas 



Ramillos 

"’ríO:! ” Avi Am^^ 


Hojas 



Total 

“ Aio A?o ^Aii A:i A4i^ ^ 


= peso seco de las distintas fracciones o total de biomasa arbórea (kg) 
¡}.j= parámetro y de la fracción arbórea i 
d = diámetro normal del árbol (cm) 
h = altura total del árbol (m) 


Tabla 3. Estimaciones de los parámetros, contrastes de significación aproximados, estadísticos de bondad del ajuste y pesos utilizados 
en la ponderación para las ecuaciones de biomasa de cada fracción arbórea y de biomasa total obtenidas mediante ajuste simultáneo 

NSUR (ver modelos en la tabla 2). 
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parámetro y de la fracción i (madera, corteza, ramas gruesas, ramas finas, ramillos y hojas). 
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Sin embargo, para el caso del roble, 
la REMC fue superior siempre, aunque los 
valores no son directamente comparables, 
puesto que los valores medios de biomasa 
fueron mayores en roble que en abedul 
(Tabla 1). Considerando las mismas frac- 
ciones, los estadísticos de ajuste de las 
ecuaciones desarrolladas en este estudio 
fueron superiores (mayores valores de y 
menores valores de remc) a los obtenidos 
en estudios previos para las mismas espe- 
cies en la región (Balboa-Murias et al., 
2005; Balboa-Murias et al., 2006), con la 
excepción de la fracción hojas en roble. En 
esta fracción el valor de fue mayor en el 
estudio de Balboa-Murias et al. (2006) -R^ 
= 0,832 frente al valor de 0,789 obtenido en 
este estudio -, pero a su vez también el 
valor de remc fue mayor (11,2 kg frente al 
valor de 2,53 kg obtenido en este estudio). 
Además, también se debe considerar que 
en los estudios previos la fracción hojas se 
ajustó con valores estimados. 

En relación con la biomasa de las 
fracciones de copa (ramas gruesas, ramas 
finas, ramillos y hojas), éstas presentaron 
un peor ajuste que la biomasa del fuste. El 
mismo resultado ha sido constatado en 
estudios similares, especialmente en el 
caso de frondosas (Satoo y Madgwick, 
1982), lo que posiblemente sea debido a la 
variabilidad de la estructura de la copa, el 
número de ramas y la variación de la den- 
sidad de la madera a lo largo de las ramas 
(Fardé, 1980). La inclusión de otra varia- 
ble, como la longitud de copa viva, podría 
mejorar el ajuste de estas fracciones, aun- 
que a expensas de un mayor esfuerzo de 
muestreo para poder aplicar los modelos. 

El ajuste simultáneo mediante nsur 
proporcionó estimaciones de los paráme- 
tros y estadísticos de bondad del ajuste 
muy similares a los obtenidos mediante el 
ajuste individual; sin embargo, los errores 
estándar aproximados disminuyeron sen- 
siblemente cuando se empleó nsur (resul- 
tado directo de considerar la correlación 


entre los errores de las ecuaciones de bio- 
masa), por lo que se logró una mayor efi- 
ciencia en la estimación de parámetros 
con esta metodología. En otros estudios 
se obtuvieron resultados similares, verifi- 
cándose también que el método de esti- 
mación NSUR mejora la precisión en las 
estimaciones (Parresol, 1999; Parresol, 
2001 ; Carvalho y Parresol, 2003; Návar et 
al., 2004). El número de condición en el 
ajuste simultáneo ponderado fue de 41 ,2 
para abedul y 16,7 para roble, con lo que 
se puede concluir que no se presentan 
problemas graves de multicolinealidad 
entre las variables independientes. 

En las figuras 3 y 4 se muestran los 
gráficos de valores observados frente a pre- 
dichos para las distintas fracciones de bio- 
masa consideradas y para la biomasa total. 
Las nubes de puntos se sitúan cerca de la 
diagonal 1:1, indicando que los modelos 
proporcionan estimaciones con bajo sesgo. 
En el caso del roble, y para la biomasa de 
ramas gruesas, aparece una ligera tenden- 
cia a la subestimación para los mayores 
valores. Esta ecuación es una de las que 
presenta mayor valor de remc (115,9 kg). 

La proporción media de biomasa de 
las fracciones madera del fuste, corteza 
del fuste, ramas (ramas gruesas, ramas 
finas y ramillos) y hojas con respecto a la 
biomasa aérea total fue, respectivamente, 
de 53,3%; 8,9%; 34,2% y 3,6% en abedul 
y de 57,1%; 9,1%; 31,2% y 2,6% en roble. 
Resultados similares fueron obtenidos en 
estudios previos en Galicia (Balboa- 
Murias et al., 2005; Balboa-Murias et al., 
2006). Ambas especies presentan porcen- 
tajes similares de las mismas fracciones, 
con una proporción de las fracciones que 
forman la copa entre 34% y 38%, mayor 
que otras especies para las que se han 
desarrollado ecuaciones de estimación de 
biomasa en la región (E. globulus, E. 
nitens, P. radlata y P. pinaster) con una 
proporción de biomasa de copa entre 10% 
y 19% (Diéguez Aranda et al., 2009). 


wtotai observado {kg) Wr<os observado (kg) Wr>2 observado (kg) Wm observado (kg) 
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Wr<05 predicho (kg) 






iVft predicho (kg) 


Figura 3. Valores observados frente a predichos para las distintas fracciones de biomasa 
y para la biomasa total de 8. pubescens (variables definidas en Tabla 2). 

Los gráficos de biomasa de ramas finas, ramillos, hojas y biomasa total sólo incluyen 
datos de los 34 árboles muestreados en verano de 2009. 
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vji, predicho (kg) 



Figura 4. Valores observados frente a predichos para las distintas fracciones de biomasa 
y para la biomasa total de Q. robur (variables definidas en Tabla 2). 

Los gráficos de biomasa de hojas y biomasa total sólo incluyen datos de los 19 árboles 

muestreados en verano de 2009. 
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En abedul las fracciones se agrupan 
en una proporción de biomasa de fuste 
de 62,2% y de biomasa de copa del 
37,8%. Estos resultados son acordes con 
los obtenidos en otros estudios para 
especies del mismo género en España 
(Gracia et al., 2002) y Norteamérica 
(Jokela et al., 1981; Schmitt y Grigal, 
1981; Maliondo et al., 1990; Jian et al., 
1998), pero difieren de los obtenidos por 
Markiund (1988) para Suecia, en donde 
no se consideraba la fracción hojas y las 
proporciones eran 79,9% de biomasa de 
fuste y 20,3% de biomasa de ramas; en 
este estudio, si no se considera la frac- 
ción hojas las proporciones son 64,5% de 
biomasa de fuste y 35,5% de biomasa de 
ramas. 

En el caso del roble, las fracciones 
se agrupan en una proporción de biomasa 
de fuste de 66,2% y de biomasa de copa 
de 33,8%. Estos resultados son acordes 
con los obtenidos en otros estudios para 
especies del mismo género en España 
(Rapp et al., 1999; Gracia et al., 2002) y 
Portugal (Carvalho y Parresol, 2003), pero 
difieren de los obtenidos en estudios reali- 
zados en Centro Europa para Quercus 
petraea (André y Ponette, 2003) y en Nor- 
teamérica para Q. alba L., Q. coccínea 
Muenchh., Q. prinus Willd. y Q. rubra L. 
(Martin etai, 1998), en los que la biomasa 
de fuste representaba 73% a 84% de la 
biomasa aérea total. 

En este estudio no se han desarro- 
llado ecuaciones de estimación de bio- 
masa de raíces debido al enorme esfuerzo 
de muestreo que supone la toma de datos. 
Ecuaciones de este tipo fueron desarrolla- 
das para B. pubescens por Balboa-Muñas 
et al. (2005) y para Q. robur por Balboa- 
Muñas et al. (2006). Pese a que el número 
de observaciones utilizado fue muy redu- 
cido (6 datos para abedul y 11 para roble), 
es la única información disponible por el 
momento. 


CONCLUSIONES 

Las ecuaciones desarrolladas en este 
estudio permiten estimar con suficiente 
precisión los pesos secos de las distintas 
fracciones consideradas y la biomasa 
aérea total. Éstas podrán ser empleadas 
con varios fines, como la cuantificación 
del potencial de las dos especies para el 
almacenamiento de carbono a medio y 
largo plazo en los montes de Galicia, o el 
posible aprovechamiento energético de 
los residuos (fundamentalmente los 
extraídos de la copa) que se generan en 
las operaciones silvícolas intermedias 
(aclareos, claras y podas) y en el aprove- 
chamiento final. 

El conjunto de datos utilizado en este 
estudio -significativamente superior al 
que existía hasta el momento- y las mejo- 
ras metodológicas en el caso de abedul 
podrían explicar los mejores resultados 
obtenidos en la fase de ajuste con res- 
pecto a estudios previos para estas espe- 
cies en Galicia. Por lo tanto, se recomienda 
la utilización de las ecuaciones desarrolla- 
das en este estudio en sustitución de las 
ecuaciones existentes. 
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Dendrochronological reconstruction 
of long-term precipitation patterns in 
Basaseachi National Park, Chihuahua, México 

Reconstrucción dendrocronológica de patrones 
de precipitación de largo plazo en el Parque Nacional 
de Basaseachi, Chihuahua, México 

Christen M. Irby\ PeterZ. Fulé^ *, Larissa L. Yocom^ 

y José Villanueva Díaz^ 

ABSTRACT 

The purpose of this study was to reconstruct the precipitation at Basaseachi National Park (BNP) in 
Chihuahua, Mexico.Tree-ring samples from pine species including Pinus durangensis, P. lumholtzii and 
P.engelmannii were collected in and near BNP and they were cross-dated with existing chronologies. 

Ring widths of each sample were measured and modeis to remove non-climatic trends from the data were 
applied. The relationship between precipitation from the weather station in Yécora, Sonora, and the ring 
width indices were modeled using a bivariate linear regression (r^ = 0,59, p<0.01). It was found that preci- 
pitation from the months of October, December, February and April was most highiy correlated with tree 
ring growth. The reconstruction extends 225 years from 2007-1782. During this time period there were 
thirteen severe droughts, ten of which affected large areas in northern and central México. 

KEYWORDS: 

Climate, dendrochronology, drought, Pinus, precipitation reconstruction. 


RESUMEN 

El propósito de este estudio fue generar una cronología de anillos de crecimiento lo más extensa 
posible y usarla como un método indirecto para el desarrollo de una reconstrucción de precipitación 
estacional en el Parque Nacional de Basaseachi (BNP), Chihuahua, México. Para el desarrollo de la 
cronología, se colectaron núcleos de crecimiento de ejemplares de tres especies de pino, ubicados 
tanto dentro como fuera del BNP. Las muestras se fecharon con cronologías existentes y cercanas al 
área de estudio. Se midió el ancho de anillos de cada una de las muestras y se aplicó un procedimiento 
de estandarización para remover tendencias biológicas y geométricas. Con los índices dendrocronoló- 
gicos y datos de precipitación de la estación climática Yécora, Sonora, se generó un modelo de regre- 
sión lineal bivariado que tuvo una varianza significativa (r^ = 0,59, p<0,01). La precipitación de los 
meses de octubre, diciembre, febrero y abril fue la más correlacionada con el grosor de los anillos. La 
reconstrucción se extendió 225 años, período 2007-1782. Durante este tiempo ocurrieron trece sequías 
severas, diez de las cuales afectaron grandes áreas en el norte y centro de México. 
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Clima, dendrocronología, Pinus, sequía, reconstrucción de precipitación. 
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INTRODUCTION 

Megadiverse ecosystems are those in 
which the magnitude of species richness 
is one to two orders higher than that of 
most températe zones (Brehm et al., 
2008; McCune, 1988). Many of these 
ecosystems are being damaged by 
human land use practices such as fire 
suppression, logging and livestock gra- 
zing. Additionally, global warming is 
currently affecting the climate in these 
regions, adding to the disturbance within 
the ecosystems and this trend is expec- 
ted to become more pronounced in the 
near future (Rosenzweig et al., 2007). 
From a conservation perspectiva, mega- 
diverse ecosystems require an increased 
level of understanding about climatic 
variability over time as well as climate’s 
influence on the hydrologic cycle and the 
growth of dominant species (Villanueva- 
Díaz et a/., 2009). 

The Sierra Madre Occidental (smo) 
in Northwestern México is an example of 
a megadiverse system (Bye, 1995). The 
SMO holds some of the most extensiva 
and diverse pine-oak forests remaining in 
the World, yet there is relatively littie pro- 
tection of these ecosystems (Bye, 1995; 
HeyerdahI and Alvarado, 2003). The 
aquifers and surface runoff provided by 
the SMO sustain life for over six million 
people living in the States of Sonora and 
Chihuahua. Basaseachi National Park 
(bnp) is of interest for conservation 
because it contains the second tallest 
waterfall in México (246 m) and draws 
thousands of tourists annually, stimula- 
ting the local economy and providing pro- 
tection for a small area of forest on the 
canyon rim. 

The study of historical climatic 
variability is important for understanding 
current and future climatic trends as well 
as their relationship with social and eco- 
nomic stability (Timmerman et al., 1999; 


Villanueva-Díaz et al., 2007). Instrumen- 
tal weather data in México is short and 
incomplete, thus it is necessary to obtain 
information about past climate from a 
proxy source. Climate reconstructions 
provide useful information on the range 
and variability of precipitation, which can 
aid in resource management, restora- 
tion, and conservation efforts (Villa- 
nueva-Díaz et al., 2007). Several 
tree-ring chronologies exist for nor- 
thwestern México, compiled by Stahie 
and Cleveland (1993), Villanueva-Díaz 
and McPherson (1996), Biondi (2001), 
Díaz et al. (2002), Cleveland et al. 
(2003), González-Elizondo et al. (2005), 
Villanueva-Díaz et al. (2009); Cerano- 
Paredes et al. (2009). 

These chronologies range from 320 
to 607 years long. Most of the existing 
chronologies in the smo were constructed 
from Pseudotsuga menziesii because 
trees of this species are long-lived and 
sensitiva to climate. However, in northern 
México, this species is restricted to the 
highest elevations, northern slopes, and 
canyon bottoms. Pine species, in con- 
trast, are widespread. By using pine spe- 
cies including Pinus durangensis, P. 
lumholtzii and P engelmannii, it is possi- 
ble to expand the number of areas availa- 
ble for climate reconstruction as well as 
provide a better regional understanding 
of past climate. 


OBJECTIVES 

The objectives of the present study were 
to (1) develop the oidest possible chrono- 
logy using tree-ring records from pine 
trees in and around Basaseachi National 
Park (bnp), (2) model precipitation to deve- 
lop long-term climate reconstructions from 
the proxy tree-ring data, and (3) evalúate 
historical patterns of drought and its rela- 
tionship to the El Niño Southern Cscilla- 
tion (enso). 
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METHODS 

Study Site 

The study site is located in a pine-oak 
forest in and around BNP in Chihuahua, 
México (Fig. 1 y 2), centered at about 
108,2° West, 28,2° North (Fig. 3). Eleva- 
tion ranges from approximately 1720 to 
2300 m. The climate can be characterized 
as arid to semi-arid with summer monsoo- 
nal rains. Average annual precipitation in 
this area is approximately 700 mm (eric II, 
2000). The majority of the soil is classified 
into Leptosols and Regosols. Species pre- 
sent in this región inelude Pseudotsuga 
menziesii (Shwerin), Pinos durangensis 
(Martínez), P lumholtzii (Robinson and 
Fernaid), and P. engelmannii (Carriére) 
(Spellenberg et al. (1996), species author- 
ship from Farjon and Styies (1997) and 
Uchytil (1991)). Much of the forest in the 
Sierra Madre Occidental has been heavily 
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logged (Cortés Montaño and Cruz Gais- 
tardo, in press). Logging started in the 
area in the late 1800’s and proceeded 
through the 20'^ century and as a result, 
oniy a small fraction of the original vegeta- 
tion remains intact (Lammertink, 1996). 
Overgrazing has aiso been an issue in the 
area (Valero et al., 2001). Frequent, low 
intensity fires are important disturbance 
agents in these forests, which burned 
every 4-5 years before fire exclusión 
began in the 1940’s (HeyerdahI and Alva- 
rado, 2003 and Fulé etal., 2005). 


Sample Collection 

Tree species used in this chronology 
inelude Pinos dorangensis, P. lomholtzii 
and P. engelmannii. Cores from large 
mature trees growing in climatically sensi- 
tiva areas characterized by shallow, rocky 
soils and steep terrain were extractad 



Figure 1. Pictures taken at the study site, Basaseachi National Park (BNP), Chihuahua. 
Top- looking out from Cascada Basaseachi, samples were taken on the rim as well as the 

canyon bottom. 
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Figure 2. Pinus lumholtzii. 


(Villanueva-Díaz etal., 2007). Using incre- 
ment borers, two cores were removed 
from each suitable tree. The cores were 
taken 40 cm from the g round. To extend 
the length of the chronology, several 
cross-sections from oíd trees collected 
during a companion fire history recons- 
truction study in the same area were used 
(Yocom et al., unpublished data). All wood 
samples were air-dried, cores stored in 
paper straws, and cross sections wrapped 
in plástic wrap for transport. 


Laboratory Analysis 

All cross sections and cores were glued to 
wood mounts and sanded with increasin- 
gly finer grades of sandpaper until cells 
were clearly visible under a microscope 
(Grissino-Mayer, 2001). Using existing 
chronologies, the samples were visually 
crossdated and with skeleton plots (Fritts, 


1976; Stokes and Smiley, 1968). Tree-ring 
widths were then measured to the nearest 
0,001 mm using a stereo-zoom micros- 
cope and a Velmex sliding-stage microme- 
ter. To check the cross dating accuracy, 
the program cofecha was used (Holmes, 
1986; Grissino-Mayer, 2001). Using the 
program arstan, non-climatictrends, were 
removed including the general decrease in 
ring width as the tree ages, and autore- 
gressive growth persistence, by appiying a 
negativa exponential growth curve (any k) 
followed by autoregressive modeling. The 
residual chronology was used, which 
represente the maximal interannual varia- 
bility after removing the influence of auto- 
correlation, for the final reconstruction. 

To determine the correlation between 
monthly precipitation and annual tree-ring 
growth, residual ring-width indices were 
compared with monthly precipitation and 
temperatura data using the program Den- 



Madera y Bosques 19(1), 2013:93-105 


97 



1:100,000 

Contour interval 50 m 


Figure 3. Sampling sites in and around 

BNP. 

droClim 2000 (Biondi and Waikul, 2004; eric 
II, 2000). This was a resampling p- valué. 
The highest association of tree ring indices 
and precipitation was Oct-Apr, thus precipi- 
tation from these months was used for the 
final reconstruction. In selecting a weather 
station precipitation data from three different 
weather stations were compared with similar 


elevations and in relatively cióse proximity 
to BNP. There is a weather station, at the 
park but the data were too short and incom- 
plete to be used. The weather station at 
Yécora, Sonora was chosen because the 
association between ring width indices and 
precipitation was the highest of all the sta- 
tions that were considered and the station 
had a longer, more complete data set than 
the others. 

A least-squares regression was cal- 
culated between meteorological records of 
October-April precipitation and the resi- 
dual ring-width chronology to find the 
correlation of precipitation with tree-ring 
width (Fritts, 1976). The precipitation data 
and tree ring indices were spiit into two 
periods: a calibration period and a verifica- 
tion period. The weather data gathered at 
the Yécora station spanned the time 
period between 1971 and 2004. The data 
were spiit into two periods: 1971-1987 and 
1988-2004. Each period served as a cali- 
bration and validation period, in succes- 
sive analyses. The alpha level was 0,05. 
The regression equation developed in the 
calibration period was applied to the ring- 
width indices of the verification period. The 
reduction of error (re) was calculated as 
well as the coefficient of efficiency (ce) for 
the model. These two statistics are used to 
compare estimated data with valúes assu- 
med to be equal to the calibration period 
mean (Briffa et al., 1988). Positiva valúes 
indícate reasonable skill in the reconstruc- 
tions (Fritts, 1976). The ce differs from the 
RE in that it is a measure of the common 
variance between the actual and estima- 
ted data over the verification period whe- 
reas the latter is a measure of variance 
over the calibration period (Briffa et al., 
1988). The average expressed population 
signal (eps) was aiso calculated. The eps 
is used to express the degree to which a 
particular sample chronology portrays a 
hypothetically perfect chronology with hig- 
her valúes being more favorable (Wigley 
et al., 1984). 
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To create the final reconstruction, all 
the data were combined and a least- 
squares regression calculated between 
October-April precipitation recorded for 
the entire period of the weather station 
record and the ring-width indices. Past 
precipitation was estimated by appiying 
the model to the full data set back to the 
year 1782. 

The effect of enso on precipitation 
was examined by correlating the ñiños 
data set with the ring-width Índex. To 
determine major drought periods, a three 
year running average was calculated 
and identified periods where precipita- 
tion was less than 86% of the 225 year 
average. 


RESULTS AND DISCUSSION 
Chronology 

Pine trees from Basaeachi National Park 
proved to be useful for developing a multi- 
century precipitation reconstruction. The 
final chronology contained 44 trees and 
the series intercorrelation was 0,61. The 
climate reconstruction was extended back 
225 years, covering the time period bet- 
ween 1782 and 2007. The eps for the 
series was 0,87. For this chronology, the 
average eps exceeds the critical level of 
0,85 determined by Wigley et al. (1984). 


Precipitation Reconstruction 

Precipitation occurring in the previous win- 
ter, during the months of December and 
February, was most highiy correlated to 
tree ring growth, with valúes of 0,43 and 
0,50 respectively. In addition, October and 
April precipitation were found to have a 
positiva correlation with tree growth at 
0,27 and 0,24 respectively. Because these 
months had the highest correlations, total 
precipitation occurring between the 


months of Oct and April was used for the 
final reconstruction. There were no statis- 
tically significant correlations between tree 
growth and temperatura. Although they 
were not significant, the valúes suggest a 
negativa correlation between temperatura 
and annual tree-ring growth (Biondi and 
Waikul, 2004). 

The calibration period (1971-1987) 
had a statistical significance of r= 0,79 
(r^=0,63) and the verification period (1988- 
2004) had a statistical significance of 
r=0,74 (r^= 0,55) (Fig. 4). The re derivad 
from the model was 0,67, and the ce was 
0,63 indicating high skill in the reconstruc- 
tion (Table 3). The re and ce calculated 
with reversal of the calibration and verifi- 
cation periods were 0,63 and 0,55 respec- 
tively again indicating high skill in the 
reconstruction (Table 3). 

The linear regression representing 
the total time period of weather data 
(1971- 2004) and corresponding tree ring 
indices had a correlation of r= 0,77 
(r^= 0,59), meaning that approximately 
59% of the variability in tree ring growth 
was explained by precipitation data. The 
linear regression used to calcúlate precipi- 
tation valúes for the total reconstruction 
was: Y = -244,85 + 499, 03X where Y is the 
amount of October-April precipitation (in 
mm) and X is annual tree-ring width Índex. 
The model parameters were highiy signifi- 
cant (Tables 1 and 2). 


Effects of ENSO 

ENSO is known to have a large effect on pre- 
cipitation in northern México (Cleveland et 
al., 2003; Stahie and Cleveland, 1993). 
Generally in this región. El Niño is associated 
with higherwinter precipitation while its coun- 
terpart La Niña is associated with winters 
producing less precipitation (Magaña et al., 
2002). No significant correlation was found 
upon comparing enso patterns with winter 
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Table 1 . Linear regression model used to calcúlate the total reconstruction model is 

Y=aQ+a^X. 


Period 

Coefficients 

t ratio 

Prob>\t\ 

Standard Error 


a = -244,85 

-3,38 

0,0019 

72,3374 

1971-2004 




0,59 


a.i=499,03 

6,72 

<0,0001 

74,2628 


Table 2. Analysis of Variance of model used for the reconstruction. 


Source 

DF 

Sum of 
Squares 

Mean 

Square 

F Ratio 

Prob > F 

Model 

1 

403998 

403999 

45,1556 

<0,0001* 

Error 

32 

286298 

8947 



Total 

33 

690296 





Table 3. 

Reduction of error, coefficient of efficiency and final model. 


Calibration 

Calibration 



(1971-1987 

(1988-2004 

Final model 


Verification 

Verification 



(1988-2004 

(1971-1987 


Calibration 

0,63 

0,56 

R2 = 0,59 

Verification 

0,56 

0,63 


Reduction of error 

0,67 

0,63 

y= 244,85+499,03x 


Coefficient of 

0,63 

0,55 

efficient 
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precipitation but enso extreme years tend to 
match up with periods of uncharacteristically 
high or low precipitation, which is consistent 
with the study by Stahie and Cleaveland 
(1993). From 1782 to 1970, there were 13 
years classified as La Niña extreme events 
(Stalhe and Cleaveland, 1993). Reconstruc- 
ted precipitation for each of these 1 3 years 
was less than or equal to 90% of average 
precipitation (Table 4). There were nine El 
Niño events during this same time period. Of 
the nine years, reconstructed precipitation 
was higher than average for all but two years: 
1844 and 1964 (Table 4). 


Drought 

This reconstruction identified three drought 
periods not present in similar studies: 
1902-1904, 1943-1946 and 1959-1967. 
One likely explanation for the appearance 
of these droughts at bnp is the lower eleva- 
tion relative to the elevations of previous 
study sites. The identification of unique 
drought periods at bnp further demonstra- 
tes this area’s sensitivity to climate patterns 


and can aid resource managers in perser- 
ving the natural resources in the area. The 
five driest years identified by this recons- 
truction in order of increasing severity were 
1909, 1982, 1798, 1887 and 2000. The 
total reconstruction identified thirteen 
droughts lasting more than three years in 
duration occurring during the past 225 
years (Fig. 5). Severe droughts identified 
by the reconstruction are listed in Table 5. 
The most prolonged drought lasted 16 
years and occurred between the years of 
1 835 and 1 851 . Most of the drought periods 
identified by the reconstruction are aiso 
present in other reconstructions conducted 
in northern México with the exception of the 
three periods listed above (Villanueva et 
al., 2007). Many of these drought periods 
can be linked to social and economic hard- 
ship in México. The first of was the period 
between 1785 and 1786, although the 
reconstruction obtained here oniy shows 
1785 to be below average. This is known 
as “El Año del Hambre” (The Year of Hun- 
ger) as it resultad in extreme food shorta- 
ges and an outbreak in epidemic disease 
throughout México (Therrell, 2005). 


Table 4. List of enso extreme years from Stahie and Cleaveland (1993) and percent of 
average precipitation from our model that fell during these years at the study site. 


La Niña Extremes 

El Niño Extremes 

Year 

% precip (of ave) 

Year 

% precip (of ave) 

1789 

45 

1792 

126 

1801 

76 

1844 

73 

1805 

52 

1869 

108 

1855 

76 

1881 

101 

1862 

84 

1919 

121 

1887 

0 * 

1931 

195 

1904 

46 

1944 

113 

1909 

13 

1958 

201 

1917 

90 

1964 

90 

1921 

79 



1925 

37 



1934 

37 



1956 

59 




Ring width too small for application of the model. 
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Figure 4. Observed and reconstructed precipitation during 1971-2004. The calibration period 
lasts from 1971- 1987 and the verification period from 1988-2004. 



Véar 


Figure 5. Reconstructed precipitation for the full 225 years from 1782-2007. The grey line 
representa calculated Oct-Apr precipitation valúes (in mm) and the black line representa a 
smoothing filter of the precipitation valúes. The horizontal dashed line is the 225-year 
average calculated from annual precipitation valúes. 


El Año del Hambre, in conjunction 
with another drought lasting from 1798- 
1810 throughout northern and central 
México, is believed to have instigated hig- 
her food prices and food shortages which 
may have ultimately contributed to the 
Mexican War of Independence (Flores- 
cano and Swam, 1995; Cerano et al., 
2009). The longest drought during the past 
300 years in Chihuahua and Sonora took 
place between the years of 1 841 and 1 870 
(Díaz et al., 2002 and Cleveland et al., 
2003). Our reconstruction from Basasea- 


chi did show a decrease in precipitation 
during this general time period but the 
drought was most intense between the 
years of 1835 and 1851. One of the most 
severe droughts pointed out by this 
reconstruction occurred between 1909 
and 1911. This drought was feit by much 
of northern México and caused massive 
food shortages, leading to social discon- 
tent. It is for this reason that the drought 
has been hypothesized to be a major con- 
tributorto the Mexican Revolution (Cerano 
et al., 2009). 
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Table 5. Droughts lasting three years or more where the average precipitation for the 
entire time period did not exceed 86% of the 225 year average. 


3-4 



>4 


Year 

%of 

Year 


%of 

1788 - 

6 

1817 - 


5 

1892 

6 

1829 - 


6 

1902 

7 

1835 - 


7 

1907 

7 

1934 - 


6 

1943 - 

7 

1959 - 


6 

1954 - 

6 

1994 - 


6 

1971 - 

8 





The 1 950’s represent an entire 
decade of drought throughout northern 
and central México and the southwestern 
United States (Díaz et al., 2002 and Cle- 
veland et al., 2003 and Breshears et al., 
2005). Drought years in the Basaseachi 
area during this decade occurred in 1950- 
1951, 1954-1956 and in 1959 (Table 5). 
This drought is aiso known to have had a 
major social and economic impact throug- 
hout México (Cerano et al, 2009). 

Finally, the most recent drought 
highiighted by the chronology occurred 
between 1994 - 2002. The year 2000 was 
the driest year in the past 225 years. This 
drought was even greater in southeastern 
Chihuahua where it lasted through 2005 
and was one of the most severe droughts 
in the past 300 years (Cerano et al., 
2009). 


CONCLUSIONS 

Fine chronologies tend to be shorter and 
less well-correlated with climate than P. 
menziesii chronologies, but they are useful 
for reconstructing climate in areas where P. 
menziesii are sparse. A winter-spring preci- 
pitation reconstruction for bnp that exten- 


ded back to 1782 was developed. Climatic 
fluctuation was associated with the positive 
and negativa phases of El Niño/Southern 
Oscillation. This reconstruction identified 
periods of drought similar to those identified 
in the rest of northern México, but episodes 
of local climatic variability were aiso noted. 
Additonal collection in high elevation pine 
forests would be useful to expand the indi- 
vidual grid in México and contribute to a 
greater understanding of climate and its 
effect on dominant species in the respec- 
tive ecosystems. 
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